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Zur Berechnung der Strömung durch axiale Schaufelgitter 


Von Wolfgang-Hermann Isay in Berlin 


Die vorliegende Arbeit enthält Zahlenbeispiele zu der vom Verfasser im Jahre 1953 in dieser Zeitschrift ver- 
öffentlichten Theorie. Außerdem wird die T'heorie der damaligen Arbeit weiter ergänzt und ausgebaut, und es 
werden eine Reihe für die numerische Rechnung wichtiger Formeln nachgetragen, die damals weggelassen wurden, 
um den Umfang der Arbeit zu kürzen. 


The present paper contains numerical, examples to illustrate the author’s theory published im this Periodical 
in 1953. The theory of that paper is developed further. A number of formulae, omitted from the pervious paper 


‚for the sake of brevity but important for the numerical computation, are now supplied. 


L’article present contient des examples nume£riques au sujet de la theorie publiee par Vauteur en 1953 dans 
cette revue. En outre la theorie de l’ancien article est encore completee et Elargie, et une serie de formules impor- 
tantes pour la caleulation nume£rique est donnee comme supplement,. formules qui ont &t& omises ü ce temps-lä 
afın de reduire le volume de l’article. - 


IIlpennaraeman ma6oTa ComMep’kuT YMCJIOBbIe IIPHMepbI K TeOPHH, OIYÖOJUKOBAHHOÄ aB- 
ToPoM B 3TOM ;kypHane B 1953 Tr. KpoMe Toro, TeopuA YIOMAHYTOH Pa60TBbI MONONHAETCH, 
IPHMBONHTCH PA PopMyıl, BAKHBIX WIA IIPAKTHYECKUX BEIYNCHEHHÄ, KOTOPbIe PaHbINe ÖBLIM 
-ONYINeHEI IA COKPAINEHNMA 06’8Ma PpaooTkl. 


1. Einleitung 
In meiner früheren Arbeit „Beitrag zur Potentialströmung durch axiale Schaufelgitter“ 


ZAMM 33 (1953) S. 397 habe ich eine allgemeine Theorie zur Behandlung der Potentialströmung 
durch axiale Schaufelgitter entwickelt. [Jetzt als Arbeit I bezeichnet, und ihre Kenntnis sei 
im folgenden vorausgesetzt; (I, 8) bedeute Formel (8) aus Arbeit I, usw.]. 


Damals war es mir aus Mangel an geeigneten Hilfskräften nicht möglich, systematische 


numerische Rechnungen zu dieser Theorie durchzuführen. Dieses konnte jedoch inzwischen 
mit Hilfe der mir in der Deutschen Akademie der Wissenschaften zur ‚Verfügung stehenden 
Rechner nachgeholt werden.!) Dabei wurden auch noch neue, bisher unveröffentlichte theoretische 
Untersuchungen zum weiteren Ausbau der genannten Theorie rechnerisch erprobt. Diese letzteren 
Untersuchungen betrafen folgendes Problem: 


Die Strömung durch Einzelgitter ist ja streng eine stationäre Strömung; dagegen ist die 


Strömung durch eine aus Laufrad und Leitrad bestehende Gitterstufe eigentlich instationär. 
Der instationäre Charakter ist dadurch bedingt, daß das Strömungsfeld je nach der momentanen 
Stellung von Laufrad und Leitrad zueinander etwas unterschiedlich ist. Durch diese Unter- 
schiede ändert sich die Zirkulation um die Schaufeln zeitlich, und dadurch entstehen bekanntlich 
freie Wirbel, die mit der Strömung hinter der Gitterstufe abfließen. Die an den Leitradschaufeln 
entstehenden freien Wirbel müssen sogar noch durch das Laufrad hindurch. Diese Strömung 
ist in ihrer allgemeinsten Form sehr kompliziert, und ihre theoretische Behandlung erscheint 
fast aussichtslos. Nun ist jedoch der instationäre Anteil relativ unbedeutend, so daß die Strö- 
mung näherungsweise als quasistationäre Momentaufnahme ohne freie Wirbel behandelt werden 
kann, wie ich es auch in Arbeit I getan habe. 


Um die Arbeit I zu ergänzen, hat der Verfasser dann eine Näherungstheorie der instatio- 


nären Strömung durch axiale Gitterstufen mit linienhaften Profilen entwickelt. Auch diese ist 


bereits recht kompliziert und erfordert einen hohen Formel- und Rechenaufwand. Aus den bis- 
herigen Untersuchungen hat sich ergeben, daß die quasistationäre Theorie der Arbeit I auch 
bei Gitterstufen recht befriedigende numerische Ergebnisse liefert, während die instationäre 
Näherungstheorie mit größerem Aufwand schlechtere Resultate gibt. Dieser Umstand ist da- 
durch bedingt, daß bei der instationären Näherungstheorie eine Reihe sehr schematischer und 
wahrscheinlich nicht haltbarer Annahmen gemacht werden mußte, um das äußerst komplizierte 
Formelsystem in eine überhaupt der numerischen Rechnung zugängliche Gestalt zu bringen. 
Deshalb wird vorläufig bis zu einer weiteren Abklärung ihrer Brauchbarkeit auf eine Veröffent- 


lichung der instationären Theorie verzichtet. 


1) Die in dieser Arbeit mitgeteilter? Rechnungen haben Frau Märzke und Fräulein Schälicke unter meiner 


Anleitung und Aufsicht durchgeführt. 
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Für eine vereinfachte Behandlung des instationären Problemes bei Streckenprofilen unter 
Verwendung von Ergebnissen aus der instationären Tragflügeltheorie möge auf zwei Arbeiten 
von Kemp und Sears verwiesen werden.? 

In der hier vorliegenden Arbeit sollen Ergebnisse numerischer Rechnungen auf Grund der 
in Arbeit I entwickelten Theorie mitgeteilt werden. 

Außerdem werden wir zeigen, wie man den Einfluß der durch den an sich dreidimensionalen 
Charakter der Strömung (d.h. durch die Änderung der Schaufelzirkulation in radialer Richtung) be- 
dingtenfreien Wirbelhinterdem Schaufelgitter im Rahmen derebenen Theorie berücksichtigenkann. 


2. Einige ergänzende Formeln zu Arbeit I 
In Punkt 8 der Arbeit I habe ich den Weg zur Auflösung des sich bei der Behandlung der 
Gitterstufe mit linienhaften Profilen ergebenden Integralgleichungssystems 11,13) 


— (Q2+6) 42 
f 1 1 
at + Mesa +, | ade) Mater 50 dk, 
— (a,.+6+20,) wo  \ (1,15) 
1 — (a, +6) 1 Az n 
Ix(X) = =, IE Rn Hy, &) döı + = 0z(&5) I Bi + H39(&, 85) | dr 
mit > re 
e Ein r (2, &) | y, sın 2. (Y, —r Un) 1 1 y. 11 
Hp) In In DE > Hp %p)= > 1+ ya ? 
(Un ESCHE BEL a 1 er Pe au i 
1 
TE HIER ZIGE u 
A Sin T- a —$) + y, sin Yp— nd) 
H,&pE)=-7 5 5 or @+9) 
Ci =) 0 7 m) 


nur ganz knapp skizziert. Da diese Auflösung einen wichtigen Teil der vorliegenden Arbeit aus- 
macht, mögen die benötigten Formeln jetzt etwas ausführlicher mitgeteilt werden. 


Wir setzen 
= —- 10sh—y— ,—Ö; = —-0C0st,; | t ) 
0 2 Ye 
= —- 4,081 —— W—Ö; = —0C05%);, er 2) 
und ferner 
| 1 
f,&,) = sin 1,9? (Eo)R: 0,(&) = sin z, BR A FL RE en 


Die Funktionen sin af: cosvr u=1,23,..,.v»=0,1,2, ....) bilden in dem Quadrat 0 < 5 <ı 
T 


ein vollständiges Orthogonalsystem. Die stetigen Kernanteile H aus (1) können somit durch 
dieses Orthogonalsystem approximiert werden. Wir wollen uns gleich mit je sechs-gliedrigen 
Fourierpolynomen begnügen, diefür praktische Anwendungen ausreichen; sie werden numerisch 
mit Hilfe einer harmonischen Analyse in der üblichen Weise berechnet; man hat also 


6 6 
a,sint, HH, („m)= 2 26,2 sin nt, -cosyr,, (Ds 
u=1lv=0 
und dabei sind die b/P,® aus folgenden Formeln zu berechnen: 
BEE Ping; a 
De 9 e2 B asin tt, H(t,7,) cosvr, + 2. 5n t, [Ad 0)+— 1)” H(t,, ))) sin ut,; 


en) 


1 b 5 
= 5,2 PX sin t; Hl) +, Sin 4; [H(, 0) + Hd, m]|sinut,; 
pe ı >| hr: VE pr ARUHBE n8 ( ©). 
lin » 3). (—2) +5 sin 5, [H(t, 0) + Hd, a)]} sin ut; 
(= 1,22 50) 
=; wert (Schrittwei 
6 3 6 ırittweite 72/6) 


°) N.H.Kemp and W. R. Sears: Aerodynamic Interference Between Moving Blade Rows; Journ. Aeronau- 


tical Sciences 20 (1953) S. 585. 


T 2: <ı N 
N.H. Kemp and W.R.Sears: The unsteady Forces Due to Viscous Wakes in Turbomachines; Journ 


Aeronautical Sciences 22 (1955) $. 478. 
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Zur Vereinfachung haben wir in (5) die Indices p, q weggelassen. Das Integralgleichungssystem 
(I, 15) hat nunmehr ne: Gestalt: 


1 sin Da 1 
gı(t) — = oT) cost —cosh + 2 2 du sin Kb L, cosv Tı dt, 
ö Tee 


1f 6 6 
Er » | @g(T,) 2 Di bi; simwihcosr 1, dr, 
/ Fe 


v=0 


(6). 


TI 


1 IR 
926) = = | @,(T1) 2, PR: bi,’ sin ut, cos v 7, de; 
6 Ken 


1 i | sin t, Pi 

+ = @z(T5) er +2 3 Das. si etaCcosı2Tz dr 

Die Auflösung von (6) erfolgt wie stets nach der Schmeidlerschen Methode über die Kern- 

matrix. Wir verwenden für die Darstellung der Funktionen g9,(t,) und der gesuchten Lösungen 
@,(t,) die beiden vollständigen Orthonormalsysteme 


ap 1 
nl) = Vz inet, a ET er Be 


In unseren praktischen Beispielen haben die Funktionen 9,(t,) die Form 


.: V2 sopsi 
lt) =), Zr sinat,, VE I) ee el: 


wobei die Koeffizienten g/P’ durch eine harmonische Analyse mit der Schrittweite 2/6 zu berechnen 
sind. 

Das ae ee (6) ist dann dem folgenden linearen Gleichungssystem äqui- 
Yalent: (a ='1;". ; 


1 
de — F (wo be + wi Do) — Dr — — >> (w5 bi 2. O) ba 2: 


. (8). 
gr — v2 a ber u an De ehe eo a S (wi be 2,1 " E op br 5) 


Das System (8) für die zu berechnenden zwölf (meist sogar nur zehn) Fourier Koeiizienten Oo; 
@% ist in praktischen Fällen durch Itration sehr leicht auflösbar, da die Hauptdiagonale bei 
weitem überwiegt, wie sich gezeigt hat. Die gesuchten Lösungen des Integralgleichungssystems 


(6) erhalten wir somit in der übersichtlichen Form 


@,(T,) = = ur +), Zr oshn, (Ve er): 
Dabei sind die auch noch in allen w$? enthaltenen Konstanten »(P in der üblichen Weise durch die 
beiden Abflußbedingungen 
ol) = 0, or) = 0 
festgelegt. 

* Damit ist die Lösung (9) des Systems (6) endgültig und eindeutig bestimmt. Aus ihr erhält 
man durch die Transformation (2), (3) auch die Lösung des ursprünglichen Integralgleichungs- 
system (I, 15). 

Ferner ist 
Ni RT A el) 


Wir wollen nun noch eine bequeme Methode zur Berechnung der resultierenden Schaufel- 
kräfte angeben. 

Wenn man auf die Berechnung der Druckverteilung am Profil verzichtet, und sich nur für 
die resultierenden Kräfte interessiert, braucht man auch die etwas mühsame Bestimmung des 
Geschwindigkeitsfeldes am Schaufelprofil nicht durchzuführen. Wir ersetzen dann alle Schaufel- 
profile durch Punktwirbel der entsprechenden Gesamtzirkulation 7’, die wir in der Mitte der 
Schaufelkontur angeordnet denken. 

21* 
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m Kutta-Joukowskischen Satz steht die resultierende Kraft K senkrecht auf 
der en Geschwindigkeit w (relativ zu dem als ruhend gedachten Profil) und hat die | 
Größe K = o w T. Tist die Gesamtzirkulation des betrachteten Profils und w die Geschwindigkeit, 
die am Ort des Profils relativ zu diesem herrscht, abzüglich der von den gebundenen Wirbeln des 

i t induzierten Geschwindigkeit. 
Be Se Fall setzt sich ne resultierende Geschwindigkeit w beim Laufradprofil zu- 
sammen aus der Anströmung (u,, v,), der negativen Umfangsgeschwindigkeit (0, —V)) und den 
vom Leitradgitter stammenden Störgeschwindigkeiten (ur, Dr,) 3), und entsprechend beim Leit- 
radprofil aus der Anströmung (U, d,) und den vom Laufradgitter stammenden Störgeschwindig- 


keiten (ur, dr,). Dabei ist en 
ee ' 
sın . (ys(0) — N —Uu— ö)) j 


1 
7 


’ 
’ 


Sa + +9) 0057 WO) mm — 0) 


in (+0 +9) 


I; 
Dr = = 
ı > 27 
Ua + +) 0 RO) RN); 
11) 
2 
sın n (Yı ga ö) N7z(0)) 
Ur, = £ 
a T > 2 
Co (@ +0 +8) — 005 (m 9) 0); 
2 WORTE | 
Sin— (y+qa+6) 
: ls l 
u : 
el 2 DATE 
TOT (a + +0) — 08m (a — a8) — 14(0)) | 
In den meisten Fällen kann man ohne weiteres 
IR 
un, ®&0; er ur, & 0; Dr, 5] F ur 


annehmen. _ 

Es ist übrigens leicht einzusehen, daß die benachbarten Profile aus dem gleichen Gitter wie 
das gerade betrachtete Schaufelprofil, keine Störgeschwindigkeiten liefern, daz. B. beim Lauf- 
radgitter 


in © 1 
2 DL i y,(0) — in(0)—-ivl. N 
ist. 

Somit ist die resultierende Schaufelkraft zerlegt in ihre &- und y-Komponente an der Lauf- 


radschaufel (Rotor) gegeben durch 


KR= — o(w— V,+vr)T;,, Ka ziolug PUR) IE SR ORT) 
und analog an der Leitradschaufel (Stator): 
K9)=—o(v, +vr)T,; K9= ol, Fur) Wera) 


3. Zahlenbeispiel einer Gitterstufe 
Als Beispiel behandeln wir die folgende Gitterstufe (Bild 1): Es sei 
2 
a: 1 =11,795; 4 2, 1, Done 


und die Kontur der Leitradschaufel 
y,(&) = 0,8 — 0,9 (x, + 2,2) — 0,4 (x + 2,2). 
Ferner sei 


Yolz,) = c + 0,92%, + 0,42 


*) Alle Geschwindigkeitskomponenten u, v sind in der vorliegenden Arbeit auf das Absolutsystem bezogen. 


5) 


a a ni a 
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Laufrad 


Leitrad 


die Kontur der Laufradschaufel, wobei wir die 


Werte 


G= 0, Cı = 0,45, Cm = 0,90, cvı = 1.35 


untersuchen wollen, d. h. vier verschiedene Stellun- 
gen von Leitrad und Laufrad zueinander, die wir 
entsprechend mit Fall I, II, III, IV bezeichnen wer- 
den. Wir stützen uns nun auf das Formelsystem 


aus Pkt. 2). 


Es ist ohne weiteres klar, daß H,, und H, 
unabhängig von der Stellung von Laufrad und 
Leitrad zueinander sind, also in allen vier Fällen 
unverändert bleiben. 
Schaufelkonturen ergibt sich nach (1), (2) überdies 
noch, daß H,,(h, ı) und H;(t,, 75) formelmäßig ganz 


Bei 


den hier 


übereinstimmen, somit also auch nach (4) 


ist. 


(1,1) (2,2) 
De Ten De 


gewählten 


In den folgenden neun Tabellen ist das Resultat der Berechnung der bP,? nach den 
Formeln (5) für alle betrachteten vier Fälle zusammengestellt. 


Tabelle 1. Br == Be in allen vier Fällen 


Rt. 0 1 2 3 4 5 6 
u | 
1 +0,1268 | 0,6075 +0,06799 | +0,0077 | 0,0271 | +0,0022 | 0,0033 
2 —0,3532 | +0,1373 | —0,1192 | +00471 | +0,0273 | —0,0155 | —0,0010 
3 —0,0682 | -+0,0188 | -+0,0723 | 0,0900 | +0,0198 | -+0,0139 | 0,0062 
Er +0,0330 | —-0,1067 +0.0580 | -+0.0221 | —0,0297 | -+0,0015  -+0,0046 
5 —0,0059 +0,02866 | —0.0554 | -+0,0300 | -+0,00831 | —0,0043 | —0,0001 
6 0,00 0,00 0,00 0.00 0,00 0,00 0,00 
Tabelle 2. 52 im Fall I 
I 
EL 0 1 2 3 4 5 6 
it 
| | 
1 —1,7863. | —0,0635 | —0,0416 | —0,0198 | —0,0078 | —0,0038 | —0,0013 
2 40.0708 |, -+0,1240 +0.0835 +0,0409 | +0,0157 | -+0,0068 | +0,0023 
3 os | 01500 | -—omT7 | 0,0864 | —0,0219: | 0.0080, | 0,0027 
4 +0,0729 +0,1331 10.0985 +0.0536 | -+0,02022 | -+0,0070 | 0,0020 
5 0.0379 | —0,007 | 00551 | —0,0315 | —0,0120 | 0,0055 | —0,0010 
6 00 | .0,00 00 0,00 0,00 0.00 0,00 
Tabelle 3. 2,2 im Fall II 
Sn ) | 1 2 | 3 4 5 | 6 
a we.‘ | 
| | 0,0016 
—1,7516 | 0,0044 | 0,0052 | -—0,oo18° | +0,0024 | +0,0085 | +0, 
2 0.0110 | 0,0162 | 0,0066 | —0,0027 | —0,0059 | —0,0069 | —0,0031 
3 +0,0420 | +0,0695 | -+0,0391 Home | +Ho0ıs | +00007 | +0,0083 
0.0660 | —0,1125 | —0,0688 | —0,0312 | —0,015 ee 
5 1.0,0513 1.0.0886 +0.0560 | -+0.0258 | +0,0112 | -+0,0066 | +0,0027 
6 0,00 | 0,00 0,00 0,00 0,00 0,00 0.00 
Tabelle 4. Be im Fall III 
1 | Sum | soon | some | som | Home | som | aaa 
0.0477 | —0,0824 | —0,05 = N, 10, 
a ee 
== NN 852 — 0,0613 | ZN z Zen er ann» 
5 0.0221 = 0418 1.0,0326 40.0169 | -+0.0040 | —0,0003 | —0,0003 
6 "00 0,00 0,00 0,00 0,00 0.0 | .0,00 
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Tabelle 5. DR im Fall IV 
a 0 1 2 3 < a B 
u E 
— 0,0249 0,0156 +-0,0061 +0,0015 +-0,0003 
- ns 0 20301 Foo220 | 0.0001. | 0,0019. | 0,0002 
3 —0,0127 —0,0152 +.0,0043 +0,0125 0,0066 —.0,0007 
4 0,0392 0,0628 0,0293 +0,0043 — 0,0011 —-0,0009 En 
5 —0,0353 —0,0589 —0,0324 —0,0099 —0,0017 —0,0012 —, 
6 0,00 0,00 0,00 0,00 0,00 0,00 - 0,00 
Tabelle 6. er: im FallI 
a N) 1 2 3 4 5 6 
u 
1 +1,7797 —0,0596 +-.0,0582 —0,0477 +.0,0278 —0,0088 +-0,0006 
3 0,0378 —0,0767 0,0775 —0,0672 0,0432 —0,0186 +.0,0042 
3 +0,0254 —0,0524 +-0,0557 —0,0523 +0,0385 —0,0225 0,0078 
4 +0,0122 —0,0252 0,0273 — 0,0270 +0,0223 | —0,0163 0,0068 
5 +0,0055 —0,0112 +0,0114 —0,0107 +.0,0090 —0,0073 +-0,0033 
6 0,00 0,00 0,00 0,00 0,00 0,00 0,00 
Tabelle 7. De im Fall II 
2 0 1 2 3 4 5 6 
u 
ii 1,7598 —0,0157 0,0030 +0,0137 —0,0241 +.0,0250 —0,0118 
2 —+0,0171 —0,0289 0,0120 0,0120 —0,0295 +0,0345 —0,0171 
3 0,0147 —0,0261 +0,0145 —+0,0038 —0,0194 —+0,0264 —0,0139 
4 0,0053 —0,0093 +0,0047 +0,0032 —0,0104 ' —+0,0141 —0,0075 
5 0,0001 +-0,0005 —0,0020 0,0041 —0,0056 +0,0061 —0,0030 
6 0,00 0,00 0,00 0,00 0,00 0,00 0,00 
Tabelle 8. De im Fall III 
a | 
SE 0 1 2 s 4 5 6 
7 
1 +-1,7269 +.0,0466 —0,0456 +0,0355 —0,0156 —0,0035 +0,0056 
2 —0,0274 +0,0562 —0,0574 +0,0473 —0,0231 —0,0020 —+0,0062 
2) —0,0151 +.0,0320 —0,0354 0,0317 —0,0170 —0,0002 "0,0038 
4 —0,0048 +.0,0104 —0,0123 0,0112 —(,0053 —0,0020 |  +0,0026 
5 —0,0018 0,0039 —0,0039 +.0,0026 —+0,0001 —0,0028 | 0,0019 
6 0,00 0,00 0,00 0,00 0,00 0,00 0,00 
Tabelle9. De im Fall IV 
0 1 22 8 4 5 6 
n 
1 —+1,7342 +.0,0272 —0,0138 —0,0036 | —+0,0144 | —0,0153 +0,0069 
2 | —0,0263 0,0470 —0,0291 0,0042 —+0,0135 | —0,0182 +.0,0089 
3 —0,0238 0,0440 —0,0316 --0,0129 +0,0024 —0,0087 +0,0048 
4 —0,0117 0,0219 —0,0170 —+0,0093 —0,0028 +0,0002 0,0001 
5 —0,0030 0,0054 —0,0040 --0,0022 —0,0013 0,0015 —0,0009 
6 0,00 0,00 0,00 0,00 0,00 0,00 0,00 


ß Die Funktionen f(x), 12’) aus (1,15) bzw. g,(t), 9.(t,) aus (6) sind unabhängig von der 
Stellung von Laufrad und Leitrad zueinander. In allen vier Fällen unseres Beispiels haben wir 
somit die folgenden Koeffizienten g, nach (7): 


IT IT 
u 1% 2 +18 u); -— y® Bu: nd 0 ne 
2 IT 7 
m = Vz (2 Do 7 2V, Te 8 U.) gm = = 0,8 Ug; ga = 0 für & >= 3 : 


NEUERUNGEN ON FORTE > 
eV U ATER 


haben, daß Re 6 a = 


de en durch die: ei und ee 


Tabelle 10 


aD wi | 


+0,055 | -+0,013 | —0,022 | +1,000 
+0,108 | —0,032 | —0,002 | +1,131 

+0,096 | —0,018 | —0,014 | +1,064 

+0,073 | +0,012 | —0,032 |. +1,082 | 


Tabelle Il 
Fall wi | wi | wi wp wi» wi | ber 
iz Is +0,063 | +0,729 | 0,028 | +0,010 | +0,015 | —1,045 —1,852 
- DI :- | +0,076 +0,764 | —0,092 | -+0,004& | +0,016 | —1,086 | —1,924 
III +-0,063 +0,727 | —0,081 +0,007 +0,018 | —1,038 | —1,839 
IV —+-0,061 +0,778 | —0,080 | +0,010 | +0,018 | —1,114 | —1,974 
Außerdem ergab sich: h 
im Fall I: v, = — 0,369 us; Velen, 
im ,Kall [12 vo = — 0,387 ug; V. = — 1,860 u,, 
. im Fall III: vd. = — 0,406 us; V „= —.1,839 u,, 
im Fall IV: DV. = — 0,344 us; V.= — 1,777 u,. 


Man erkennt, daß die einzelnen Werte von v., Vo, /7, Z5 für die verschiedenen Stellungen 
von Laufrad und Leitrad zueinander gewisse Unterschiede aufweisen, die dadurch bedingt sind, 
daß die Strömung in der Gitterstufe nicht genau schaufelkongruent verläuft. Diese Unterschiede 
zeigen den an sich instationären Charakter der Strömung, doch sind sie nicht sehr bedeutend, so 
daß die Behandlung der Gitterströmung als quasistationäre Momentaufnahme gerechtfertigt ist 
(vgl. unsere Ausführungen in der Einleitung). 

Für die weitere Untersuchung der Gitterströmung, d.h. insbesondere die Berechnung der 
Schaufelkräfte müssen wir natürlich in allen vier Fällen eine bestimmte Umfangsgeschwindigkeit 
und einen Wert von v, zugrunde legen. Als solche wählen wir die arithmetischen Mittel der 
obigen Werte 

V,=— 1812 us; D, —— 0,377u,. 
Verzichtet man also jetzt auf eine genaue Erfüllung der Bedingung des stoßfreien Eintrittes und 
rechnet mit dem eben genannten Mittelwerten für V, und v,, so ergeben sich nunmehr die Lösun- 
gen &,(T,}) und w,(r,) mit den in Tabelle 12 und 13 enthaltenen Koeffizienten. Sie unterscheiden 
sich nicht wesentlich von denen aus Tabelle 10 und 11. 


Tabelle 12 
Fall | wi | wid | Me) | wi P w@ R get) 2 
I —0,019 | —0,717 | -+0,062 | +0,009 | —0,021 | +1,033 | +1,831 
II —0,135-) —0,772 | +0,106 | —0,039 | +0,003 | +1,110 | +1,967 
II —0,061 | —0,730 | +0,101 | —0,021 | -—-0,013 | +1,101 | -+1,951 
IV —0,071 | —0,780 | +0,073 | +0,009 | —0,030 | +1,051 | +1,863 
Tabelle 13 
Fall wi? | wg» j wo m ao iD ie 
I +0,093 | +0,736 | —0,077 | +0,010 | +0,014 | —1,013 | —1,795 
II +0,027 +0,752 | —0,094 | +0,003 | +0,015 | —1,142 | —2,024 
III +0,034 | +0,721 | —0,082 | +0,006 | +0,018 | —1,071 | —1,898 
IV | +0,092 +0,782 | —0,081 +0,011 +0,018 | —1,081 | —1,916 


Um nun mit diesen Werten von I, I’ do, v, die resultierenden Schaufelkräfte nach Formel 
(12), (13) berechnen zu können, fehlen uns noch die Störgeschwindigkeiten ur, Ur, Vr, Dr, aus 


ispi i ier Fä j die vereinfachten Formeln 
11). Bei unserem Beispiel kann man in allen vier Fällen ohne weiteres nfachte 
1 benutzen, wie man leicht erkennt. Dann ergibt die Kraftberechnung das in Tabelle 14 ent- 
haltene Resultat (in Einheiten von o u,). 
' Tabelle 14 


K®) K® 


R Y 
Fall K&) RN) arctg_—Y_ K(® 1 ) arc tg 5) 
% 


S 
K®) 


I +1,606 | +1,831 —+48,8° +1,661 _ —1,795 —47,2° 
II +1,851 | +1,967 +46,8° +1,796 | —2,024 —48,4° 
II +1,767 | +1,951 -+47,8° +1,692 | —1,898 —48,3° 
IV, | +1,696 | +1,863 +47,7° +1,755 | —1,916 —47,5 
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Man erkennt: Die Schwankungen der resultierenden Schaufelkräfte bei den verschiedenen 
Stellungen von Laufrad und Leitrad zueinander betragen etwa + 5% um die Mittelwerte 


K@® = + 1,730; K@ = + 1,903; K® = + 1,726; KR) — — 1,908; ? 
und bei der Kraftrichtung maximal + 1° um die Mittelwerte 
Ks) e K® % 
arc tg zo — + 47, 8°; arctg KR = — 47,8. 
T T 


Diese geringfügigen Unterschiede rechtfertigen den sehr großen Aufwand einer instationären 
Theorie nicht. 


4. Räumliche Ergänzungen zur ebenen Theorie 


Wir haben bisher die Strömung durch axiale Schaufelgitter als zweidimensionale Strömung 
in abgewickelten Zylinderschnitten behandelt. Es ist bekannt und auch unlängst wieder be- 
stätigt®), daß die ebene Theorie eine recht gute Näherung für die axiale Schaufelgitterströmung 
darstellt. Da zudem die Formeln der zweidimensionalen Theorie für alle Untersuchungen sehr 
übersichtlich und angenehm sind, liegt es nahe zu prüfen, ob man 
Effekte, die durch den an sich räumlichen Charakter der Gitter- 
strömung bedingt sind, im Rahmen der ebenen Theorie berück- 
sichtigen kann. 


So sind die Schaufeln in der ebenen Theorie als unendlich lang 
in radialer Richtung angenommen, was ja nicht der Fall ist; die 
Schaufelzirkulation ist im allgemeinen in radialer Richtung ver- 
änderlich, und dadurch entstehen hinter einem Schaufelkranz 
freie Wirbel, ähnlich wie hinter einem Propeller, wenn auch von 
geringerem Einfluß. Diese freien Wirbel induzieren an den Schau- 
feln eine zusätzliche Geschwindigkeit. Wenn man grundsätzlich 
die ebene Theorie der abgewickelten Zylinderschnitte beibehalten 
will (was wir tun), so interessiert von dieser zusätzlichen Geschwin- 
digkeit die Axial- und die Umfangskomponente. Hierfür werden 
wir eine einfache Näherungsdarstellung ableiten. 

Wir betrachten einen axialen Schaufelkranz mit n Profilen, 
die wir durch tragende Stabwirbel der Zirkulation /‘(r) ersetzt denken. (Bild 2). Die x-Achse 
unseres Koordinantensystems falle in die Richtung der Schaufelkranzachse. Es sei 


z2=Tc0so; $ =r0 008.9; y=rTsing; n=osiny. 
(R = : = R,) R, ist der Außenradius desSchaufelkranzes, R;ist der Radius der Nabe.5) In jedem 
Punkt 9 der gebundenen Stabwirbel wird dauernd ein freier Wirbel der Stärke — = do indu- 
ziert, dessen Wirbelachse annähernd in die durch do 
cosß-i + sinßcosy-j— sinßsiny-f 


angegebene Richtung fällt‘), und der dann mit der Strömung abfließt. Für die Berechnung der 
von den freien Wirbel induzierten Geschwindigkeit ersetzen wir die n rotierenden Stabwirbel 


J „A ckeret: Nabeneffekte bei gebundenen radialen Wirbeln; Vortrag GAMM-Tagung Berlin 1955. 

) 516,0, Y sind mit ©, y,2,r,gp an sich gleichbedeutend und dienen in der üblichen Weise als Integra- 
een hler — eg ist natürlich jetzt keine Dichte! i,j,f sind die Einheitsvektoren der kartesischen Koordi- 
naten. 


6) tg entspricht in der Bezeichnungsweise der ebenen Theorie d i i j i 
ee N g S e dem Anstiegswinkel in der Mitte der 


Kai 2u8,r u, 
N 


# 


DE sc EZB ed. EEE an Zu DE inc 
VRR TAN. es 
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n 


durch eine äquivalente tragende Stabwirbeldichte der Stärke 5 
0 


_Ie) über den ganzen Um- 


fang des Schaufelkranzes. Dann ist der gesamte Raum hinter dem Schaufelkranz kontinuierlich 


mit freien Wirbelelementen der Stärke — 3 = do ausgefüllt, deren Achsen die oben an- 
0 


gegebene Richtung haben. Die von diesen freien Wirbeln induzierte Geschwindigkeit ergibt 


sich nach dem Biot-Savartschen Gesetz zu :”) 
(naar Ex 
2 dydo. 
\ ri ). ds 


(S) 


g=uyti+vy.j+tw t=— 


Ra 

ER 
4 

e=R, v 
Dabei ist 

d3 = ds (cosß -i + sinßcosy-j—sinßsiny-H =de(i + tg ßcosy-j—tgßsiny-d 
der elementare Richtungsvektor der Achsen der freien Wirbel, und 

t=(@ — Hi+(rsing—osiny)-j+ (rcosp— ocosy)-E. 


Für die weiteren Untersuchungen benötigen wir die Geschwindigkeit c, bei der hier ver- 
wendeten Näherung nur in der Schaufelkranzebene, können also x = 0 setzen. Außerdem kann 
wegen der Zylindersymmetrie des Problemes ohne Einschränkung der Allgemeinheit nochg = 0 
gesetzt werden. Die beiden uns interessierenden Geschwindigkeitskomponenten sind dann gegeben 


durch 
Ra 2n © r ? 
1 | | | nd (r cos y — o) tg 
Aadı rd —_— -dE-dy-do . . . (14), 
er v=0 80 2ndo Ye@+ r2+ —2rocosy)® 
und 
R, 2a © er 
# er. | n al geosyor+ätgßsiny ae a ee 
Eu > EN, ae Ve@+r+ 0°—2rocosy)® 


Die Integration von (15) ist ohne weiteres möglich, wobei wir auf eine Reihe bekannter Integral- 


u se F 
formeln verweisen.®) Man erhält, wenn noch durch die Relation er = l die Gitterteilung / 
eingeführt wird, endgültig 2 


1 Ar 7 
ur FR] ee) rl 


?) Das Vorzeichen im Biot-Savartschen Gesetz ist so gewählt, daß einem im Uhrzeigersinn drehenden 
Wirbel eine positive Zirkulation zukommt, wie auch in Arbeit I. 


8) Es ist 
[0,0] 
de 1 (18) 
ee ee a de nr ee 0 On en od > 
s Ve +r +0 —2rocosy) Pre Tg Son 
oo 
& JE 1 
RE ee 2 = ea = a nn Det ur, aa este ee ae (19), 
; Ve +1? +0?—2rocos w)® Vr: +02 —2rocosy 
27 DEzrun or 
r — 0 C08 y 5 5 
Sr — — l ) = 32. 77 Fr re 20 
| r? to —2rocosy Ba r fürrsor en 
( 
27 2; 
“ 0 —rcosy a, e>r 
| - as Sr sen dıy —— N a (21), 
Pe 2n 0008 y 0 für o<ı 
2 an : 
| sinydy °—  _ 0: | 1 ‚sinyd iR An (22) 
ö r? +02 —2rocosy ; ö Vr: +0 —2ro cos y 


Zum Beweis bemerken wir: Die Formeln (18), (19) ergeben sich ohne weiteres nach den Regeln der elementaren 
Integralrechnung. Formel (22) gilt, da die Integranden ungerade mit 2 x periodische Funktionen von y sind. 
Die Formeln (20), (21) erhält man mit Hilfe der Residuenmethode durch die Substitution z = eiv und Integration 
über den Einheitskreis der komplexen z-Ebene. 
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i i i im Zylinderschnitt r = 
16) stellt die y- oder Umfangskomponente der von den freien Wirbeln im Zyli 
en Geschwindigkeit dar. Aus (14) erhält man nach Integration über & und y: 


Ra 
n 1 dl, 


Näherungsweise kann man wegen der bekannten Verwindung der Turbinenschaufeln in radialer 
Richtung (s. Abb. 3) annehmen, daß 
/ ’ 0, 
igB. u Wo) ae 
8 8 5 ” [1 

d. h. unabhängig vom Radius 0 bzw. r ist. Damit kann die letzte Integration sofort ausgeführt 
werden, und wir erhalten endgültig ai 

| = Try. 
u,(r) u ») \r) [/r) ( a)] ( ) 
(17) stellt die x- oder Axialkomponente der von den freien Wirbeln im Zylinderschnitt r = konst. 
induzierten Geschwindigkeit dar. 


Die in Punkt 2) der Arbeit I entwickelte ebene Theorie zur Berechnung der Strömung durch 
ein axiales Schaufelgitter kann jetzt leicht vervollständigt werden. 
Die Strömungsbedingung lautet: Längs der Profilkontur y(x, r) muß 
Don) — Vor) + 2m) + Dym y@); 7) 
uo(e) + ul) + u, (x, y@), r) 
sein °), und zwar jeweils für den betrachteten Zylinderschnitt r = konst. Für jeden dieser Zylin- 
derschnitte ist also eine Integralgleichung der Form aufzulösen: 


2 ul) ya) LO — TRIER 1) + 204) —2 Vol) 


DU) = WR, 1) = 


210.,.2&%) 
\(r) Urieh \ 


[ven 4 + H@8&n dE 


N 
BUN 


+4 


mit 
a 
Ir) = | o(&.r)dE. 

Die Auflösungstheorie von (23) bleibt gegenüber der in Punkt 6 der Arbeit I angegebenen 
unverändert. Wegen der auftretenden Werte /(R,) und I'(R,) muß (23) natürlich zunächst für 
diese Zylinderschnitte (außen und an der Nabe) gelöst werden, bevor andere Zwischenwerte 
r = konst. behandelt werden. 

Für die Berechnung der resultierenden Schaufelkräfte in den einzelnen Zylinderschnitten 
unseres Schaufelgitters kann die in Punkt 2) entwickelte Methode mit Formel (12), (13) sinngemäß 
angewendet werden. Man hat also jetzt 


= ee 

u. Wr Re el ar (24) 

K,= Kart y‘(0, r) y'(0, r) er . 
vr Urn IK) In) —-75 IK) I(R,)| I(r) 


5. Zahlenbeispiele eines Einzelgitters unter Berücksichtigung der räumlichen Ergänzungen 
zur ebenen Theorie 


Wir betrachten ein axiales Schaufelgitter mit einem Radienverhältnis R =0,6 und 
wollen die drei Zylinderschnitter = R; r—= R,undr = 0,8 R, untersuchen. (Beispiel 18) 
Die Schaufelkontur im Zylinderschnitt r = 0,8 R, sei y = 0,422 + 0,9 x; (also wie diejenige 


der Laufradschaufel bei unserer Gitterstufe in Pkt. 3). Die Schaufelkontur am Außenradius 
r= R, sei!) s. Abb. 3) 


022-1 
und die an der Nabe r = R, sei y, = 0,5892 20.5 
°) Das gilt für ein Laufradgitter; bei einem Leitradgitter ist V, wegzulassen. 


10) Um den Gedankengang hier nicht zu unterbrechen, sollen die Überlegungen, auf Grund derer die Kon- 
turkurven %, und y; gefunden wurden, in Pkt. 6) am Schluß der Arbeit gebracht werden. 


PEOBEE VRR SERUEAEEN: — 


- 
3 


TERN 


na ı 
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7 Ferner sei für alle Zylinderschnitte a = 1, und iz, konstant unabhängig vom Radius, außer- 
em ; 
ale) Er 


0,8 R, Ku ee 


Per 
Ay, 
Für den Kern der Integralgleichung (23) können wir im Zylinderschnitt r = 0,8 R, die b,, aus 


Tabelle 1 beibehalten; dagegen ergab die Berechnung nach Formel (5) für die be(r= R,) und 
bu,(r = R;) das in den Tabellen 15 und 16 enthaltene Ergebnis: 


Tabelle 15. 0 (Beispiel 1) 
u ed per 2 Be u 1 Er EFT PET a rn re .. 


v 
0 1 2 3 
ie 4 5 6 
1 0,0937 +0,1604 0,0704 +0,0736 —0,0347 —0,0043 -+0,0054 
2 —0,1806 +0,1407 —0,3162 +0,0493 +0,0757 —0,0161 —0,0094 
3 —0,0460 +0,2130 +0,0832 —0,1725 +0,0214 —+0,0175 —0,0008 
4 +0,0317 —0,1148 —+0,1438 —+0,0249 —0,0336 —0,0018 —0,0030 
6) —0,0078 —0,0292 —0,0732 +0,0450 —0,0120 —+0,0184 —0,0028 
6 0,00 0,00 0,00 0,00 0,00 0,00 0,00 
Tabelle 16. 0 (Beispiel 
v 
0 il 2 3 4 5 6 
En 
1 +-.0,0322 —+1,1169 —+0,0118 —0,0526 —0,0060 +0,0109 —+0,0001 
2 —0,6816 +0,0236 -+0,0747 +0,0191 —0,0164 —0,0067 0,0078 
3 —0,0245 —0,1555 —-0,0280 —0,0216 —0,0070 —-.0,0020 —0,0035 
4 +.0,0676 —0,0271 —0,0266 +0,0082 —0,0150 —+0,0023 —+.0,0007 
5 —0,0037 +.0,0661 —0,0201 —0,0073 0,0060 —0,0058 —+0,0006 
6 0,00 0,00 0,00 0,00 0,00 0,00 0,00 


Wir lösen nunmehr die Integralgleichung (23), d. h. das ihr äquivalente Gleichungssystem 
(1,25) [Sonderfall von (8)] für die drei betrachteten Zylinderschnitte. Mit den angegebenen 
Zahlenwerten wird 


9% 75 175 
re 9, = 20,(0,8 R,) — 2 V (0,8 R,) — 1,8 u, — (0,9)? is (og — a) — je (0 —w); 
= 1 
jr 9 = 0,8 u, + 0,4 0,9 — (w, — o®); ET RR 
7 Yr 


Er Read" Ve: 
v2 g® = 2 v.(R;) — 2 Vo(R,) — 0,3 us— (0,15)? 3 (vH — wi); 


2: DIA MN, 
"= 9% = 1,066 u, + 0,533 : 0,15 nz (N a): N 
= 1,75 4 | 
y2 ga— 2 v,(R,) — 2 VolR,) — 3,06 u, — — 5 (wi — 8); 


er ga— 0,64 us; Bd, =, 
7 
Dabei haben wir gleich mit (10) d. h. hier 
Ve hn)=!Yros: IXRB)=Vro®9; 
die /' durch die w, ersetzt. 


Die Berechnung der Koeffizienten &,, 5, o@®(ß =1,... 5) macht wegen des Überwiegens 
der Hauptdiagonale in den betreffenden Gleichungssystemen keine Schwierigkeiten. 

Die Koeffizienten &,, ©, &® wurden aus den Abflußbedingungen o(a) = 0; wa) = 0; 
@,(r) = 0 bestimmt. Ferner haben wir für alle Zylinderschnitte stoßfreien Eintritt zugrunde 
gelegt, d.h. den Wert von 

v.(0,8 R,) — V (0,8 R,) vo(R,) — Vo(R) vDo(Ru) — VolRa) 


Ug Ug Ug 


IR) = iz a® 
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aus den Bedingungen (0) = 0; ®;(0) = 0; @,(0) = 0 ermittelt. Das Ergebnis dieser Rechnung 
zeigt Tabelle 17. Dort sind die Koeffizienten @5, @p', @$” der Lösungen @(r), ®;(T), @.(T) gemäß 
Formel (9) angegeben. Und zwar alle in Einheiten von u,. Ganz rechts stehen die nach (10) aus- 
gerechneten /-Werte ebenfalls in Einheiten von u,. 


Tabelle 17 (Beispiel 1) 


r=(0,8h, Ya — R; r= fer 

ß @B or wg 

0 | —1,092 | —1,148 | —1,075 Ru 
1 | +0.065 | +0025 | +0,05 7108 R.) 1,937; 
2 | +0,765 | -+0,746 | +0,789 et : 
3 | —0,085 | —0,029  —0,080 I(R;) 2,035; 
4 | -+0,008 | -+-0,066 | —0,029 Be. 
5 | +0,018 | +0,004 | -+0,025 IR.) = 


Außerdem ergab sich: 
fürr=R;: 1, — V, = + 0,156 u,; für r= 0,8 R,: 1, — Vo = + 0,926 us; 


ı 


tür Br, Vo 1,5640... 


Mit den erhaltenen Werte können nun noch die resultierenden Schaufelkräfte nach (24) er- 
mittelt werden. Man erhält 


ER 
fürr=R: KR,-+0380ow, K,=-2000u; areıe 7a = -— 81,1°, 
ZT 
R K} R 
ier=0sH.R,= 1, ou, RK, 1200 are Kae 41,85 
. - K, 2 
fün = DW uK, = #2, 94.0 ER, = 00 ars 5 = — 33,1°. 


T 


Um ein Urteil über den Einfluß der freien Wirbel zu erhalten, rechnen wir nun das gleiche 
Beispiel noch einmal ohne deren Berücksichtigung durch; d. h. wir behandeln alle drei Zylinder- 
schnitte unabhängig voneinander nach der rein ebenen Theorie. 

An die Stelle von dem in Tabelle 17 enthaltenen Ergebnis tritt dann bei sonst gleichen Be- 
dingungen das folgende aus Tabelle 18: 


Tabelle 18 (Beispiel 1) 


’ Ir=08R, reR;, eg 
ß @B | op | wg 
0 —1,102 — 1,156 | —1,075 } en r2, 
1 14 4+0.067%)%4-0,038. 10.086 Da 
2. 10850,771 1.5.0, 7810=:°0:289 ar HE 
3°) 0.085.200 0.0 
4 | +0,008 | +0,066 | —0,029 LE e 
5-1: +0,018;.| 50,004 00a wu 
Außerdem ergab sich: 
für r= R:9, —V,= +0,1570,; fürr= 0,8 R,: u — Vo, = + 0,906 us; 
fürr=R,: u - V,=r154u, 
Für die Schaufelkräfte erhält man: 
R 4 Red 3 K 
fürr=R: K,=+0320u; K,=—2049 ou}; areig 7a = —. 81.17, 
T 
fürr=0,8R:K, = 1,769 o 2; K 1 ar =, 7 
=0,8R,: RK, = + 1,76% o 13; \; = 1,953 2 us; Farce 27 —= — 47,8, 
2 
fürr—=R: K= 2.992 0 uR: 7 2 ı K, o 
ER .,.=+2220wW K,--- 19050; "arete Kr 3» 


T 
Ein Vergleich der in Tabelle 17 und 18 enthaltenen Resultate für die Zirkulation und die 
sich aus ihnen ergebenden Komponenten der Schaufelkräfte zeigt: 


1. Sa 


DErTATER 


VRR ERENTO DE 
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Die Unterschiede zwischen der 
Rechnung unter Berücksichtigung der 
freien Wirbel (räumliche Ergänzung) und 
der Rechnung nach der rein ebenen 
Theorie, d.h. also der Einfluß der freien 
Wirbel, sind hier bei diesem Beispiel 
ohne Bedeutung. Der Grund hierfür ist, 
daß wir (wie in Pkt. 6 gezeigt wird) die 
Innen- und Außenprofilkontur y, und y, 
nach den Formeln der technischen Tur- 
binentheorie unter Verwendung des Flä- 
chensatzes, d.h. eben konstanter Schau- 
felzirkulation in radialer Richtung, kon- 
struiert haben. Die Rechnung nach der 
richtigen Schaufelgittertheorie liefert 
nun zwar eine geringfügige Abhängig- 
keit der Schaufelzirkulation I/' vom 
Radius r, die darauf beruht, daß die 
Gitterströmung nicht genau schaufel- 
kongruent verläuft. Jedoch ist dieser 
Effekt nicht groß genug, um freie Wir- 
bel von merklichem- Einfluß hervorzu- 
bringen. 


Zylinderschnitt r=R; 


Zylinderschnitt P-Q8R, 


Zylinderschnittr=R; 


Bild 3. Schaufelkonturen in den drei betrachteten Zylinderschnitten unserer 
Beispiele aus Pkt. 5 
(ausgezogene Kurve: Beispiel 1, gestrichelte Kurve: Beispiel 2) 


Als nächstes (Beispiel 2) betrachten wir ein axiales Schaufelgitter, das in allen übrigen 
Daten mit dem vorigen übereinstimmt, nur ist jetzt die Außenschaufelkontur y, = 1,3 x + 0,352 
und die Innenschaufelkontur (R, = 0,6 R,) y; = 0,4% + 0,52, d.h. willkürlich angenommen 
und nicht wie bei Beispiel 1 unter Verwendung des Flächensatzes konstruiert, 

Dann erhält man für die 5) (r = R,) und die b\, (r = R,) das in den Tabellen 19 und 20 


enthaltene Ergebnis 


Tabelle 19. a (Beispiel 2) 


Er 0 | 2 3 4 5 6 
BT | | 
1 +0,1154 +0,2319 +0,0730 +0,0512 —0,0296 —0,0023 —0,0040 
2 —0,1821 +0,1455 —0,2284 +0,0426 -+0,0460 —0,0125 —0,0042 
3 —0,0633 +0,1496 +0,0674 —0,1090 +0,0124 +0,0089 +0,0013 
4 +0,0195 | —0,1076 | 0,0881 +0,0173 —0,0155 —0,0009 —0,0028 
5 +0,0040 | —0,0111 | —0,0435 +0,0214 —0,0086 | -+0,0106 —0,0019 
6 0,00 0,00 0,00 0,00 0,00 0,00 0,00 
Tabelle 20. DD (Beispiel 2) 
es - 
E 0 il 2 3 4 > 6 
| 
1 +0,0794 +1,1307 +0,0350 —0,0506 —0,0160 +0,0104 | —+0,0003 
2 0,6534 +.0,0700 .0,0438 0,0441 —0,0076 | 0,0169 | -1.0,0076 
3 —0,0549 | —0,1491 +0,0642 —0,0373 +0,0202 +0,0073 | —0,0099 
4 +0,0594 —0,0726 —0,0062 +0,0234 —0,0264 | 0,0056 | 0,0034 
5 —0,0102 +0,0695 —0,0520 -+0,0022 +0,0154 ‚| —0,0106 +0,0018 
6 0,00 0,00 0,00 0,00 0,00 | 0,00 0,00 
Ferner ist jetzt 
2. | 1,75 1.75 
y2 9, = 2v,(0,8R,) — 2 V,(0,8 R,) — 1,8 u, — (0,9)? — (@, — ©) — —— (@ — N); 
Tre r VYr 
2 1,75 
y2 9, = 0,8 u, + 0,4 0,9 x Vo; 9-9 = =0. 
10 7 


V- ga» = 2v,(R)— 2 Vo(R,) — 0,8 u, 


„1,75 4 
= (0, 4) Yr 3 


(v9 — 0W); 
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% DO, ’ RE 
2 = 1,0 u, + 0,5-0,4 nn 3 (eo); u ee VE 
7 
2 1,75 4 Re 
VE so=20.(R9 2 VuR) — 26 — Vz (a: 
v2 3 = 0,70 us; u 


Der weitere Gang der Rechnung ist genau wie vorher (s. S. 331/32). Als Ergebnis erhalten wir: 


Tabelle 21. (Beispiel 2) 


ß r=(0,8R, r=R; Rn 

@B a |. WB 
9%. 1,114 Ua ar: 
1 | +0,068 | -+0,055 +0,087 T(0,8R.) 1,975, 
2 | +0,780 | +0,733 | +0,82 Br ue 
3 | —0,086 | —0,064 | —0,089 I(R,) 1,969 , 
4 | +0,008 | +0,052 | —0,016 rer 
5 | 40018 | +0,009 | +0,02 Ra) 2,02 


Außerdem ergab sich: 
fürr=R;: v,— V, = 0,418 u,; fürr= 08 R,:v, — V, = + 0,95 us; 
für = Ra = u > 122u,% 


Die resultierenden Schaufelkräfte ergeben sich nach (24) zu: 


K o 
Arne: K,= + 0,824 o us; K,= — 1,984 0 u; arc 7 = — 67,5°, 
T 
K 0 
für DNS. Re 1,8108: K,= — 1,997 o us; arctg = — 47,8, 
T 
K 
fürspenr K,= + 23,634 0 u; K,=— 2,0210 u; arc tg = — 37,5°. 


T 
Auch bei dem zweiten Beispiel behandeln wir nun noch alle drei Zylinderschnitte unab- 
hängig voneinander nach der rein ebenen Theorie, um im Vergleich mit den Resultaten aus Tabelle 
21 ein Urteil über den Einfluß der freien Wirbel zu erhalten. Das Ergebnis dieser Rechnung 
enthält Tabelle 22: 


Tabelle 22. (Beispiel 2) 


2. r=0,8R,| r=R, | BE: 
OR RR OR 
DE NE 00a ER) = — 1,953, 
en ee nes 
5 Tools 10.008 on FR) = — 3,021. 
Außerdem ergab sich: 
fürr=R: w—V,= +045u,; fürr=0,8R,: u — Vo = + 0,906 us; 


fürr=R,: v,— V,= +1,303 u,. 
Und die Schaufelkräfte: 


tim ’R.: K,= +0,8110 u}; K,=-—-19%3 013; arc tg = = — 67,5, 
T 

fürr=0,8R,: K,= + 1,769 0 u2; R,=—1l930U; arctg = = — 47,8, 
T 

fürn: K,„= + 2,633 0 u2; K,= — 2,021 0.3; arc tg = = — 37,5% 
T 


Man erkennt, daß auch bei Beispiel 2 


Mar ‚ bei dem der Flächensatz der technischen Turbinen- 
theorie nicht erfüllt ist, 


die abgehenden freien Wirbel kaum einen Einfluß ausüben. 


a 


ER 


BRENNEN NT NONNEN N DELIONGRRU VERUTTENN EN TNORONSTE EUIRRENIETEHRTETTEN 


ei en für die der lächensätz der technischen urBinentisone nicht erfüllt 
irkt der Effekt der Winkelübertreibung der Entstehung freier Wirbel hinter dem Gitter ent- 
gegen. Der Einfluß der letzteren dürfte nur r bei ganz unverwundenen Schaufeln eine größere 
Br Rolle EBEN 


6. Zusatzbetrachtung zur Verwindung der Schaufeln in radialer Richtung 
Zum Schluß wollen wir noch kurz angeben, wie die Profilkonturen y, und y„ des Beispiel 1 


° aus Punkt5) gefunden wurden. Dazu sind wir von der Voraussetzung ausgegangen, daß das 
 Mittelprofil r = 0,80 R, zum Laufrad der in Punkt 3) behandelten Turbinenstufe gehöre. Rechnet 


man nach den technischen Formeln der schaufelkongruenten Strömung, so sind für den Zylinder 
‚schnitt r = 0,8 R, mit den Schaufelkonturen auch die Richtungen 

tg wı(r); tg w.(r); iger);  tgcr) 
der Ein- und Austrittsgeschwindigkeiten 

wm = {wu V}; mein VE; 0-05 9={uD} 

bekannt. (Abb. 4). Außerdem gelte der Flächensatz | EEE: mg 

i Zzp (eo) =.const;z „rp,(r), — const; 

und es sei u konstant über den ganzen Radius. ") 


Dann läßt sich in elementarer Weise die Richtung 2 
der Eintrittsrelativgeschwindigkeit an der Nabe durch 
die bekannten Größen aus dem mittleren Zylinderschnitt r 0, 


darstellen, nämlich 
R, R, NS 
m) (3 7) + tem). sth 


Eine genau analoge Formel gilt auch für tg iv,(R, 
(Vertauschung von , durch ,) und für den Außenradius R, 
wobei R, durch R, zu ersetzen ist. 


Bei der Turbinenstufe in Punkt 3 war 
nl. tem) Ol; tele, = el,t, ts, () = —Ol. 


i R e ! $ 
Ferner ist bei Beispiel 1 aus Punkt 5) = — 10.79: = — 1,25. Damit ergibt sich 


SS 


Bild 4 


y 1) =tgmwi(R) = — 0,916; tg wı(R,) = y;(l) = + 1,217; 
y, = 1) >= tg w,(R,) =+ 0,89; tg 1,(R,) = y„(l) =+ 2,17 . 
Beschränkt man sich nun auf Profilkonturen der Form Ax + Bx?, so erhält man sofort die 
in Punkt 5) angegebenen Kurven y, und y,. Mit der bekannten technischen Formel für die Zir- 


kulation 
IT,=—u:1(tgw, —tgi;) 


bestätigt man sogleich, daß I, unabhängig vom Zylinderschnitt ist, wie es nach dem Flächen- 
satz sein muß. 

u) Vgl. z. B.: C. Pfleiderer, Dampfturbinen, Wolfenbüttel 1949, 8. 137. 
Eingegangen am 21. September 1956. 


A physical system is suposed to exist in one of a discrete set of mubually ‚xchus 
en time from one state to another are allowed. The “evolution” of the sy 
studied. S Rear re Er 
N ö 6 & n . Fe Fe 6 antile E 
On suppose quun systeme physique existe dans un &tat faisant partie d’une quantile 
celuant en et que des transitions accidentelles d’un &tat ü Vautre sont admises dans 


L’ „evolution“ du systeme et ses characteristiques ergodiques sont Etudies. An id ‘ 

Ilpennosaraerca, YTO HEKOTOpan Pusmyeckan CHCTeMA HAXONUTCH B OMHOM { 
HOTO MHOKECTBA UCKIOYAWIIMX ApyT Apyra CocTosHNÄ, HPHUEM B TeueHme BpeMeHH p: 
INAIOTCH CHYYAÜHLIE EPEXOoNEI 13 ONHOTO COCTOAHNA B pyroe. Vecnenyerch pasBuTue 1 
CHCTEM N e& IPToNnyeckMe CBOÜCTBAa. 


| 1. Introduetion ur 
rn m. A rigorous treatment of the so called ergodic problem in classical dynamics is comprehen- Rn 
y N sible only to those trained to the discipline of mathematical analysis. The ordinary physiit 


approach, the ergodic properties of some simple stochastic systems can be studied in physical 
terms and the essential difference between stochastic and deterministic systems explained 
within the frame work of elementary probability theory. The author is of opinion that the 
phenomenological approach is helpful to the physicist who is more concerned with the under- 
standing of the phenomenon and does not want to be troubled or diverted by too many analytical 
difficulties. An attempt has been made, however, to translate these analytical difficulties into 

- “physical” language so that the physicists may also be aware of the nature of the mathematical 
problems that arise from a rigorous treatment. 


2. Description of a simple stochastie system 


Let us consider a system which is capable of assuming one of a finite number of states 
5» 89 ...,S, and evolves in a Markovian manner with respect to a one dimensional 
parameter, Att=t, we assume that the system is in the state S,.. We define the function 
II(Qli; 1, t,) as the probability that the system is in a state S, at i, given that it was in S, at 
h(& >). If we further assume the process to be homogeneous in t, IT is only a function of 
l, — t, and so without loss of generality we can set, =0,1,—, = tand write II(jli; t) instead 
of IIjli; tt) where ,—t, =t. The evolution of the process as i tends to infinity can be 
studied if we know that the limit /Z(jli;t) as t—0 exists. For the sake of simplicity we shall 
assume that 


nn A, ee 


Ai; d—RG)A+A) if ji 
1— A DR(KI) — 0A ker Fee Eee 
k,k=+i : 


and 


as 
re >08 


Phenomenologically, this means that if a system is in a state S;, at/=0, in an interval A it 
will “jump” to another state with probability proportional to A or will continue in the same 
state 5, with probability 1— A 3 R(kli). 
k, ki 

Once we assume the existence of this limit, the ‘forward’ differential equation of the process 
is obtained as 
oll’jli;d) ir 
= Mili;d DR) + 3 Mais d AGO. 

k,ks+j k, kaj 


*) Read at the GAMM Conference in May 1956 at Stuttgart, W. Germany. 
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where IT is a column vector whose j-th member is IZ(jli; D)!), I the unit matrix and R the matrix 
with nondiagonal elements R(k|l) and the diagonal elements 1— S R(k|l). In the above 


'equation the initial condition at £ = 0 is defined by BB 
IG;0)=0, he) vr 
Hd; 1 in 


For the validity of equation (2), this restriction on the initial conditions is not necessary. We 


can assume an arbitrary “injection spectrum” at?=0i.e., the system is in state S; at!= 0 
with probability ®(i). In that case we must define //(j; f) as the probability that the system is 


in jatt, and note that 
IMS DT Er 
[2 
IIG ; 2) satisfies the same equation but at 2=(0 we must remember to impose the condition 
II(|0) = ®G). Br 


We note that since R(j|jk) A is a probability magnitude R(j|k) > 0 for all j, k, +k). 
Also if R(j|k) > 0 for any particucular j, k say j* and k* then IZ(jt|kt; 2) > 0 for all finite 


t>0. For the purpose of our discussion we shall assume one of the two following restrietions 


on the non-diagonal terms of the matrix R. 

1) R(ilk)>0 for allj, k, (j=+k). In this case we say that the states are completely 
connected in an infinitesimal transformation. In physical terms, this means that a 
transition of the system from every state to every other state is possible in an infinitesimal 
interval of time A, with probability proportional to A. It is of course trivial to note that 
II(jli;t) > 0 for all j, i for finite t>0. 

2) R(j|k) = 0 for some j’s and k’s but for every such pair of states, say S;« and S», it is 
possible to find a sequence of states S,, S%,,... , S7, such that each of the following terms 

Ru|R*), Rd), .., Rlolinı), RG*|l,) 
is greater than zero. In such a case we shall show that /Z(j*|k*; 2) >0 for all finite {>0 and 
hence I/(j|k; i) > 0 for all j, k, for all finite £>0. In this case we say that the states are 
completely connected for finite transformations though not for infinitesimal 
transformations. 

The proof of the above property of // is immediate if we adopt the following device. We 
divide the finite interval (0, fZ)inton +1 finite intervals (0, i,), (kı, fo)». - . » (n-ı, £,), (,, &) and 


' use the Markovian property of the process to write 


HIG*\k*;) = 3% IKU|R*; t) IR; a —t) alon;  —) A Inst) >. (6) 
where /,%,,...,/,„ represent intermediate states and $ denotes summation as each of the 


l;’s goes over the whole ‘gamut’ of states. Since I, l,,...,/„ represents only one such sequence 
of intermediate states and each of the terms on the R.H. S. is positive or zero 

II(j*|k*; 0) > II(U |R*; 4) li; a —t) Il nn — if) A) N. 
Each of the terms on the R.H. S. is greater than zero since its corresponding R term is greater 
than zero. Hence I/(j*|k*;) >Oi.e Ali; >0O for all j, k, forali>0O. 

Without bothering with the various definitions of stationarity we shall confine our dis- 
cussion to the stationary nature of the distribution function /T. If I/(jli; it) tends to the same 
limit IZ(j) as t—> oo for all values of i, we shall call this distribution, the stationary distribution 
of the process. Physical examples can be given to demonstrate that the process may not have 
a stationary character phenomenologically even ifthe distribution is stationary— a feature which 
is more striking when we consider the total number of states m to be enumerably infinite. 
However, at present, we shall understand stationarity to be a mathematical property associated 
with the Z/ function. 

A simple consequence of the existence of //(j) is that if in (5) we assume &(i) = I/(i) then 
II; t) = II(j) and independent of t. This result is obvious if we conceive of II) ati=0 to 
be the stationary distribution of a process started at — oo. This result will be used in the course 
of our study. 

* We have now an adequate description of the process and its analytical features to formulate 
the “ergodie problem” for such a process. | 
> 1) In some papers on stochastic theory IZ(j|i; t) is called a transition probability. The author is of opinion 


that this notation is unsatisfactory from a phenomenological point of view since the state j might be entered 
prior to t. However R(j/i) : 4 can be called a transition probability for the infinitesimal interval A, or R(jli) 


the transition probability per unit time. 


22 


PR 


that ati=0, the system is in $,. It is well known tl 
can oceur in a finite interval of time and a typical realis 


Er 


We Shall consider a particular rea 


on 


n«ı 
known 0) I 
on of the proces 


follows. = ee 
The system continues in S, from t=0 to! = Li at i, it Jjum 
that state till {= 1, when it jumps to $,, and so on..., jumps to 
in that state till {. We can ascribe a measure to a class of such ‘realis 
transitions occur in the intervals BE 
td), Gate dids... lat din) 


exp [- {ei)h + ek M+ te. 
Rilio) Rliz|i) - - - Rlinlin) dt, dt, .. . di, 3 
where nr Q\ 
ei) = 3 RG). ea 
We must remember that the above equation does not forbid the recurrence of any state. VRR 
In ihe above expression any i, can be identical with any i,. Hence we can speak of the “first 
realised state” S;,, and “second realised state’”’ S;, and so on. (An n-th realisation is also termed "m 
an (n— 1)th recurrence.) From the above expression it is possible in principle to obtain the 
probability P(r, t; jjk) dr that the system spends in S,a ‘time’ of magnitude betweenrand 7 + dr 
in the interval (0, £) given that at 1= 0 the system was in $S,. We must remember that 7, the 
time spent in state S, is the sum of the times spent in the first, the second, ..., and n-th realised 
state j. (n can be made as large as we please.) - 
The computation of this probability is obviously a very formidable task. It however 
Gl; 
t 


respectively as 


en, 


seems interresting to define a stochastic variable x(jli; t) = as the fraction of the time 


spent in state S, in the interval (0, f) given that the system was in ,at!=0. ( i;) is the 
stochastic variable representing the time spent in state j and its probability frequency function 
is of course P(r, t; jli). It is clear that the study of x is more informative than that of r as 
! — oo since the moments of x may tend to a limit while those ofr may not. A comprehensive 
study of x(jli; {) should deal with its frequency function and so we define P(z, t; jli) dx as the 
probability that x + de > x(jli;t) > x and study the function P(x, 1; jli) as ! tends to infinity. 1 
It is quite possible that in some physically interesting cases this tends to a limit and it is to 
such cases we turn our attention. This is precisely the study of the ergotic properties of process. - 


4. Ergodie property of the process 


Since it is more convenient to deal with the moments of x than with its probability fre- 
quency function, we shall state the ergodic property of the system in the following form, 


If as >00 li; > LH for. all ran 
then 

etz(li5d} ID an. u Sr 

e {2770} — [ZZ] - > 2 Re De re 


where e denotes the expectation value. This is equivalent to stating that the stochastic variable 
x(jli; t) “converges” in probability to //(j). Note that & {x@/li;d} and /I(li; t) are not related 
in a simple manner for finite {. 

Before attempting to prove this, we first realise that it is a long and tedious task to obtain 
the moments of x from its frequency function. On the other hand, we adopt the expedient of 
expressing the moments of x in terms of the J7 function. ‘To obtain the first moment we shall 
use the simple and self-evident theorem that the mean of the sum of random variables (corre- 
lated or uncorrelated) is equal to the sum of their means, in its most general form to include 
the case of a continuous infinity of random variables. To obtain higher moments, we apply a 
generalised correlation theorem explieitly stated and established by the author recently, (Rama- 
krishnan, 1954) but probably known earlier by implication in the works of others. 

We shall divide the interval (0, 2) at los... , in into N intervals A,, As, ...,Ar- 
By making N very large, each of the A’s can be made as small as we please. If we consider a 
state 5, each interval A, is either spent entirely in a state S, or not at all. This 
follows from the property of the stochastie system already referred to; in an infinitesimal inter- 
val A the probability of transition is O(A) while the probability of continuance is 1—0O(A). Thus 
to each interval A, =, —l;_ı we can assign a stochastie variate st) = rd; th) — TS) 


result _ 
4 sul er 
a or wartet; >= Lett) = Fett; n} = Jaounar 2.2... 
2 j b 
% singe we are more interested in sch: f) the Gacion of the u spent we note that 


ta) =4 | mucnar. re Ba ta ide 


Be: Ari 


To obtain the higher moments of r we Sette the correlation theorem referred to before 
in the following form = 


e elng; = ll. Je rl; i)örsh)-- RE, SEEN 


- Note that in (16) the variables of integration 4, t,,...,t, need not be ordered; they are varied 
over the whole range of integration. The right hand side is a formal expression and an exact 
expression can be found only if we know the nature of correlation between the various ör(j; t;). 
This is possible if we assume Z,,t,,...,t, are ordered and make use of the Markovian property 
of the process and equation (11). Eu <t,:---...<d, then 

ee förl;t) kit)... rt); = gli) Ali)... ls — tn) di di...d, (17). 
In (13), integration with respect to each of the variables is to be done over the whole range 0 
to t. Therefore we can also rewrite (13) in the form 


e{"Gl;d;-n! far fat. .. ‚Soli 1) Ali; —t)... Gi, — tn) dt, (18). 


ori 


_ 


5. Proof of the ergodie property 


The proof of (16) now consists of four steps. 
E- 1) Let Fl; dt, (FÜ respresent the probability that the state S, is entered for 
2 the first time between ? and ! + di given that the state at {= 0 was $,. We shall prove that 


[Fülse) ar 1 a5 | >09 -Bonlall 7 


ö 
provided that R(j|k) > 0 for all j, k i.e. the states are completely connected for unfinitesimal 
transformations. This condition is sufficient but not necessary. It will be shown later that 
the necessary condition is that the states should be completely connected for finite transfor- 


_ mations. 
E Defining K(ilj; ) as the probability that the state jis not entered at all between O and f, 


; we note that _— 
m... OLE a 3 1 20113207 Zee} 


We shall show that K-—0 as t— oo without to compute K, by adopting the familiar 


mathematical device of “majorisation’. 
For this purpose we build up a “comparison process” in which there are only two states A 


and B and we set 

NE oe RR a ih 
where the right hand side refers to the original process and R(jle)is the least R for the given]. 
Note that A corresponds to and B corresponds to the entire aggregate except j which for 


convenience we shall call “non j”. In this comparison system the probability of a jump from 
22% 


7 


\ 


1 äh ü a ana a a ie ne ah PRLPLATTEN ' 
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a state non-A (i. e., B) to A is less than the jump from non- to j in the original system. We 
immediately see that K(A|B; t) the probability that A has not been entered at all in this 
system is greater than K(jli; 2) for all values of t. Though this is “intuitively obvious”, for 


the sake of completeness we shall ““demonstrate” its truth by the following argument. We, 


define for formal purposes, a function R(j|non j) A which represents the probability that given 
the system is in a non-j state it jumps to the j state in A. Of course to express R(j|non j) in 
terms of the R(j]k) we must know: the “distribution of the system in the aggregate of non ] 
states”. However, dividing the interval (0, £) into N intervals A,, A,,..., Ay and making 
each of them infinitesimal by letting N tend to infinity, we can write K(jli; i) formally in 
the form of an infinite product. 

KGli;) = [11 — RG|i) A] II — RÜjlnon j) Az] 11 — RGlmon j) A,] . . » 22). 
The R(j|non j) terms on the R. H. S. are not the same but are merely formal expressions. In 
the “comparison system” since non A=B 

K(A|B; i) = [1 — R(A|B) A,] [L — R(A|B) A] "11 — R(AlB)A,]. - -- » 23). 
Thus K(A|B;t) >exp [— R(AlB){] as N—o0. We observe that R(j|nonj) > R(A|B). 
And so since each of the terms in the infinite product representing K(jli; i) is less than the 
corresponding term in the expression for K(A|B; t) we obtain 

Kil;d)< K(AlBED << exp RAD) IE vr (24). 

Since we have assumed R(A|B) = R(jle)>0, K(jli;tl) >—0 ast—oo. Hence 


t 
[FGli; öde —1. as) >00. Sr (3). 


0 
2) The second step consists in proving that e {x(jli; {)} — L(j) independent of the initial 
state iast—oo. This is equivalent to stating that 


efaGli;d} and etaGli;d)} 
tend to the same limit. We first write the well known stochastic integral equation 


gli) = | Foli; ey gli; — 0) ar N I 


which is self-explanatory from the definitions of // and F. We can now write 
t t v 
2 | 4 Mr 
e{z(li;D} 3 IIGli; tt’) dt [a [raisı INS E —IN)dEs 2202 
Ö f) () 
We split the integration with respect to f’ on the R.H. S. into two parts 


1f% t 
ai +) tor IIrcd tg a. (28). 


Inverting the order of integration in each of the above integrals, w i 

| ‚we can w 

equation (27) as a sum of three terms, 5 ne 
Jı ’ Js D Jg . 


to t, 


1 |. 2 
= | Folseyar [ Bone dt 
u v 


where 


2 


1 SEN 7 j lt 
2-7 | Foist ) dt [av —enar 0 
to; 


t t 
1 L ee 1) 2 4 
nen r] Fli;t’) dt | Adi; — r)ar 
% v2 
By taking t, suitably lar tti 
gt ge and letting 2—oo we can make J, and J 
a 
near and simultaneously make J, aproach the product of the ee ; 


t 
ii F Ole, N 
lim J td, en ee (30) 


ero arbitrarily 


t 
18: 
a Ft Il): f ! 
im | GlisE) di en ne ee re 
0 


t> 00 


VEREINE 


LAU 
x 


are 


2) 
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 Remembering that the first limit is unity, and the second is lim & {xz(jlj;D}, we have 
t>o0 


Im net 2; im ET Le (32) 


i. e., the limit is independent of the initial FE &. 


3) The third step is to adopt the following simple device taking came of the result 
obtained above. Since. we have proved that Am e{z(jli; {)} is independent of i we are entitled 


to write 
Im. 220: Di 5 Dii).limie LElflid) Verer (33) 
—00 v t>00 

provided that p2 D(i) =1 i.e. we are at liberty to choose any “injection” spectrum &(i) for 


the probability distribution in various states at {= 0. So, in particular, we choose the “sta- 
tionary’’ distribution N and write 


als DI 20) Im ’etadit, u ll] BE 1 [mo 2) dla 38% 


Inverting the order of summation and integration, we have 
t 


al nn BROT EdL 2. a 0) 
Remembering i 
Selm Z2los LE De EHEN SO) 
we have 
LU) =lme{zGli;d)} = lm 70 AR LT m 
t> 00 t>00 


4) The final step consists in computing e{a"”(li;t)} as t—oo. To do this we make 
use of (18) and prove that, 
i in—1 


lim a [aa [mawo mann. HG; m — ni) de, 


rot 
n—1 


t 
= Im z [Zul !)dt' In mann Fr A TE 
0 


Since each of the limits on the R.H.S. has been 4 to tend to the same limit //(j), the dis- 
cussion on the ergodie property is complete if we prove (38). 

To accomplish this?), we make use of a lemma, which is very similar to a classical theorem 
on integral resultants (e. g. Theorem 42, G. H. Hardy’s Divergent series, Oxford). 


Lemma 
IE ©(t') and W(t’) are non-negative and integrable over any finite interval (0, ?) where 
t>0andifk>0 
t 


e D(t) dt is bounded for > 0 and tends to Mas L—x 
ö 
s 
z Y(t)dt' is bounded for £>0 and tends to Nas L—oo 
M) 
then 
t et 
im Gr | Par [arm 2. 


t>o0 rl 


To avoid any interruption to our discussion the proof of this lemma due to C. T. Rajagopal 
and Padmavally is given in the form of an appendix. 


al am grateful to C. T. Rajagopal and Miss K, Padmavally who helped me to prove (38). 
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“ The application of this lemma to our problem is immediate. 


We first establish (38) in the special case when n = 2 i.e. when 
i 


t 
| Zi; Aus) En (40). 


Interchanging the order of integration we have 
t 


wi | dh | Mas ag — wa 


ö 


{ ; 
= [Miisnau | main — ak 
Ö m 
t N 
= [Polsa | maysnar. 0 
0 


In our lemma, take : ” 

; k-1, =), 7O=Udli;d. 

The conditions on ® and Y are obviously satisfied by // since it is a probability magnitude. 

Noting that M = N = L(i) we have | > % 
im {2,49 = Le [UND Fe 
t>00 


To prove (38) for general n, after changing the order of integration, we write (18) in the form 


t 
! en 
tel) = [Mol 
() 


en 1 15 |» r »|.* Mi [2 .|. t ’ 
x n di, ae Ex [at HGl; t,) Id; ar Il; An | ee RE (43). 
0 0 0 
Writing for convenience the therm in the square brackets on the right hand side as /„ _ ı(E — I) 
& ” 
we note that /„_ı(f) corresponds to f D(t') dt’ and n to k in our lemma, and so assuming the 
0 


existence of the limit 
t 


| 1) El re 


0 


1 


m 
which we identify to be 
lim & {z(\j; ) Pr 
t>o°' 


N 
we have 
n! 
. Ri een (il; der RE EN AA 
lim e {Gl D) a Bel 2 (44) 
In a like manner we have 
himetau(i dr Li) 0a e{ x(lj; 917 RR MR 
t>o0 t>o0 > 
Thus we obtain 
lim 8 {2*(/\1,9,: = [LO =. II) en (46). 
t>00 


2 
One important corollary can be obtained from the result | F(j 


;t)d' —1ast—oo; 
i.e. every state is entered “sooner or later” with certainty. This immediately implies that 
every state is entered at least N times “sooner or later” where N can be made as 
large as we please. This result is obvious if we divide the interval (0, i) into N parts and 
let each part tend to infinity independently. Since in each part as it tends to infinity 
every state must occur at least once, it must occur at least N times in the aggregate. In other 
words an event which is certain to happen will occur infinitely often as too! This can be 
called the “Fundamental periodie property” of the stochastic system. The demonstration 


ii 
. 


a 
u re ee ee 
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t 
of the ergodie property rests on proving that [FUENar>1 as {—oo. For this purpose 
0 


we assumed that the discrete set of states was completely connected for infinitesimal trans- 
formations. This assumption is sufficient but not necessary. The neccessary condition is that 
the set of states must be connected for finite transformations. This is easily shown by an argu- 
ment that is now familiar. ; 

If the system is connected for finite transformations we know that //(jli; t) > 0 for all 
J,i for all finite £. So we divide the interval (0, 2) into N intervals (0, 4), (4; bo)». (nd) 
andgexpress I1(1jE, Din terms of. FL(Rli; 4), INBIRSE, 1)... NGe 15: 6 lade 
scheme I/(l; |; 1; ,— t;_ı) correspond to the R’s of the system which was completely connected 
for infinitesimal transformations. Hence we apply the ““majorisation” to these /]’s and demonstrate 


t 
that Ki; ) = 1— [ F(jli;t’) dt’ is less than the product of N terms each of which is less 


0 
than unity. As N tends to infinity the product vanishes. We must take care that as N tends 
to infinity each of the intervals remains finite in length. Note therefore that the infinite product 
cannot be expressed as an exponential. 

The above work was done during my stay at-the Yukawa Hall, Kyoto, Japan, in the 
spring of 1956 at the kind invitation of Professor H. Yukawa. I wish to express my thanks 
to him for his hospitality and to the Asia Foundation for its generous financial support. To 
Professor Bartlett, I owe the suggestion to prove (38) for general n. My discussion with 
Dr. T. E. Harris in Los Angeles encouraged me to publish this work which is a “simple minded’”’ 
but I hope useful, approach to what is normally regarded as a difficult mathematieal problem. 
Dr. Harrisin a private communication has pointed out that the results of this paper are implied 
in Birkoff’s proof of the ergodic theorem if suitably interpreted for the stochastic process 
considered here. 


Appendix 
To prove the lemma we first set 
1A 
od) = [Dlt)ar. Re el ee er rd area et) 
0 
so that (39) can be written in the form 
t 
lim ID —-P)E=MN TR AIR AT), 
t>0 


0 
Given any small positive numbers 6, e by the hypothesis on ® we can find Z, such that 
(MH -fE <oßlt—!)<(MTgHlt—t)% for I-!>ä>ö,.. (A). 
Also by the hypothesis on &(t), (tl) we can find constants M,, N, such that for all2> 0 
t t 


®] St) dt’ <M,, = U ERN eg a eat. (AA): 
Ö 
Now 
Ä k 1 fa-9t ! 
ET ve) D(—r)di =" f veiyoßd—)a+ f Pie) D(t—t) dt |= 1,+ 1, (say) 
je+1 e+1 5 ao) 
ö (A—5) 
where using (A—4) we find that 
= en MEN JE=(k+IM) MN. 2.2... . (A6). 
Further by (A—3) As 
k al 1 { 7 r ’ 
<(M-+e) en Pr) (E— V’)* dt 
hi Kl ei 
| >(M=3) a | pr) (te VE de’ 
l 0 
while 
rn k+1 
k n NK AH If! Ned —- 5 = 
en Pe) (red | DO ar]... 


Ö 0 
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where t 
+ (naar <EFL r ee 
as) 


Using (A—8) and (A—9) in (A—7) we get 
5 - t 
k 1 ß ’ 1 ’ 
tz a tyat | 
I \ = A 
Al t ug 
| > [7 (E — Ped! —(M — e) (k +1) N. 6% 
\ B L ” ” .ı. 
Now our assumption of (c, 1)-summability of Y(f) to N implies (ce, k + 1) summability of 7) 
toN,i. e. 
DT 


t 
In Sn | Per N are ne ee 
i> 


0 
Using (A—6) and (A—10) in (A—8), then letting {—o0o and using the above relation, we see 
that (A—5) leads to the lemma in the form (A—1) since ö and e are arbitrary. 


Remarks on the lemma: 


The existence of limits M and N regarding ®(f) and (tl) concerns large as well as small 
(positive) values of {. Hence the condition of boundedness at 1— +0 is also imposed. The 
assumption of ®(f) and W(f) being positive is used in obtaining the inequalities of the above proof. 
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Eingegangen am 24. September 1956. 


Abschätzung der Lösung der ersten Platten-Randwertaufgabe 
nach der Methode von Maple-Synge 
Von Rüdiger Nicolovius in Hamburg?) 


Unter Verwendung der von Synge eingeführten Hyperkreismethode wird eine Abschätzung für die Lösung 
der ersten Platten- Randwertaufgabe hergeleitet, die der von Maple für die erste Randwertaufgabe der Potential- 
theorie angegebenen Abschätzung entspricht. Die Methode wird an einem Zahlenbeispiel erprobt. 


Using the hypercircle method introduced by Syngean estimate is obtained for the solution of the 
first biharmicon boundary value problem which corresponds to Maple’s estimate for the solution of the first 
boundary value problem in potential theory. T’he method is applied to a numerical example. 

En appliquant la methode d’hypercercle introduite par Synge une evaluation pour la solution du premier 
probleme marginal de la plaque est derivee correspondant & l’ Evaluation pour la premiere probleme marginale 
de la theorie potentielle indiquee par Maple. La methode est &prouvee au moyen d’un exemple num£rique. 


BpiBonuTca IpMm HOMONM TAMIepKPyToB CHH3KA OMeHKa Is pemmennst ITepBoli KpaeBoü 
3alayy B TeOpUM INNACTUH, AHaIOTUUHan ONeHKe peimenns IIepBofi kpaeBoi sanaum) ns 
ypasHuenna Jlarsıaca, ykasannof Mannem. MeTon WINOCTPMPyeTcH IPHMepoM. 

1. Definition der Aufgabe 
In einem zusammenhängenden, beschränkten, offenen Gebiet B der euklidischen x-y-Ebene 
mit dem stückweise glatten Rand I’ genüge die Funktion w den Bedingungen 


AM er) Ina, 


w= [($) 
on Auer ee a ne de 
A 965) 


Dabei bedeuten s die Bogenlänge auf I’ und d/On die Ableitung nach der äußeren Normale. 
A ist der Laplace-Operator: 
2 02 0% 0: 0% 
= an PR als / = —_ 2 = 
en 


!) Institut für Angewandte Mathematik der Universität Hamburg. 


BEIcH Er ei 


2. Durchführung r Hyper 


an das Verfahren von Synge [5] er nd dient zur \ 
nitt. Es sei F der lineare, reelle Funktionenraum, der von allen r" 

wird, die in B definiert und dort zweimal stückweise stetig partiell diferen erbar 
d. Definiert man als Eee u 


(9 y) = 2/4 Ap: Aydedy = (pp), @, Bo) = Koll (2), 


| so sind Assoziativ- und Distributivgesetze erfüllt, jedoch IM die Metrik nur positiv semidefinit, 
3 da jede Potentialfunktion einem Vektor aus F mit der Länge Null entspricht. Die Gültigkeit der 
für das folgende wichtigen Schwarzschen Ungleichung ist aber auch im semidefiniten Falle ganz 
 elementar beweisbar. , 

Durch die Forderungen Bu 


Ali DB 
Dee 


ER | WW auf /'; v zweimal stetig differenzierbar in B 

on 
werden im Funktionenraum F zwei Teilmengen gekennzeichnet, deren Durchschnitt die Lösung 
w von (1) ist. Mit Hilfe des Greenschen Satzes 


n a ae a P) 
Wa 2, ee a vr Av)ds.. . A 


berechnet man nach (1) und (3) 


- 


ee 


F“ wu =, | wm Arnd + 


4 eben ee a N  ı 


on 
Diese Beziehung läßt sich en als ges einer Hyperkugel in F schreiben: 
lv—z|?= R2 mit 2 = (u +2) »undn Ri= = kun elek 


E Um die Schranken des Ergebnisses (22) möglichst eng werden zu lassen, ist es nötig, u und 
2 v so zu bestimmen, daß ||u — v|| klein wird. Wählt man zunächst beliebige u= u, undv =D, 
E die (3) erfüllen, und dazu Funktionen u, und v,, die den homogenen Bedingungen 


Miro "mB ER 
0) 
z me — c ? B a Re lee 7 
ed 03H auf 7 re Be (7) 


v, zweimal stetig differenzierbar in B 
genügen, so illen die Funktionen 
uv-Ww+au und v-ew+bVU,..... Le)! 
die Bedingungen (3) für jede, Wahl der Konstanten a, und b,. Über auftretende In- 


dizes ist hierbei und im folgenden zu summieren. Sind die Funktionen u, und v, orthonormiert, 


ist also 
(u,ı,) =6;,. (dd) = Öj1 a ea Sn 


{Nach (4) ist überdies immer (1,, v,) = 0}, so wird 
ur |w vv) +, — bull... 22022... 
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zum Minimum für 
d; = — (u;, Ug ng Do) ’ b, = (v; Ug en D,) ne @ E Er (10) 


und nimmt dann den Wert 
M N 
4R = |u—ov|P = u, vl? — 2 (u. u, — do? — 20 es N 
il = 


AR Damit ist die Spitze des Vektors w € F auf einer verkleinerten Hyperkugel mit dem Mittel- 
punkt (vgl. 6) 


2 Lu +0) = 5 In + de + |. 3 


festgelegt. Nach (4) ist ferner 
(w — z, u) = (w— v,u)=0 ((=1,...,m),] ie, 
(w — 2, v,) = (w— u„D,) =0 Gel..n) j 


d. h. die Spitze von w liegt außer auf der Hyperkugel (6/12) auch aufm+n Hyperebenen, die 
durch die Spitze von z gehen und orthogonal zu den u, bzw. v, sind. w liegt also auf einem Hyper- 
kreis des Funktionenraumes F [5]. Dieser Hyperkreis läßt sich auch in der Form 


Dez Rd. iger ae gu)=@v)=0...(d9) 


darstellen. Danach gilt für einen beliebigen, aber bekannten Vektor p aus F 


w—2,)=-RGM)=-RGQp) mit 9-9 -Ww)w— wg); . . . (AD), 
also nach der Schwarzschen Ungleichung 


WEI -R@AISRIFl mit Tl = el — wm’ — 2 Wo. - (16). 


3. Aufstellung von Näherungswert und Schranke 
In Analogie zu den Betrachtungen, die Maple [3] für die erste Randwertaufgabe der 
Potentialtheorie durchgeführt hat, wird die ‚dritte Greensche Identität‘ herangezogen. Ist 
P ein fester Punkt von B (im Inneren von B + I!) mit den Koordinaten x, und y,, und ist 
w,» der Wert von w in diesem Punkt, so gilt bekanntlich 


i oh 0) | ( 0 ow 
= | | = F 
WW» = h AAw dx dy - | = Auen ar Aw \ds - w an Ah ar Ah) ds (17), 
B r 2 
worin h die Elementarlösung oder Grundfunktion zum Punkte P bedeutet 2), nämlich 
2 
== 25 -Inr mit = (t— Hp? +lYy—y)® ....% MB). 


Berechnet man das Skalarprodukt 2(h, w), so ist der Greensche Satz (4) nur dann anwendbar, 
wenn zuvor ein kleiner Kreis K und P ausgeschnitten wird, weil Ah = = (1 + Inr) in diesem 
7 


Punkt singulär ist. Jedoch kann man nachträglich den Radius og dieses Kreises gegen Null 
streben lassen, so daß gilt: 


2(h, w) = lim hp dx dy — lim Ki Aw—h > Aw )as -H & Aw — > Aw) ds (19. 
0 00 5 on on Y on on 
Dabei wird /”, der Rand des Kreises K, in mathematisch positivem Sinne durchlaufen. Man 
überzeugt sich leicht davon, daß das zweite Integral in der Grenze verschwindet. Das erste 
Integral bietet keine Schwierigkeit, da h in P stetig ist, und so gilt die formal direkt zu bildende 
Gleichung 


N 9] 
2(k,w) = hp dx dy - Kia Aw— h 2 Aw) ds. Se 2 
B T 


Setzt man dieses und (1) in (17) ein, so erhält man 


2 (h,w) = w — | I 
Tr 


0 
5; An gAn)as Be Re 


?) Die verwandte Methode von Diaz und Greenberg [2] benutzt hier die Greensche Bunkkan statt Ah. 


En RN 
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Abschätzung nach (16) mit 9 = A liefert das Ergebnis 


Balz gAn) ds — (hu + u) < Ju —ol- | 


mit IA = IrlP — 2. 9° — 3 @, Mn 
= = 


N 


N — 


Zwei Dinge sind an diesem Resultat bemerkenswert: Während bei den Verfahren von Ritz 
und Trefftz vom Abstand auf die Elemente von F geschlossen wird, was eine positiv definite 
Metrik erfordert, treten in (22) außer w, nur Skalarprodukte und das von der Wahl der Ansatz- 


funktionen unabhängige | (1 Ar g An) ds auf. Daher ist hier nur eine positiv semi- 
; I; 


definite Metrik erforderlich. 


Der Vergleich mit der Arbeit von Maple [3] zeigt ferner, daß hier ||h|]®? konvergiert, im 
Gegensatz zur ersten Randwertaufgabe der Potentialtheorie. Es ist also nicht nötig, wie bei 
Maple, einen Kreis um P auszuschneiden. 


4. Beispiel 


Es sei die Aufgabe gestellt, die Auslenkung w, des Mittelpunktes einer allseitig einge- 
spannten Platte (siehe Bild) unter konstanter Flächenlast in Schranken einzuschließen. Die 
Bedingungen (1) lauten hier 


AAw=]1 Ina eh, 


Di E 

BU 2, art | unaret23 
w_) Au J «] ( ) r 
a | el=1, ylsl 


Die Ansatzfunktionen stehen in Tabelle1. Sie erfüllen die Be- 
dingungen 


AAus=ıls AAu,=0 inB (aber Au, #0), el 


CR : - 
ar le Bu auf ’(w,=0), == 1,243 


SS? 


Dabei bedeutet das ° über den u, und v,, daß diese Funktionen noch nicht orthonormiert sind. 
Außerdem enthält Tabelle 1 die Grundfunktion A (18) für den Ursprung. 


Tabelle 1. Ansatzfunktionen 


Funktion {17} h Ap 
ae 212 | Done 2 
U = at) rg (+?) 
= + y 
a, =—y | 12 —y?) 
rer —5iryt + y 20 (# — 622 y? + y) 
vo, = (2 + y—.2) (2? — 1)? | 6038 + 72 0 y? + 12.22 yt — 196 & — 120 2? y® — 4 y! + 17622 + 
1 32 y2 — 32 
d, = (ea +2 — 2 (ae — 122 | 94224 12420 y? + 30.04 yt — 352 2 — 264 ty? — 2422 y1 441024 
+ 132 22% + 2% — 1522 — 87 +8 
= er? — 22 (a 1 | 22848008 y2 + 90 at yt + 122248 — 1226 — 276 2% y2 — 156 22 y! — 
| — 46 + 260% + 2762? y? + 50 yt — 240? — 72 y? +8 
IR 2 as N B => R Aa Fr X 5 
1 1 2 2 5 
h= =—-rılnr (1-+Inr) (rt = x? + 2) 
87 27 


Die Hauptarbeit bei der Durchführung der Maple-Synge-Methode bildet die Berechnung 
sämtlicher Skalarprodukte (Tabelle 2). 
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Tabelle 2. Skalarprodukte 


ü, %, Q; 
JRR ER BEL DEE EEE BEER VERRRERSBEEE or mre  n 
u +41,13274.12287 — 32,00000.00000 —  21,33333.33333 
2 —32,00000.00000 +176,39644.73723 — 164,57142.85714 
Üg —21,33333.33333 —164,57142.85714 -+1245,05568.089 


JEDER RER EZ EB a een lien un 


“= 


2 


® 


>} 


FR 


vg 


Ö, +388,14316.536 —27,63306.3470 + 69,81385.7337 
OB — 27,63306.3470 -+47,19882.4140 + 11,60135.1563 
® + 69,81385.7337 +11,60135.1563 +176,65636.398 


Up 

Ug -+0,16978.85828 -+-0,09256.66367 

h -+-0,09256.66367 --0,05520.57865 

% +2,23746.29935 -+1,24889.93043 

Ülg —3,20000.00000 —1,56333.71080 

üs —1,90476.19048 —0,99452.53742 x 
DM —+2,71246.79278 -+2,00000.00000 

d, --0,29930.18549 0 

23 +.0,73223.34287 0 


Nach dem Schema von Peach [4] wurde die Orthonormierung vorgenommen. Tabelle 3 


enthält die Orthonormierungskoeffizienten. 


Tabelle 3. 


I 


U, 


Us = + 0,04542.76597 


+ 0,15592.15607 - 
%s = + 0,06320.53656 - %, 
», + 0,03749.12898 - ü, + 0,03135.20109 - %3 


Orthonormierungskoeffizienten 


An 
u, 


ü, + 0,08124.40609 - 


v, = + 0,05075.79447 -%, 
v, = + 0,01058.56220 - 
2, = — 0,01638.30138 - 


v 


©, + 0,14868.91545 - 


ö, — 0,02914.44885 


Ö, 
-d, + 0,07954.87252 - 9, 


Mit Hilfe dieser Orthonormierungskoeffizienten werden die Skalarprodukte mit u, undh 
umgerechnet (Tabelle 4). 


Tabelle 4 Transformierte Skalarprodukte 


ug h 
Ban +0,16978.85828 +-0,09256.66367 
h +0,09256.66367 +.0,05520.57865 
u | +-0,34886.87220 1.0,19473.03287 
| —0,11856.13284 —0,04807.47181 
u | —0,07804.75359 —0,03305.73224 
Yu +0,13767.92971 +.0,10151.58894 


-+0,07321.60999 


| +0,02117.12440 
Vz -+0,00508.68368 | 


—0,03276.60276 


Die in (22) auftretenden Werte sind hieraus nach (11) und (12) leicht bestimmbar. Es 
ergibt sich, wenn man von den u, und v, nur u, und v, benutzt 


lv, — 0,03860.795| < 0,01425.815 . 


Unter Verwendung aller Ansatzfunktionen erhält man 


| lv, — 0,03171.151| < 0,00271.548, | 


d.h. der Fehler ist kleiner als 00. 


Wenn auch dieser noch recht ungenaue Wert nur durch relativ großen Arbeitsaufwand 
gefunden werden kann, so sollte doch betont werden, daß Methoden zur Fehlerabschätzung bei 
der Plattengleichung bei weitem nicht in wünschenswerter Anzahl existieren. 


Abschließend danke ich Herrn Prof. Dr. L. Collatz für die Förderung dieser Arbeit und 
Herrn Dr. J. Albrecht für wertvolle Ratschläge. 
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U Advancement of a Compressible heat conducting fluid 
3 | over an infinite flat plate | | 2 
By A. S. Gupta (Indian Institute of Technology, Kharagpur) I 
Die Bewegungsgleichungen des instabilen Falles, wenn eine komprimierbare, zähe Flüssigkeit unter dem | k 


Winkel null eine unendliche, ebene Platte trifft und auf ihr mit verschwindendem Druckgradient fortschreitet, AN: 
werden integriert, wobei sich Semi-konvergierende Ausdrücke für den Fall ergeben, daß die physikalischen Kon- j 
stanten von der Temperatur abhängen. Einfache Ausdrücke werden für die Verteilung von Temperatur und Ge- 
schwindigkeit in der Grenzschicht, den Oberflächenwiderstandsbeiwert und ihre Abhängigkeit von den physi- 
kalischen Konstanten erhalten. 


The equations of motion in Ihe unsteady case when a compressible viscous fluid with vanishing pressure 
gradient flows past an infinite flat plate have been integrated in semi-convergent expressions when the physical 
constants depend on temperature. Simple expressions are obtained for temperature and velocity distribution 
in the boundary layer, the drag coefficients and their dependence on physical constants. 


Les equations de mouvement du cas instable sont integrees quand un ligquide compressible et visqueux (sous 
l’angle zero) rencontre une plaque infinie et plane <—) ei y avance avec un gradient de pression minime; pour le 
cas que les constantes physiques dependent de la temperature des expressions semi-convergeantes en resultent. On 

" obtient de simples expressions pour la distribution de la temperature et de la vitesse dans la couche limite ainsi que 
pour le coefficient de resistance superficielle et leur dependance des constantes physiques. 


YpaBHeHHA NBIHKeHNA UHTETPNPYIOTCAH MIA TOTO HeyCToiYHBOTO CIIyYaA, KoTNa BA3KAA 

_ GYEHMACMAaH FKHNKOCTb BCepeyaeT IIIIOCKyIO ÖeCKOHEeYHYIO IIIACTHHKY ION HYJIeBbIM YTJIOM M 

Tey&T. IO Heü IPM HyJIeBoM TpannueHTe naBıeHun. IllosyueHHble BbIPasReHHA CXONATCAH AacHM- 

ITOTHUYeCKH, KOTNaA PH3NYecKMe KOHCTAHTEBI 3ABHUCAHT OT TeMmeparTypbI. IlosyueHlo IIPocTbIe 

3 BbIP&#KeHUA MIA PacıpeneneHnn TeMIepaTypbI U CKOPOCTeÜ B IOTPAHHYHOM CJIOE, IA JTO00- 
E- BOTO CONPOTHBAICHUA U IIA UX 3ABUCHMOCTH OT PUSHYECKMX KOHCTAHT. 


Introduetion 


The flow of a compressible heat conducting fluid in the boundary layer has been investigated 
by Emmonsand Brainerd [1], Crocco [2]. Nigam [4], on the other hand, has considered 
z the advancement of an incompressible fluid over an infinite flat plate at zero incidence without 
4 taking heat transfer into account. In this paper the corresponding problem of a compressible 
fluid with heat transfer has been considered. The boundary layer equations of the unsteady flow 

with vanishing pressure gradient have been transformed with the aid of suitable similarity 
variables and then integrated in semi-convergent expressions following Meksyn [3]. The form 
of the differential equations obtained are similar in structure to those of Emmons and Brai- 
nerd [1] but the argument here has a different meaning. 
It has been found that when Prandtl number is taken unity, temperature becomes a 
quadratic function of the veloeity, a result similar in the steady case given by Crocco. The 
- essential conclusion of this paper is the growth of the temperature profile in a direction opposite 
to the advancement of the fluid over the plate. Expressions for the velocity and temperature 
distributions in the boundary layer are obtained for the case when the physical constants depend 
on temperature and the Prandtl number is not unity. Finally the drag coefficient and heat 
transfer have been calculated. 


Equations of Motion 


An infinite plate defined by the equation y = 0 is placed on the zx plane. The axis of x 
is taken along the length of the plate and the axis of y perpendicular to it. At the time ! = 0 the 
fluid occupies a vast expanse of space, for which x < 0 and the front surface of the fluid, the 


RENNER TUR TUNER EN TERM UF MIT SPTENENEE TO, 


2 plane x = 0, starts moving with a constant velocity u, in the x-direction. After an interval 
. of time { the stream has advanced a distance x = u, ft along the plate. The front surface at any 
E instant is defined by the equation x = u, i which is always kept moving with a constant velocity 
+ 


Sa 


} De: Z. angew. Math. Mech. 
350 Gupta, Advancement of a Compressible heat conducting fluid ete. a. 37 Nr. 910 Sept./0kt.1957 


ı.. The components of velocity are non-existent for x > u, { because the colum of fluid has not 


advanced so far. 
Approximating on the basis that ö, the maximum value of y in the boundary layer is <x 


and that v < ın, the equations of momentum in the boundary laver in the absence of external 
forces, are 


ou ou en & al nr N 
o(% + u or -- D oy = dy u ou a Yale 1 ( ” 
=0[0)<: . .... 
oy 3, = 


The equation of energy gives 
GE oT = E 2 Ban 0 ( ) 
- —+-—|) = x 2 Kae 
el tus +0 2, Er F\oy herr oy 8) 
where c, is assumed constant and x is the thermal conductivity. 
The equation of continuity is 


oo, 2 0 
tet en rear En 


The equation of state is 

pP Hol, ee ee ee 
Here, (2) gives the order of magnitude of Op/oy, that Op/ox can be taken as negligible is an assump- 
tion appropriate to the problem at hand. Using the fact that p is constant across the boundary 
layer, (1) and (5) give 


oT oT oT & = 
z a ed ee (6b) 
ol Nr | oy O2 0 (6) 
The boundary conditions are 
oT 
— _uzut I), — —0 D3 
ale N N (a) dy or (b) || m. 


at eye Nora 2 ne 11, 


In what follows the subscript 1 refers to conditions at infinity, O0 to conditions at the wall and 
dashes denote differentiation. 
Let us define the dimensionless quantities 


_ di SrU Bere 4 I ra 2 ER 
“ 0 a u en (nt Ay h 5 9 U, j “= Ta = 
We also take after Emmons and Brainerd [1] 
En N 
% Kı 
Making the substitutions (8), the equations (1), (3) and (6) become 
dr du dU 
20,0) +@ ib ar a ee; 
dU av. 2nder de dU\? 
ars DE er | 3% 
" dn ir < dn o „® n ne 2%) > au! 
dO dU dV 
Pe a a a 
( +nLl n) u (n in = n) (0) . (12) 


respectively, where o is the Prandtl number u C,j% Y = Cy[C„ My, is the free-stream Mach 
number u,/a, , a, being the acoustic velocity in the undistrurbed stream and b — ( —1)M*. 
Introdueing a variable &=n— n U—2V, (10) — (12) become 


d/.dU dU | 
20,0} es 
er Re O2 ee 16 
2 a . Er 2 dE 
Rn A) v a AN TE er 
er: er 10 Re 
de dO 
fe T EEE 
Ab= L ) E Pe cn 


a 
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The boundary conditions now are 
VeVz0 an) Geh); 
0-1 ası, 4,» 00%, 


I 


a3 nm, 


el 
9=.09, DER, U ale 


I 


Integration of equations 
We consider the general case of heat transfer when is ©, constant on the boundary. The me- 
thod of integration has been explained by Meksyn [3]. 
The function U is.to be found from (13) written as 


dad 4 = | = = 
.eE +599 Mn: —() EN A ee (17): 
The solutions of the equations (13) — (15) are found in series involving the two parameters 


& (= u,) and ©, by expansion of U, © and in Taylor series about n = 0 and using the boundary 
conditions at n = 0. The successive er Ures are obtained by differentiating the equations 


(13)—(15). 
The series solutions are 

eO,n? 28 Bar a—& m 
, bo 9,0? 22 

BE: 2 NEIT a  ; (19), 

0 

[OR 7? 

I Bi er TR EL or re a ze 

7 9, 9 (20) 


For evaluation of parameters we retain only the first term of the three series in (18)—(20). Thus 
we have 


ea 
3098 20H 


Hal): 
From (17) and (21) we obtain 


n 
2 
U EN D, | eXp - oD dn . 
0 


The exponential under the sign of integral is a rapidly decreasing function whereas ® is a slowly 
varying function; we can take accordingly ® = ®, under the integral sign, whence finally 
n 


1 \d 22) 
U=o| eXp 4095 an EN ee Set ) 
0 
and from the boundary condition U(&) = 1, we find 
1 
ra RENTE INTER 29): 
208 ee (23) 
From (13)—(15), we get 
202.0” do (G) do. 
- 1—o) — en DU ih. (28). 
Integrating (24) under the boundary conditions 
+=0, o= 09 0707 


a, NO] 


where, however, either ©, or ©, can be taken arbitrary, we find as Meskyn [3] 


[e0) N 


e=1— 0 |w dn +-bo (U’) dn GESTERN RABEN. 
8%, ; \ 
N N 0) 
This gives the temperature distribution when the velocity distribution is known. Substi- 
tuting the value of U from (22) in (25), we get 


oo N 


! ü 0 i a | u) 
0-1 | ex 5 a) uerEze | exp | non) a Gr 
: : 


n 7 


dn 
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Replacing the inner integral of the third term on the righthand side of the above equation by 
the error function erf x defined by 


% 
erf = [ern 
ym 
0 


we get after using (23), 
[6°) 29 3/2 
f. j n en BRZER j = en Be EC a. ;: 
9-1) ol 40 Eee exp 10 | 7: LE | 


n 
Since either ©, or ©% is arbitarry, the relation between them is obtained by putting n = 0 
in (26). Assuming that the plate is insulated (©, — 0), (26) gives 


ra 2 2 —— fr 1/2 
9=1+ Es Fer) et A 
en I 20: | 
Hence when n = 0, ©, is given by (when ©, = 0) 
1/2 “ u 1/2 
%=1+ m) | exp (— £) - erf f E = a N (28) 
0 


sirg (23). The integral in (28) has been evaluated numerically for different valus of Prandtl 
ımker. The following table gives an idea of the dependence of ©, on o for the free stream pro- 
pertyb=4 when &, =. 


Table 
o % 
‚25 1.659 
.133 2.247 
1 2.482 
1.2 2.629 
Again differentiating (27) with respect to n, we get 
do boa 70 n(2 — o)!" 
— 5 eXD _ )- ent = Eu ee 
dn (2 — 0) 4 9 9 Or 


which is negative for n > 0. Hence the temperature decreases from the plate for given values of 
oand ©. 


Also from (14), we have 
n n do do 
; DB — — B,|-— 
1 dU dU\: i dn Ab ) 
0 


[07 


which gives V when U and © are known. 
When o = 1, (25) admits of direct integration and it gives, using the boundary conditions 


(16), 


9, b 
et TEEN 3... ul 
Putting 7 = 0 we find 
oe ee] 

Een 5 Mi. a ee 

Eliminating ©,/x between (31) and (32), we get 

y—l1 —1 

=1+!5—M, 9=1 +7 Mi. se 


The skin friction r,, is defined by 


(na 


. (84). 


N 


z 

2 
” 

.. 
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The drag coefficient C, is given by 


v, | 12 


Now since u, £— x has the dimension of length, we can define a local Reynolds number Re 


et u, (ut—x) 
v ” 
Thus 
C)= 2D,e (Reale. Re) 


In conclusion I thank Dr. G. Bandyopadhyay for suggesting the problem and his kind 
help throughout this work. 
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Zur Abschätzung des Restgliedes der Mittelwertformeln 
zur genäherten Quadratur 
Von H. Schlechtweg in Kettwig 


Zur Abschätzung des durch die erste Ableitung der zu integrierenden Funktion dargestellten Restgliedes werden 

die Teilpunkte x,der Quadraturformel zusammen mit den Nullstellen der periodischen Bernoullischen Polynome 
B, (&, — x) und B, (z, — x + 1) herangezogen. 

The points of division, x,, in the formula for quadrature are used in conjunction with the zeros of the periodic 
Bernoullipolynomüals B, (x, — x) and B, (x, — x + 1) in ordert o obtain an estimate for the remainder given 
by the First derivative of the function to be integrated. 

Afin d’evaluer le terme residuel represente par la premiere derivee de la fonction a integrer les points partiels x, 
de la formule de la quadrature sont employes ensemble avec les points zero des polynomes periodiques de Ber- 
noulli B, (u — x) B ( v —c-+1]). 

Ina OMeHKU OCTATO4HOTO yUIeHA, CONepzKamıeroO IIPOM3BOAHYIO UHTETPNPpyeMofü DYHRIHN, 
HUCIHOJAb3YIOTCH BMECTE C TOYKAMH AECJIEHHA B KBANPATypHoA BopMyale x, HyJIH HEePHOANYecKuUX 
DONAHOMOB BepHyasu B, (, —x2)u B (u —cH1]). 


Stellt man den Wert eines bestimmten Integrals näherungsweise dar durch 
al N 
Ri I) de z 2 AsTE) ea ae e e- er e i 


wobei die x, gewisse im Intervall O<x<<1 liegende Abszisse sind, so läßt sich bei im Inte- 
grationsintervall stetig differenzierbarem f(x) das den Unterschied zwischen Näherungswert (1) 
und exaktem Wert angebende Restglied in der Form darstellen?) 


N 2 Se 
Be SA STOSS) 
v—1 0 


wobei B,(x) das erste periodische Bernoullische Polynom ist. Die auftretenden Integrale 
lassen sich zwar mit dem ersten oder zweiten Mittelwertsatz leicht abschätzen; jedoch wird 
hierdurch in praktischen Fällen zuweilen eine geringere Genauigkeit der Näherung (1) vorge- 
täuscht, als sie tatsächlich vorliegt. Im folgenden wird daher das Restglied (2) in etwas schärferer 
Weise abgeschätzt. N 

Da die periodischen Bernoullischen Polynome B;,(x) aus den Bernoullischen Poly- 
nomen B,(x) durch die Festsetzung 


B,(x) = B,(&) für 0<ı<1 
By(& +1) = B,(«) für alle & 
hervorgehen und 
Bar > 


u 


)) vol. z. B. B. W. Tollmien, Z. angew. Math. Mech. 29 (1949), S. 193; daselbst weitere Literaturangaben. 
23 
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ist, so läßt sich das Restglied (2) in der folgenden Form schreiben, in der das Argument von 
B, stets zwischen 0 und 1 liegt: 


-5 ZBNE: Ma) Ba 2) + [def Br + Y. 


was gleichbedeutend ist mit 
f | or 1 
= 3 8 | dx f(x) (.—2— =, + ja: 6) ( —ı +5) : 
N () %, 


ist 2 = für» <n’, so sind in 


r, | 2,+5 i 
RA [dee 4,)@W+ 1 (v4 1a) r@- [ dee, 5)r@ 
BE art: 
+34, Mn ro j del@— 2,4 3 ar jeß+s -2)F@r @. 


die im Integranden mit /’(x) multiplizierten Faktoren stets positiv. Jedes der hier auftretenden 
Integrale läßt sich zwar sofort mit dem ersten Mittelwertsatz abschätzen. Jedoch soll hier die 
Abschätzung etwas schärfer durchgeführt werden. 


1. Große Anzahl Teilpunkte bei beliebiger Intervallteilung 


Die Restgliedformel (3) läßt sich bequem noch genauer abschätzen, wenn man die dort 
vorkommenden Integrale zerlegt in Teilintegrale zwischen den einzelnen x, und in diesen zur 
Abschätzung die Differenzen der Funktion /(x) benutzt. Es sei x,, der letzte Punkt der zur 


Näherungsberechnung des Integrals benutzten Intervallteilung vor x, + = ausschließlich 
dieses Punktes selbst; analog sei x,_ der letzte Punkt der Intervallteilung vor x, — 5° aus- 


schließlich dieses Punktes selbst. Bezeichnet man die sechs in der Restgliedformel (3) vorkom- 
menden Summen, abgesehen von den dort stehenden Vorzeichen mit Ry Ru Ri; R 


20> Ry 
R;,, so hat man 


n'’—l v—1 
ee ‚2 de (8 — 3, En re, 
= a 
2+5 
i n’—l b, 1 
Ru=Rr+24 | del, +3 —2)f@), 
v=1 N d 
41 
—1l »„t—1 f 
Men AaN | au (r + — :)7@), 
v—i A=V 2 
R el ] 
En | du (x u 5) re ZA, | de (0a BR): 
= : ö ag 9 I) 9 
2. t5 ®n 


7 
A 
u” 
di 
ar 
ne 
F 
e 
3 
er 
w. 
I 

2 


vw 


DEREN 
! 


in 
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1 
4,5 
2. 1 
Rg= Run + 24, | (5 — 2)/6), 
v=n! 2 
%,— 
ZAHL 
n v—1 
Ru= 3A | 5-2) r@), 
vn" i=0 
2 
%y-+1 
u ji 
Rua=RyvH+t DA, | dx az, + = Io) 
v=n 
AR 


+1 
v—l 


Rr= 343 |de-n+s) ro, 


v=n I=r1,+1 
7 


1 


Ra=Rr+ 234, [del +3 —2)f@), 
“+1 
n—1 1 
Rr= ZA, 2 da (x, + a F.@)* 


2 


Die hier vorkommenden Integrale lassen sich mit dem ersten Mittelwertsatz bequem abschätzen 
wobei dann f. (x) auftritt an gewissen Zwischenstellen, die von den jeweiligen Integrationsgrenzen 
abhängen sowie außerdem noch von den im Integranden als Parameter vorkommenden x,. 
Bei der Abschätzung kann man bei numerisch nicht genauer bekanntem /’(x) in jedem einzelnen 
R;ı die im gleichen Integrationsintervall vorstehender Integrale genommenen Zwischenwerte 
von f’(x) für alle Parameter x, des Integranden durch die gleiche Zahl abschätzen. Dieser Fall 
möge hier interessieren; es darf dann also im folgenden die Abhängigkeit von den genannten 
Parametern x, des Integranden außer acht gelassen werden. Allerdings dürfen auch dann noch 
gleiche Integrale, die in verschiedenen R,, vorkommen, im allgemeinen nicht durch dieselben, 
aus Zwischenwerten von f’(x) entstehenden Zahlen abgeschätzt werden, da die R,, mit teilweise 
verschiedenen Vorzeichen in den Rest R, eingehen. Aus diesem Grund werden die zur Ab- 
schätzung von R, zu verwendenden Schranken von f’(x) in verschiedener Weise bezeichnet wie 
mit Di 341, Di,a+ı usw., wobei dies Schränken von f’(z) sind jeweils im Intervall 2, <x< x.ı. 
Die zu entstehenden Schranken von R,, und R, werden bezeichnet Rf, und R}. 


In der obigen Darstellung von R;,; ist, falls n’” als Teilpunkt der zur numerischen Inte- 
gration benutzten Intervallteilung nicht existiert, das Glied » = n’ wegzulassen, während in 
dem zweiten Summanden von R,, und den Summanden von R,, dann n’” = 0 zu setzen ist. 
R,, Ri bis Ry} lassen sich nun in folgender Weise umordnen nach Anwendung des ersten Mittel- 
wertsatzes: 


n’—2 | nn .n x 1 
Ri = 2 Dan au) 2 2 + an 
al 2 
Die Schranke 
= I RE» x 
=: nal) An 


läßt sich für hinreichend große n, die die Bedingung 1 > n’— 1 erfüllen, in der Form schreiben 
w—lj—1 


1 DE Bi ar (+1 — 3) Yıa 
4=1 


mit 
A 
> Als + : m für A.=1,2,...,n" —2 
vl 
RA ME: 5 
Ss A,[® \ u füri=xt, »r+1l,....,@&+Dr—1, 
ah 2 = wenn #0, 1, 2, Henne 


23* 


1 
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wobei 0+ = n’— 1 gesetzt und «= 0 wegfällt, wenn 1* = n’—1 ist. Entsprechend erhält 
man 


n—1 


Rı= 3 Dr (Urı — U) Yına 
A=1++1 


Se Ts =) für a = ar ae Ne aa 
=1 


2 wenn <=1,2,...,n" —2 


a ze 
= on ee 3 für 1 (mn 11 Sn 
11 
und 
n—1l 
An > Di, 141 (Ir — %) YIIIA 
i=0 
San - rt = für =0,1,2,...,n7—1 
ven = 
u ilesı Au 2) für 1 +5 red 
Re z = wenn <=0,1,2,...,n—n’ —1. 
Ferner wird 
n—1 
Ry= 3 Dir &rı — %)yıwa: 
A=n +1 


Wegen(n a) = (in 7x: 1) undn <n,d.hnen lot ie a <n. 
Bezeichnet dann Z,+..- die Vielfachheit, mit der ein festes (n’ + x) in der Reihe aller (n’ + x) 
vorkommt, so wird 


f n'+% R Pr 
Zus 2 AR, 4 L) für aa HH la 42h, 


wenn <=0,1,2,...,n—n’—1 
ME Se ur ) rien art 
3 a4) Tür an nt er ann 
ver 2 2 


Dabei sind Summen, deren höchster Summationsindex kleiner ist als der unterste, wegzulassen. 
Ferner wird 


n—1 N N £ N r 
n; 2 Du (ar — U) yvas wı= 2A, (# us nn eh 


2 1 
Da &+.1 entweder gleich x, w oder der nächste Teilpunkt darüber ist und ui < 3° 
wird —1<1* +1 und somit n"<1+ +1. Ferner ist zw_1y)_ı der zweite Teilpunkt vor 
Inı + a Ist n genügend groß, so kann der in folgenden zugrunde gelegte Fall betrachtet 


werden, daß &w_1)+—ı > Xı++11st. Dann beschränkt man die Betrachtung auf Fälle, in denen 
folgende Ungleichung gilt: 


"<it+1l<a'—Ddr—1. 


1 I! 
Aus. <a, a Ss folgt n”"<n’ und demnach"n An aa %,- der nächste 
Bi 1 
Teilpunkt unter mn und x%„- der nächste Teilpunkt unter ee ist, wobei n’ <n, 
so ist n"<n” und für geeignet großes n auch n’- +1<n-—1. Damit erhält man 


in dem hiermit festgelegten Bereich von n für die Schranke R* des Restes (3) die folgende 


Mittelwertformel, wenn 9(a, b) eine Funktion bedeutet, die für a<x <b den Wert1 hat 
und sonst verschwindet: 


_ 
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: Rı= Doı "u (Yo —%0) + 23 — Dy,,41 (41 — 2%) 90, g(L,.n’ — 2) 


Z: en — 3) yvg(l, ("’— 1) —1)+ au + (e + E g(l,n’—1) 


a 
2 


- 1\2 A, ; 
a ®,, v+1 (G41—%,) Yır,g((n 1% N “a a) 594, n —n 


+ ds or —u)ym,g(l,n —1) + en, et | — = — I. g(n’, n) (1.1). 


2 
— ei (&,+1 — %&) yıvvg(n' "+1, n—1)—® | + En Ar m, n) 
x V— zer 2 2 


+ D, 41 (1 Ey %,) ver) + 


N \ n—l = 
+ 9,11 —2,)3 A, (r — =) — dı1l —x)3 Au: —) 5 
v=n' u 4 


ir 
Dabei bedeutet ®,,ı bzw. ©} | eine Schranke von /’(z) im Intervall 2, <ı<1. 


2. Das Newton-Cotes-Verfahren 
Hier werde zuerst der Fall ungeraden n, dann der Fall geraden n betrachtet. 
A.) Im Falln=2 m — 1, also 2, = ee ist n = m. Aus (3) erhält man, wenn 
man die Integrale zerlegt in die zwischen je zwei aufeinanderfolgenden x, erstreckten Teilintegrale: 
A+l x A+1 


2 m—2 2 m—2 


m—l Be mv v+m—3 
Rma= Ai | + a)f@a+ 3 
= , =v—1 


v»+m—2 y 
2m—2 | man oz 


SS 
3 


a 2 „+ In ; 
— IB: eo f(&) dx 


A=v+m—2 = 
2 m— 
A+l A+1 
2m—2 2m—2 
2m—l v—m—l zT w—m v—2 z m 
Be A | | =, '()de— 3 (x) dx 
en - 5 “ 2m--2 "@ A=ym a 
Im 2m—2 
A+1 
2m—2 
2m3 [( +m—2 } 
- ? —— —z)f(x) de 
2 m—2 


wobei der Akzent an dem Summenzeichen bedeutet, daß im Fall» =2m—-1 die über A 
erstreckte zugehörige Summe durch Null zu ersetzen ist. Schätzt man im Intervall 


meer. a die Funktion f’(x) ab durch B,.,ı und — f(x) durch — 94,441, SO 
2m— 2 2 m—2 

liefert die Anwendung des ersten Mittelwertsatzes, wenn man die vorstehend vorkommenden, 
positiv genommenen Summen über » der Reihe nach mit S, bis S, bezeichnet, die folgende 


Schranke für Rem: 


TE ESSEN SET STE Se. a. ee te (el) 
mit 
= Bea, 2 me An Pi,AH1 > 53 =; An : un Te = Bi 2 Dur 
Se > IA ar ar Paa41, Sa =) A, aA Eee I Di, 


Z. angew. Math. Mech. 
958 Schlechtweg, Zur Abschätzung des Restgliedes der Mittelwertformeln „4.37 Nr. 9710 Sept. [0kt.1957 


Die Umordnung der letzteren Summen nach den ©, ,;ı und 9,7+1 ergibt 


m—3 2(m+Hi—N)+1 (0.2) 
En A, 2 ee .2). 
= 2pamı 2, 8 (m — 1): 


Bedeutet Max (a, b) die größte und Min (a, b) die kleinere der beiden Zahlen a und b, so erhält 
man 


A Rn 2W+m—i—-)—1 
= Pr en »=Max (1, A+3—m) ; 8 (im— 1% 
Unterteilt man die letztere Summe in 
m—2 a+1 2 +m—1—2)—1 
v Arm v=Max NET = ini 
ü 2m m me 
A, > 
T a Paar »=Max PER 8 (m Ay 
so wird ss 
M— —1—2)—1 
HH 2 tm —A—-3—1 mi - 2, +m—4 
= 79m 24 sm I + 3 Pan 23 nn gm TE CB). 
m wird h DB Min (m—1, 242m) 2(A +2 —v— m)+1 
re rt 2, ei 8 Gm? | 
oder r er ne Er | 
2 m— —m Sm 
= — A, ie (2.2c) 
ee 8 (m — 1): 
sowie re %K »—_1 
m—2 M— Ye sme zu 
ES SA en 9 2.20 
Si = IP 2, a 8 (m — 1)? ( ) 
Die ana Wan der in folgender Form darstellbaren Summe S* 
2m—3 Min (2 m—1, A-+m) 9 (A m — 9») — f 
Se — A, 
: Pa ar: v»=Max (m, A+2) Slam 1). 
in die Einzelsummen 
m— 2 Min (2 m—1, A+m) 2(1 — m—v) IL 1 2 m—83 Min (2 m—1, A+m) > (A + m—v) _Q f 
—— N A, 
u „2 Soear en s(m— 1} 
liefert 
m—2 m+i 2A+m—)+1 2 m—3 2m— 2A +m—)-+1 
5; = 2, marı 2 As 8(m = 1)% u ea pr, Ko. , As 8 (m = 1) . (2.2e). 
Entsprechend wird 
2m—3 + 2e+m—i—)—1 
= >’ A, u. u nd 2.2.8. 22H: 
Sa R er Pr 2 2A, 8 (m — 1)? ( ) 


Einsetzen der Ergebnisse (2.2) in (2.1) führt auf die folgende Mittelwertformel für eine 
Schranke des Restgliedes: 


2 m—3 
R: nm = 2 X Dar = Pr DER EI RE (2.3) 
mit ai 
et 2m 2 —m—A)—1 
Hal De ne ya Page Be 
mw ba 8 (m — 1) am (m m a Er 
= 
ae 2 Hm del 
I A 4 r PERLE Bi 
BE az ) frrysem—ln ra 2m (2.3a). 
A-+-m © er 
Sm (Be = Ze) +1 fir an rare 
RE a 8 (m — 1)? 
A—(m—2) 9) (A- ae m—v)+ 1 a 2) (A + m—») au 1 
>: A, = BE X + Ü — = 
= 8(m— 1)? re A 8(m—1)? für m % m, m am 3 


Für die Gewichtea,, ß, der Fehlerformel (2.3) gilt die Symmetriebeziehung 
% =, für A u = n — 20 (2.4). 


san ı, 


a a ER Bal 1 önn 1 3 3a m nenn © u nn 
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Ist in dem hier vorliegenden Falln =2m— 1 nämlichA<m— 2, so betrachte man Pau 3 
welch letzteres man aus ß, für > m—.1 erhält, indem man } durch 2m—3— 4 ersetzt; dies 
ergibt 
nF 2 (mn —))—1 2m1 2(3m—A—»)—5 
Be Ar A, 
reg 2 AT 
Nun gilt für die Gewichte A, der Integralformel (1) in dem Fall des Newton-Cotes-Verfahrens 
die Symmetriebeziehung 


A Aa RE A EE (2.5). 


Setzt man dies in fgm-3_, ein und ersetzt dann den Summationsindex » durch 2m —», so be- 
stätigt man in der Tat (2.4) fürA<m—2. Für «, mit A> m-—-1 betrachte man f, mit 
A<m-—.2; dann wird mit (2.5) 
er 2—m+A)+1 
Sn A 2 
Ersetzt man wieder » durch 2m —v, so erweist sich hiermit (2.4) auch in diesem Fall als be- 
wiesen. 


Di Palln= Zmistr, = und n’= m +1. Die Ausgangsformel (3) für den 


2m—1 
Rest liefert dann, wenn man die Integrale nach allen Teilpunkten x, und allen Punkten x, + D 
unterteilt: 
Ram —— 
21 et 
5 A Tl L Dim ap i| Pa)dx Te a2 m=—3 e) ach 
ver n i=0 2 (2 mil Nee > (2 Mi 1) 
2m—l 2m—l 
2v+2 m—3 vtHm—1 
2(2m—1) 2m—1 
| er am ) | ( ep. 
, { 3 (@m— 1) / ns Om) 
2m— 2 (2 m—l) 
A+1 
2m—l 
2m—2 a 
= TI.dr 
ri | | 2(2m—|]) 7 
2m—l ) (2,6), 
Ari 2r—2 m—1 h 
2m—l 9 1 2 (2 m—1) 4 R ' - 
2m v—m—2 oe dm — A ln A N 
+ A, R= | De, («)dc+ | | e) f(&) dx 
; en en 2 2 (2m FE 1) v—m—1 = z Su 2 
mi 2m 
vm A+1 
2m—l 2m—l 9 | 
i re en | | 29» —2m— } 
= (x) de — X = x) dx 
! ( Te Ne 2 Sana 
2v—2m— N 
2 (2 m—1) 2m—l 
A+1 
zm—l » 3 
ee. 2» +2 m — j 
2et PHlie)dE,, 
2m—1 J 


wobei der Akzent bei der ersten Summe zum Ausdruck bringt, daß im Fall » = 1 nur das Glied 

— () zu benutzen ist, der Doppelakzent, daß im Fall» = m + 1 ebenfalls nur das Glied A = 0 
zu benutzen ist, und der dreifache Akzent, daß im Fall v = 2 m die zugehörige Summe über A 
durch Null zu ersetzen ist. Die hier auftretenden Integrale haben zum Teil wie im Fall ungeraden 
n als Integrationsgrenzen die in der Mittelwertformel (1) für das Integral benutzten Teilpunkte; 
zum anderen Teil treten als Integrationsgrenzen aber noch die Halbierungspunkte der 
Teilintervalle der Interpolation auf. Es mögen zunächst die ersteren Integrale betrachtet 
werden. Bezeichnet man die in (2.6) vorkommenden Summen der Reihe nach mit 5, bis 51, 
und wendet den ersten Mittelwertsatz an, so lassen sich die genannten Summen in folgende 


Ex 
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Weise umordnen, wenn man gleich die sich ergebenden Schranken hinschreibt: 


& m 2m +Hi—v+1 I I A ma Nee 
S] = 24,2 "Omi Pr,A+1 en Pirat 2, em @m— 1) | ! 2.72). 
5 m »+m3 m —A+v— 2 Du ” Min(m,A+1) m-—= + ye 
5 = 2 Ay. " Om— ii, D,41 = 2 ER re | 
Da für < m — 2 stets Max (1, + 3— m) = 1 ist, so läßt sich S$ in folgender Weise schreiben: 


m—2 A+1 m—A+v—2 "os = A mA Tr (2.7b). 
an Aufn 2,4 Female ne, ua Kal Mi RER Des % 


Ferner wird Min (mat m) +2 
m 2 m—2 mo) ra TE Pal ee 
Ze ea a 


—ı A=v+m—1 


oder 
e Zm—2 Sn ) my + 7 
Na en en ET 
s Br vn vl "(2m — 1)? 


2 TEE | 


In der gleichen Weise erhält man 
my m-iA—]| se = v—m—)—]1 


= 2 AS Tom Ban ID 


v»=m-2 i=0 v=/+m+2 


2m v—2 m+A—v-+]1 2 m—2 Min (2 m, A+m) m ti—r 1 


Ken ” — Te 
= 2 A 2 amt, alla SM ae a 
oder 
m—l mA m = A en, - 1 2 m—2 2m m + y} een = 1 
k z — > = 3 A, > . . 2.7e 
55 va Pr, A+1 a A, (2 m — 1)? - Par 2, (2 re 1)? ( ) 
und 
sr SA a era, = % m—/A+r—2 (2.70) 
Te ee ee na Bi 
Der verbleibende Teil von (2.6) ergibt zusammen: 
j 1 m 
ö ar 55 Si 5 Ss = 8 (2 m Er 1% = A, (Pins, »+m— == P,+m- >, 


2m 


a 2 A, (2, m—l, v—m oc m en ar 


v=m+1 
Dies läßt sich, wenn man in der ersten Summe » durch v— m + 1 ersetzt und in der zweiten 
Summe » durch v» + m, in die Form bringen 


1 2m—l 
er re >" y,[® — a e 
3 8 (2 m — 1) ;Zı ’ | 1, 9.) 23) 

mit 
Ana türgape tl arm 
Am+ R A 
A, lün wem mel, 2m 
Dabei bedeutet der bei der Summe in (2.8) stehende Akzent, daß das Glied v» = m zweimal zu 
nehmen ist gemäß (2.8a). Da wegen (2.5) die Beziehung A„,, = An-ıı gilt, so findet man 
füry =2m-) mi —1,2,...,moe Symmetriebeziehung 


Vom ii: > 5 Tal, u a ee (2.9). 

Aus (2.6) gewinnt man dann mit (2.7) und (2.8) unter Berücksichtigung von (2.3) für jedes n 

die allgemeine Mittelwertformel für eine Schranke des Restgliedes Te 
n—1 


1 = 
Er er ee) 


Y» == 


n—2 
N 

)=0 
wobei für ungerade n die a, ß, durch (2.3a, b) gegeben werden und y; = 0 ist, während für 
gerade n = 2m die y, durch (2.8) gegeben sind und 


I m—i+Vv—2 An v—m—/A—] a h 
ZT LE | 2 Ay ee ÜL — De) 
; v1 (2 m — 1)? a en ' 41,2, » » 
tn == 
; Ve u (2.11a), 
an -—— für = m mm See 


BEE ee 
w=i4+3—m (2m — 133 
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en a | 

Pa = Sara A—m—r+2 2m mHtA—v-+] 

ar En Terme 

Die Symmetriebeziehung (2.4) zwischen den a, und f, weist man wieder in der Weise nach, 

daß man zunächst a, für A<m—-2 und ß, für A > m betrachtet; dann wird für A<m-—2 
mIA m—v—/A 


Zn s3sm—A—v—1 
m > A, ee er A, 9 
Pam _ 2 (2 ll 22 rg (2 N — D- 
woraus man sofort unter Benutzung von (2.5) dann die Gültigkeit von (2.4) für A<m—2 
nachweist, wenn man in der zuletzt erhaltenen Darstellung von ßen _2_, den Summations- 
index » ersetzt durch 2m +1 —»v. Fürm— 1</A<2m-—2 erhält man dann aus der ersten 
Formel (2.11b): 


ame 0 are mass T, 
(2.11b). 
ir Am, mel oa 


3m—2—A LT EEE ER A+1 m-—A-+v—2 ' 
a — As +1 — 4 FPeren: = &. 
P ; N a (2 ZZ u): v=A-3—m : (2 U I 
Für kleine n erhält man aus (2.3a, b) und (2.11a, b) die folgenden Werte: 
= 6) y=4 
£ Be de 8A 72a 728 576y 
x x 
0 1 ) 1 
1 D) ) 2 3 3 
> 4 Bla 
2 2 1 N) 3 3 3 2 
ı 1 = 3 
3 3 2 3 
1 l 
7 6b 
90.4 1440 & 1440 ß y 288A 7200 72008 57600y 
x x 
i 2 37 54 ) : 19 151 200 50 
% 32 = 75 
55 34 0 5 188 150 75 
= 12 5 50 
ni 34 55 N) i 175 175 38 
B 32 — 50 
4 54 37 0 5 150 188 75 
1 7 4 2 
er 75 
5 200 151 50 
1 19 
| 
——| 3404 60480 60480ß y 
x 
- _ 1085 1433 0 
= 216 
i 1542 978 0 
4 27 
. 865 1657 0 
— 272 
z 1657 8365 0 
= 27 
5 978 1542 ) 
Be 216 
6 1433 1085 0 
1 41 


Die Anwendung der Methode dieses Paragraphen auf andere Integrationsverfahren liegt 
auf der Hand. 


Eingegangen am 22. Oktober 1956. 
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Seen ee 


The Forces on a Thin Aerofoil in Slightly Par: 

E By E.E. Jones in Nottingham, England BY 
werd j dıy ischen Bigenschaften eines dimmen gebogenen Tragflügels bes: 
ne a ns De in einem inkompressibelen nicht-viskosen 


5 | ‚ strömt wird. Unter der Annahme, daß die Abweichungen von der Grundströmung klein sind, wer 
Br! Effekte erster Ordnung abgeleitet, wobei das Geschwindigkeitsprofil als parabolisch angenommen wird. _ ; 


u The aerodymamic characteristics of a thin cambered aerofoil are determined, when the aerofoil is set at 


n cidence in a shear flow of incompressible inviscid fluid. Fürst order effects are deduced assuming that the de- 


viation from uniform flow is small, the velocity profile being parabolic in form. 


On determine les qualites adrodynamiques d’une voilure mince et courbee pour le cas que la voilure est l’objet 
d’un soufflage par un courant de cisaillement dans un medium incompressible et non-visqueux. En admettant 


que les variations du courant egal, soient mimimes, les effets du premier degre& sont derives, tout en considerant 
le profil de vitesse comme parabolique. 

VIccnenylIora a9poNMHaMHYecKuUe CBOÄCTBA TOHKONPOPHIIBHOTO Kpklla MIA TOrO CAyYan, 
Kora KPbLIO B He3ikMUMaeMON H HeBABKOH CPete HoNBepraercH KocoMy TeyeHnmo. IIpuunman, 
yYTO OTKJIOHeHUA OT PABHOMePHOFTO TeYeHHA MallbI, BBIBONATCH ABJIEHUA BTOPOTO IIOPANKA, 
IIpryeM IPHHHMaeTchH, YTO IPOPHNIB CKOPocTeü ABIACTCH TapaoolImyeckuM. 


1. Introduetion 


There are many examples in practice where a parallel flow of fluid is non-uniform to the 


extent that the velocity of the stream parallel to the flow varies in magnitude according to 

"position measured transversely to the stream. Such a flow exhibits shear characteristics, a 

typical example being the flow in the presence of a plane boundary, the velocity increasing with 

i distance from the boundary. In order to determine the flow it is necessary to prescibe the velo- 
city variation of the fluid, or the shape of the velocity profile. 

The case of linear variation in velocity, with a straight line velocity profile, has been investi- 
gated by several authors, in particular by James, [1]. This is an example of constant vorticity 
distribution throughout the fluid. An exact analysis is possible based on the Laplace differential 
equation, and a stream function exists. 

Another form of shear flow, more realistic in some cases, occurs when the velocity variation 
transverse to the stream follows a quadratic law, the velocity profile being parabolic in form, 

.e.g. such a. flow exists in flows through a tunnel or duct. In this case first order effects due to 
the variation from uniform flow can be determined provided it is assumed that the variation is 
small, so that the profile form is a parabola of large latus rectum. This problem has been in- 
vestigated by Nagamatsu, [2], who attempted to determine an expression for the stream 
function of a shear flow past a circular cylinder and a flat plate. It was assumed that the parabolic 
velocity variation of the stream, and the disturbance due to the body set in the flow were small. 
By using these approximations the analysis was again based on the Laplace differential equation. 
The solution derived does not however satisfy the hypothesis that the perturbation due to the 
presence of the body is small, in fact this assumption can hardly be valid in the case of a circular 
cylinder, and could possibly apply to the flat plate only for small values of the angle of incidence 
to the main stream. The conclusion drawn is that the problem has been over-approximated. 

In the present paper the sole restriction is that the velocity deviation from the uniform is 
small, and in this case it is found that the differential equation controlling the flow can be solved 
in terms of Mathieu functions. Similar types of differential equations occur in the theory of 
bending of thick plates, (3), and in the determination ofthe drag on a cylinder in a viscous fluid, (4). 
The results of (2) can not only be corrected, but extended to include cylinders having section 
profiles in the form of aerofoils of small thickness and camber, and the more important charac- 
teristics of the aerofoil, such as lift, drag, couple and position of the centre of pressure, can be 
determined. 

The basic theory as developed in (2) will be briefly repeated here, with suitable modi- 
fications. Orthogonal stream axes O'X, O'Y define the two-dimensional flow plane, the un- 
disturbed flow being parallel and symmetrical with respect to the X-axis. If the stream velocity 


profile is slightly parabolie, then the velocity components relative to the respective stream axes 
at the point (X, Y) are 


Yy2 
vr VD=0 we 


Here V is the velocity along the axis of symmetry of the flow, qis a non-dimensional constant of 
the first order of small quantities, and his a standard of length reference. 

‚If (u, v) are the velocity components of an incompressible inviseid fluid in steady motion 
relative to any two mutually orthogonal rectilinear x-, y-axes, then there exists a stream function y 


ah 
wi 


j 


Lau I wenn u nt 4 hr 
r j EA RTE y 
> a u 1 Zaun) 
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such that 
Y oYy 
uU=-, te Sa al a a Nee a 
oy 0x ® 


These satisfy the continuity equation, and hold in both the undisturbed stream and the stream 
disturbed by the presence of a cylindrical body. If y, is the stream function of the undisturbed 
stream, it follows from eqs. (1) and (2) that 


Y3 
w-—Vlr+l,) Pe u (3). 


As a direct consequence of Helmoltz’s law for the rate of change of vorticity in a fluid, it is 
known, (5; pp. 103, 109), that the vortieity is constant along a streamline of flow, and that 


a I a ee 


- where y is the stream function. It is possible here to introduce a perturbation stream function y,, 


where y = y, + y,, and at large distances, since y, and A%y, tend to zero, then 
Ay=f (Y) . 


qVY 
IE ee Sa ee (>); 


From eg. (3) it follows that 


hence from eq. (4), 
Ay = fly + yı) — po) > 
and the right hand side can be developed in a Taylor’s series, leading to 
2 
++ Be ii): 

In particular, using eq. (5), | 

OR q 02f Kr 

% (M+gYP)’ 02 V(R®+ qy2Pp)®' 


ee Yaz i 3 5 E 
If it is assumed thats is small compared with unity, and powers of q higher than the first 
DETSLTE IE 
ar lected, then — & —-, —, = ı— 
e neglecte en di? ya pay 
it will be assumed that y, is O (1), hence to the first order in q under these circumstances the dif- 
ferential equation of motion takes the form 


etc. By comparison with irrotational flow, when qg =, 


Ay, = U ua Le eignen at. Der, heise alte, Ma ee (7). 


If this differential equation is assumed to be the controlling equation of the flow everywhere 
in the fluid, then it is noticed that A%y,, a measure of the vorticity, and y,, tend to zero together 
at large distances as required, since there the effect of the presence of the cylinder is negligible, 
and the flow has the same characteristics as that of the undisturbed stream. A solution of eq. (7) 
will give an accurate representation of the flow pattern in the vieinity of the cylinder if situated 
near the X-axis, but although the boundary conditions at infinity are satisfied, in the intermediate 
zone the solution gives only an approximate picture of the flow. It is noticed here that y, need. 
not be small in order to produce a solution of eq. (7), as was assumed in (2). 


2. The Flow Equation 


An origin O is located inside the right section of a cylinder, with a profile in the form of an 
aerofoil, and orthogonal rectilinear cylinder axes Ox, O y are located in the plane of this section. 
The x-axis is inclined at an angle « to the direction of the X-axis, and if O is a distance d from the 
X-axis, then 

ee a Le Be or er (:) 


where (x, y) and (X, Y) are the components of the same point relative to the cylinder and stream 
axes respectively. This equation enables y,, from eq. (3), to be expressed as a function of x, y. 

The aerofoil contour is of the Joukowski type, defined by the curve corresponding to 
&E=0(in the net 


b? 
z=2+he- 


| (st+he)" 
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= [ = in, with (este 2m. The origin of coordinates O in the 
nn N Be of ne circle of ed of radius b, and h is a Be 
of the transformed circle, whose centre is at the point (X% Yo)» where 2 = % + Io. 4 ir : 
taken as the standard of length mentioned previously. The thin cambered ST ER is forme 
when x,, y, are small quantities, such that their squares and product are negligib B ara 
The aerfoil has a sharp trailing edge at the point corresponding to E= 0, n= m. e. 


{ dz A 
&=&,=in,, where (; is a zero of DE 0. Hence this leads to 
»„the=—b, | 
and by considering its conjugate, and eliminating £,, to the first order in x,, Y,, it can be shown 
that b=h—x,. The second zero of -— 0, givenby „the =b is not required and 


corresponds to a point inside the aerofoil contour. The trailing edge thus corresponds to al, 
defined by 


ee) ee ee 
The transformation of eq. (9), under the assumption of small z,, takes the form 
zen Ehe (k-20)er? men a Bee (41, 
ide] 


and corresponding to eq. (11), a function H = can be determined, where to the first order 


Inwr,,ts; IdZ 
H? = 2 h? (cosh 2E— cos2n) —Ahzye”: +2he: (z, ei" + 2,e®'”) 
— aeg er tze nz ra 
Under the transformation of eq. (11), the differential equation of eq. (7) takes the form 


Ben | 
Ge on men. u age 


A solution by separation of the variables is sought, and this is achieved if it is assumed that g, 
X Yo, are of the same order of smail quantities, so that products and powers higher than the first 
are neglected for first order effects. In this case, using eq. (12), eq. (13) becomes 


yı , Pyı 
oe ee KT A See 
and a solution will be formed of the type y, = g(n) G(£). The equations when separated lead to 
d2g 
an + (ar 20 00s 29 =0 2 ei ee le 
d?G 5 pe Ä 
ge (a tr2gcosh22)G=0 „u... 2 205 


and these are respectively the ordinary and modified Mathieu differential equations. Solutions 
are required for g(n) periodic in n of period 27 for a discrete set of eigenvalues of a, (6, p. 174). 
Such functions are the ordinary Mathieu functions of integral order denoted by ce,„(n), 
ce n+1(N), Seen+ı(M), Sean+2(n), defined in (6, pp. 21,22), and referred to there as ce,, (m —N): 
etc. These functions involve coefficients A7, B", which are known functions of q. Solutions are 
required for G(£) which decrease monotonically to zero as &—> -+ oo, and these are the mo- 
dified Mathieu functions Fek,(&), Gek,(£), defined in (6, p. 248), and referred to there as 
F ekn(8, —q), etc. These functions can be expressed as Bessel function product series converging 
uniformly in any finite region of the E-axis, including the origine =. 
The required solution of eq. (14) is thus 


Y = Z1C„Fek,(e) ce,(n) + D, Gek,t(e) sem) (17), 


Nn=0 
where C,, D, are arbitrary real constants, and DEU 


3. The Boundary Condition 


It is now required to determine C, and D, such that » is constant on the contour of the 
aerofoil section, where Sa) Firstly yo given by eqs. (3) and (8) must be expressed in terms of 
& and. This is achieved by deducing x and y as the real and imaginary parts respectively of z 


; 


MEERE 


\ 


NENNEN 


re “ FEN A 
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in eq. (11), and substituting in the expression for y,. Itisfound possible, using the stated approxi- 
mations, to write y, in the form 


; = Ffunfe) net (18), 
where 9, — O, and 

Mm=—V dt ausina + meosa + Fu [de + 31 eosh 25 — cos 20))], ] 
= —V fon+ 2 + zu cosh&sinx — 22, e*sin« 

+ 1 (cosh 3E£sina — cosh & sin 3 o)), 
= V(&sina +y,coso)e?:+ 5 qVd(cosh 2 Ec0s2x—]), 
Ha Savh (cosh 3& sin 3 — 3 cosh & sin a), (19). 
y -— ven u. u sinh &Ecosa + 22, 7° c05& 

+ a (sinh 3 & cos& — sinh &E cos 30)| £ 
»= V(ysina—ı, cosa) e”?® 5aVdsinh 24 sin 2&, 
DB Savh (sinh 3& cos3 x — 3 sinh & cos o) | 


When & = 0, it is known that 9, + y, = a constant, and if u,(0) is absorbed into the constant, 
and the new constant called p, then 


'T f1,(0) cosnn + »,(0)sinnn} + y C„Fek,(0) ce.) + D, Gek,(0) se,)} =.Pp.. #020): 
n=1 n=0 


In order to solve for C, and D, it is necessary to introduce normalised orthogonality con- 
ditions, viz. 


fi en) din = fs Endy-n, 
for all m, and 
f eo) cn) dn = fs en) sen) dn = 0, 
nm, and 
fi e,(n) sen.) dn = 0, 
for all n and m. 


Using the series expressions for the ordinary Mathieu functions, (6, p. 21, 22), it can be 
shown that 


2n 2n 
f ce,(n) sinmndn = [se,n)cosmndyn=ß®, 
0 0 


for alln and m, and 


2n 
(i) S ceon(n) cos mn dn — (), if mis odd, 
0 
ER Dr ANNEm='0, 
m 
Bauen h e 
Hl) aA it mi is.even. 
ar Ä 
(ii) f ceanzı(n) cos mn dn = 0, if mis even, 
0 


m—1 


+ - 1 © 0 
B-r31, if misiodd. 


unyC ( N 
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2 z . “ ” 
(iii) [ seon+ı(y) sin mn dn = 0, if mis even, 
a m—i 


= a(-1)" ® Adntı, ifmis odd. 


} 2n r R 
(iv) f sean+2(n) sin mn dn = 0,if m is odd, 
0 


m—2 


n+ UE: 
=n.(—1) * Bert iimiseven; 


By multiplying eq. (20) by the appropriate cosine or sine function, and by using the above 
integral formulae, the following results can be deduced: 


(a) Con Fekg„(0) = 2p(— 1)" AS" + 1)" AZ” u(0), 
v1 h 
(b) Can+ı Fekgn+1(0) = & (1)rtrHl Bit Bar+ı(Ü), En 
a . 
(c) Dan+ı Geken+ı(0) = PA ai Adrıı ver+ı(0), 
= 
(d) Den+2 Gekgn+2(0) = (—1)"+! B3”+? 9,(0) 


These equations define the constants C, and D,. 


4. The Tangential Fluid Velocity 


In view of the fact that the aerdynamic characteristics investigated in this paper depend 
only on the value of the tangential fluid velocity at points on the contour c of the aerofail, the 
actual expression for the stream function y is not determined. The main interest concentrates on 
the value of Oy/d& when & = 0, since the tangential velocity on the aerofoil contour is given by 
oy/H d£. From egs. (17) and (18) it follows that 

n=3 > > r ’ 
&), = 21 un(0) cosnn + »z(0) sinn} + 2 {Cn Fek,„(0) ce,(n) + D, Gek„(0) se,(m)} - - 22), 
with 9 = D,—=0. u 

In all the formulae which follow only first order quantities in g, x,, and y, are retained. 
It..is? known ıtthat A327 NAD BE (nu >00) are Ok is small enough, 
hence for given n to the order prescribed in the present analysis only terms withr=n,n+1 


Fek,„(0 
are retained!). In conjunction with these results, it can be shown, (6, p. 382), that F, = For 
, Gek,(0) . n 
is O(1/log g) when n = 0, and is O(1) when n >1, also G, = Gek,(ö) is O(1) for n >11, 


It will be shown later in this section that the constant p occurring in eq. (21) is O(log g), 
and if this is assumed here, it is possible to show to the prescribed approximation with 
1 > |glogg| >g, that 


h) n=3 * 
+ Ze Henn. Wa 
where =&-+in,&=0, and 

% = 2 p(Ap)? Fo; 

20 = u — MB? FR — iv, +iv(Ah2G,, 

20, = m —2pAABFR— AEQpA + mAD FR —in + iv (BI G,, 

203 = in BIBLF, + Bit, Bi 1 BI) F, — in, 


In evaluating the coefficients it is assumed that & = 0, and it.is then noticed that &%, and 
&, are O(l), &, is O(qlog g), and a, is O(g). 


(24). 


Hoe’ Hm 


The velocity at any point in the field of flow has components ( a = je and 
these are infinite when H = 0. This occurs at points where = —0(, and it is known that 


the only zero of dz/d? in the field of flow is at the point on the contour of the aerofoil referred 
to as the trailing edge. The corresponding point in the Z-plane is Z,, given by eg. (10). 


!) The author is of the opinion that the statement made by MeLachlan, (6, p. 46) should be restated as 
in the text. 
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The tangential fluid veloeity on c is given by Oy/H 0&, and this is finite when & = £, if 
Oy/o£ has a zero at this value of £. The condition for finite fluid velocity at the trailing edge 
of the aerofoil is thus from eq. (23) given by 


a +2 Io, (1 r IR) I ( = nn ) nee A (5). 


The eirculation round the aerofoil is defined by the integral 


. [| Op 
erde 

| ; (ä 
which from eq. (23) has the value 2a,. When q— 0, i.e. when the flow is irrotational, the 
circulation is a finite constant, hence the constant a, given by eq. (24), must be O(1). But F, 
is O(1/log q), hence it is necessary that p be O(log g), the assumption introduced earlier in this 
section. 

It is convenient here to introduce the values of the coefficients appearing in eq. (24). 

For small enough g, it follows by a method used in deriving formulae in (6, pp. 46, 47) that 


1 
1 
A}= bi=-<q. Again referring to (6, p. 248) it is possible to show that to the usual order 


8 
of small quantities 


1 1 
Bee len), 


where o=1+ 5 qlog q, and y is the Euler constant, such that y— log2 = — 0.1159 - - -. 
Also from (6, p. 382) it is known that to the given order F, = —— 
form il, \ £ ps g sanderı re on 


If p is eliminated from «, in eq. (24) by use of the expression for x, and known values 
of the other coefficients are introduced then 


2% =qVd(ei« —1)+ ee] te ya re (26). 
2 


This equation for a, is substituted into eq. (25), and the resulting expression solved for a,, leading 
to 


2h d? 1 ; 1 : 
= ren 2gh(4 +9 —1og2)} sina + zqn V sin 30° . (27). 
— qVd(cos2&x — 1) —2V y,cos«& 
If eq. (27) is substituted into eq. (26), to the usual order there results 
1— h 
20, = hVl SP _ gleina + g Vale) WERE: (28). 


From egq. (24) it can also be proved that 


> Beuel: 
2, = 2 V —)+g@—logdlsina+ivian+t Zur ziVgheis (9), 


and 


2 = zigqVhleis —eie) ERGO 
5. The Forces and Couple 
(i) For rotational flow the pressure equation, (5, p. 110), for a fluid of density o is 
ea a (u + v2) + fly) + constant . . 2.2 202.0 338), 


where f(y) is a function of the stream function y when the aerofoil is present. If F, F, are 
the forces acting on the aerofoil in the directions of the x-, y-axes respectively, then 


F\+iFr=/P«, 


. 2. . Math. Mech. 
368 Jones, The Forces on a Thin Aerofoil in Slightly Parabolic Shear Flow „4.37 Nr. 9/10 Sept. [Okt.1957 


where z= x +iy, and cis the contour of the aerofoil described in an anticlockwise direction. _ 
The contour c is a streamline, then f(y) is a constant on c, and since 


[E=0, 


then 1 
F,+IF, -—Io| + v2). dz. 


But for the fluid velocity measured on c, it follows that 


op \” 
2 Dee, 
u® + D ne: 
and since &=0 on c, then 
Pa An ee. ag 
de . dd 


Hence in all the force components are given by 
2arn 


1 Oy\” (dz\" 
B+ir= ie] (%) (®) dn ee er ee (32). 


If the lift force in the direction of the Y-axis is L, and the drag force in the direction of the 
negative X-axis is D, then 
L>=iD=ei(B ti es 20 ee 

The line integral of eq. (32) can be evaluated by considering the contour integral 

1 | oy\’/de\" _, 

dd, 

ee * 
where ©&=&-+in, taken round the rectangle R with corners at (0,0), (0,2), (&, 2), (&, 0) 
in the £-plane. If &£ >0, and sufficiently large, then (dz/dö)”" can be expanded in powers of 
exp (— £), by use of eq. (11). The integrand of the above integral is analytic everywhere inside 


and on the contour R, also the integral along the side „ = 0 and n = 2 cancel each other, 
hence it follows that the lift and drag forces can be deduced from the equation 


eu ie Oy\” dz\" 
L-iD=Sge /&) (2) den ir RE 


where c’ is the straight line joining (£, 0) and (&,2r), with & >0. From eg. (11), 


de\ IT 22 I 
en il BE E08 
(2) h ( h ). Te la 


mr ZT. X, ge 
en EM u Aero an et Me. —2ınd 
tortherfirst order in 2,,0,. 
From eq. (23) the following result can be obtained 


oy\” n—=4 u e R 
a)=” +2 (une ee ee 


where to the usual order of small quantities 
y=nt 20% +2, |] 
Yı = 2%, + 20105, | 
Ya = 2009 A AG ach 52 1 ee Se 
Y = 20903 + 20,0% , | 
Yyı=a + 20,0, J 

en y. are O(1), y5 is O(qlog g), and y, is Of(glogq%. 

serting the expressions for (dz/d£)! and (Aw/d&)? into eq. (34), the only contribution 


> the aus of the integral is the constant term in the integrand. Evaluation of this leads to 
1e result | 


L—iD=1i200u oe Zone (37), 
where only the first order terms have been retained. 


£ 
4 
E 
E 
£ 
E 
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The values of y, and y, are deduced from eq. (27), (28), (29) and (30) in conjunction with 
eq. (36), and then from eq. (37), )» (29) (30) i ] i 


L se 1 T d2\) . 


- (88), 
—y,41h@—y +log2)sin3a +2, cosa 
and 
D : : 
a nerEaloe 2) np ee) 


(ii) If the moment of the anticlockwise couple on he aerofoil taken about the origin 0 of 
the coordinates (x, y) is /', then 


PR ARE Wr 2idz; 


where P is given by eq. (31). 
On c it is known that /(y) is constant, and since 


Re |z&= 4 [aan =0, 


1. fay\t (dat 
nee (ae)2(z) EN an: 


The integral in eq. (40) can be evaluated as in (i) along the line & = constant, O<n< 2n. 
On this line 


then as in (i), 


— © (6%) 
z GE er ea re en a Lubee & (2n — 1) er+dE — ei, DE er 
de n—ı h h n=1 n=1 


hence using eq. (35), the constant term in the integrand, leads to 
TZRe2 zehn! Be 33 2 : nl. 


By use of eqs. (27), — , (30), the moment of the couple reduces to 


2: .! 1 4.d? 
Er DR si 
ae h e + zit 49 +log16 + fE Jin 24 


1 . (41). 
—„gh’sin do +2gdhsinasin2«— hx,sin2&«—2hy,sin®a 


(ii) The chord of the aerofoilistheliney=0, and the centre of pressure is the point where 
the line of resultant force intersects this chord. 
The line of action of the resultant force is given by 


N le 


and this meets the chord at the centre of pressure P, (X>, Y,), where y, = 0, and y = FE 
y 
— = The distance of P.from the trailing edge is thus 


- (Leos& + Dsin.o) ' 
2(h—%) + Xp; 


or 


134, +. zaleggqtyz g1)+ 20dsina 30 —2ukin 2 lest: 


6. Conelusions 


The main interest is concentrated on the perturbation effect of the curvature of the velo- 
city profile of the flow on the aerodynamic characteristics of the aerofoil, and to this end only 
terms involving gq are investigated, as these will give the variation effect over and above that 
due to the thin cambered aerofoil in a uniform stream. For small g, the variations of the lift 
and couple increase with increases in g,& or d, the drag variation correspondingly decreasing. 
Although the drag becomes negative, this is not an unusual state of affairs for perfect fluids, 

24 


not been accounted for here. For 


_teristies will be functions of the eirculation round the cylinder. 


eo 


distance of the centre of pressure from th ling 
to a minimum, then increases to positive values. It remains ı 
d = 0, and change in d with « = 0, however an increase of « with d 7 
of the distance, and so does an increase of dwitha#0. 

With a slight modification in the analysis the equivalent r 
in the shear flow can be deduced, although no Joukowski condition exist 
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Axiale Rückwirkungen von Überschallschaufelgittern 
Von I. Ryhming, SAAB, Linköping/Schweden ö 
(Aus dem Institut für Theoretische Gasdynamik der Deutschen Versuchsanstalt für Luftfahrt e. V., Aachen) 


Bei Anströmung rasch rotierender Schaufelkränze mit Unterschallgeschwindigkeit.kommt es zu Überschall- 
strömungen, welche mit Rückwirkungen stromaufwärts verbunden sind. Solche Strömungen werden für stationäre 
Fälle mit linearen und exakten Methoden behandelt. Bei Gittern ohne Leitapparat erweist sich die M achzahl 
der Anströmung als abhängig von der Umlaufgeschwindigkeit. Gitter mit Leitapparat zeigen verschiedene - 
Strömungstypen. Bei einigen von diesen unterscheidet sich die exakte Lösung von der linearen selbst bei 
kleinsten Störungen entscheidend. | 


In the case of rapidly rotating bucket rings a subsonic afflux may lead to supersonic flows or currents reacting 
on the incident flow. The stationary case of such currents is here treated both in the linear approximation and 42 
with exact methods. For cascades without guide vanes the Mach ratio of the afflux velocity to the velocity of 
sound turns out to be dependent on the rotational speed. Cascades with guide vanes show different types of flow. 


For some of these the exact solution differs essentially from the linear approximation even at the smallest per- 
turbation. 


Lors du soufflage de couronnes & palettes en rotation rapide avec vitesse subsonique des courants super- 
soniques se presentent etant lies ü des reactions en amont. Pour des cas stationnaires de tels courants son, 
traites a Vaide de methodes lineaires et exactes. Quand il s’agit de grilles sans distributeur le soufflage g’avert 
comme dependant de la vitesse de rotation. Des grilles avec distributeur presentent des types differents de courante 


Dans le cas de quelques-unes de celles-la la solution exacte se differencie distinctement de la solution lineaire. 
meme quand il y a des perturbations minimes. 


IIlpı TeyeHnn, HOCTyIalImmM Ha CKOPO BpamaroINmech NUCKM TYPOHH CO CKOPOCTLIO HHKE 
3BYKOBON, O6Pa3ylWTCA CBePX3ByKOBble TEeYeHUA, CBABAHHbIE C OÖPATHLIM BOBMeÜCTBHEM, 
HAIPaBJICHHEIM IIPOTUB TeyeHuun. Takme TeyYeHNsA B CTAUMOHAPHEIX CAYYAaAX 06PabaTbIBamTcH 
JIAHCÄHBIMH U TOYHBIMM MeTolaMmm. Y JIONACTHBIX PelIöToR 6e3 HanpaBısnmımero anmapara 
OKAa3bIBaeTcH, YTO uucıo Maxa IA TeYeHNA 3aBHCHT OT CKOPocTu Bpamenun. Pemeörkmm c 
HaAIpaBJIAIOINMM AlmapaToM BbI3bIBAIOT PA3HbIe TUIIBI TeYeHNMÜ. V HEKOTOPEIX U3 HUX TOYHOE 
PellleHme AaKe IIPH MallehllmX BO3MYIHCHNAX CYINECTBEHHO OTKIIOHAIOTCH OT IIMHEHHOTO. - 


1. Einleitung 


Der durch den Titel abgegrenzte Problemkreis wurde zuerst von K. Oswatitsch in einem 
Vortrag auf der V.D.1.-Fachsitzung in Karlsruhe 1953 behandelt. Auf die seinerzeitigen Resultate 
wird im Text hingewiesen. Einen weiteren Kurzvortrag über den derzeitigen Stand der Arbeiten 
hielt der Verfasser auf dem 2. Europäischen Luftfahrt-Kongreß in Scheveningen 1956. Die 
ziemlich umfangreichen Ergebnisse werden in zwei Veröffentlichungen mitgeteilt. Die vor- 
liegende behandelt vorwiegend stationäre und quasistationäre Fälle, während eine weitere Arbeit 
auf die wesentlich instationären Erscheinungen eingehen wird. 

Den folgenden Betrachtungen liegt ein in üblicher Weise abgewickeltes, gestaffeltes Flügel- 
gitter einer Turbine oder eines Axialverdichters zugrunde. Dieses werde unter einem zunächst 
nicht näher festgelegten Winkel mit Überschallgeschwindigkeit angeströmt. Die Machzahl M 
der Anströmungkann dann in eine Normalkomponente zum Schaufelgitter M, (Axialkomponente) 
und eine Tangentialkomponente M, (Umfangskomponente) zerlegt werden. Je nachdem, ob die 


BR 
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 Kanalströmung mit Charakteristiken-Methoden 
berechnet werden. Der erstere Fall der axialen 
"Rückwirkung (M„< 1) ist problematischer und 


nicht unnötig zu komplizieren, seien dieSchaufeln 
als so tief angenommen, daß keine Machlinie 


en den een rs =. eine ebene 


sei im folgenden behandelt. Um das Problem 


vom Schaufelende vor das Schaufelgitter gelangen 
kann, wie es in Bild 1 dargestellt ist. In diesem 
hier ausgeschlossenen Falle gibt es Rückwirkun- € 
gen durch das Schaufelgitter hindurch. Er kann 
nicht gelöst werden ohne nähere Angaben über 
die Strömungsbedingungen hinter dem Gitter [1]. 


Das formulierte Problem interessiert vor 
allem für das Studium der Strömungsverhältnisse 
am Eintritt oder Austritt von Laufrädern hoher 
Umlaufgeschwindigkeit, beispielsweise bei Cur- 
tisturbinen. In der Praxis kommen noch Leit- 
räder oder Lavaldüsen hinzu, was von ausschlag- 
gebender Bedeutung ist, wie gezeigt werden wird. Mm 
Die Strömung sei als quasistationär ange- 
nommen, d.h., die Strömung wird für jede Stel- 
lung des Laufrades zum Leitrad als stationär be- 


handelt. Diese vereinfachte Annahme hat für eine M,>1 Me 

erste Behandlung den großen Vorteil, schnell all- i  Anströmung 

gemeine Ergebnisse zu liefern als Ausgangspunkt er 

für genauere und speziellere Rechnungen. Bild 1. Flügelgitter ohne Leitrad bei verschiedenen Machzahlen 


2. Gitter ohne Leitrad 

Das Gitter sei hier wie im folgenden als Strichgitter paralleler Geraden und das Koordi- 
natensystem sei in diesem Abschnitt als ruhend bezüglich des Gitters angenommen. Die Normal- 
komponente W, der Anströmung ist dann dieselbe wie in dem bezüglich des Einlaufes ruhenden 
System. Die Tangentialkomponente W, der Anströmung im gitterfesten System ist gleich jener 
im Einlauf festen Systems, vermindert um die Umlaufgesch windigkeit des Gitters. Im „Normal- 
falle‘ weisen die Gittergeraden in Richtung des Vektors der Anströmgeschwindigkeit im gitter- 
festen System und vom Gitter gehen keine Störungen aus. In dem hier zu behandelnden Falle 
axialer Rückwirkungen jedoch unterscheidet sich der Gitterwinkel + e vom Winkel des 
Geschwindigkeitsvektors der Anströmung: 


WW teh Le re Pat, 


wobei e stets als klein vorausgesetzt sei. Es ist dabei grundsätzlich dasselbe, ob man eine Un- 
genauigkeit im Gitterwinkel, in der Anströmrichtung des Einlaufes oder in der Umlaufgeschwin- 
digkeit des Laufrades annimmt. Wesentlich ist nur, daß die Bedingung für den Normalfall nicht 
genau erfüllt ist. 

. Das kartesische Koordinatensystem sei so gedreht, daß die x-Achse in Richtung der An- 
strömung und die y-Achse normal dazu weist. Der Gitterwinkel erscheint in diesem System als 
Winkel e. 

Für die Behandlung der Rückwirkung eines Laufrades ohne Leitrad sei zunächst eine dre- 
hungsfreie Strömung angenommen, die jedoch „schwache“ Stöße aufweisen darf. Darunter sind 
Stöße zu verstehen, bei welchen die Stoßpolare durch die Charakteristik im Hodographen er- 
setzt werden darf. Es ist bekannt, daß dies bei schiefen Stößen, welche nicht auf Unterschall- 
geschwindigkeit führen, in einem weiten Bereich erlaubt ist. Unter dieser Voraussetzung gelten 
auf linksläufigen und rechtsläufigen Machlinien (Bild 2) folgende zwei Bedingungen: 


EEE RE RE SIEr ER Le) 


Hierin ist 9 der Strömungswinkel und Ch(W) eine monotone Funktion der Geschwindigkeit, wie 


sie etwa durch die Prandtl-Meyer-Expansion gegeben ist ([2], S. 306, Gl. (83)). 
In Bild 2 ist die Lage der Punkte 3 und 4 dadurch gekennzeichnet, daß von ihnen eine 


linksläufige Machlinie zu den Punkten 1 und 2 der unteren Gittergeraden und eine rechtsläufige 


24* 


meer A nichts Neues, Die er 
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ini | j ührt. Daß durch diese Annahme 
Machlinie zu denselben Punkten der oberen Gittergeraden fü r 
die Base der Punkte nur indirekt festgelegt ist, tut für das Folgende nichts zur Sache. Gemäß 
Gl. (2) gilt beispielsweise für Punkt 3undl: 
Ch(W;) — 9; = Ch(W,) — ds; 
Ch(W;) + 9 = Ch(W,) + %ı- 
— il die F i ber monoton ist 
D folst: 9, = 9, und Ch(W,) = Ch(W,). Weil die Funktion Ch(W) a s 
folgt eh W EW. Derselbe Beweis kann für Punkt 4 und 2 und für alle auf den Verbindungs- 
\ linien 34 und 1—2 gelegenen Punkte geführt werden. 
D.h., auf der Linie 3—4 herrscht derselbe Konstante 
Strömungswinkel wie auf der Linie 1—2 der Gitter- 
geraden. Da die Linie 3—4 aber auf rechtsläufigen 
Machlinien aus dem Anströmgebiet mit dem konstan- 
ten Zustand W, 9, erreichbar ist, muß auf 3—4 wie auf 
1—2 die gleiche konstante Geschwindigkeit herrschen. 
Sie ist wegen 9, = 9%, = 9, = d, gegeben durch 


Diese Überlegung gilt auch dann, wenn zunächst 
ein von der Vorderkante ausgehender Stoß gemäß 
Bild 2 angenommen wird. Wenn aber beiderseits des 
Stoßes der gleiche Strömungszustand herrscht, so be- 

Bild 2. Charakteristiken und Stöße im Gitter steht auch keine Ursache mehr für einen Stoß und es ist 

nicht einzusehen, wie Verdichtungsstöße in der Strö- 
mung entstehen sollen, wenn sie nicht durch die Vorderkanten der Gittergeraden ausgelöst 
werden. ; 

In einer stoßlosen Strömung jedoch geht von jedem Punkte der Gitteroberseite, welcher 
vor der zur darübergelegenen Gitterkante führenden Grenz-Machlinie liegt, sowohl eine links- 
läufige als auch eine rechtsläufige Machlinie ins Anströmgebiet. Damit gilt für jeden solchen 
Punkt der Gitteroberseite: 


Ch(W) — #-= Ch(W) — 9% , 
Ch(W) + 8 = Ch(W) + dx ) 
te ee 


We Was 


d. h., der Anströmzustand dreht sich unter der Einwirkung des Gitters in dessen Richtung; das 
Resultat ist eine ungestörte Parallelelströmung mit W„ und 9% im gitterfesten System. Bei 
M„< 1 kann folglich nicht Wx» und 9 vorgeschrieben werden, wie bei einem hyperbolischen 
Anfangswertproblem, sondern nur eine der beiden Größen oder eine Kombination von ihnen. 


Zum gleichen Schluß aus physikalischen Überlegungen kam Kantrowitz in einer Geheim- 
arbeit aus dem Jahre 1946, die im Jahre 1950 veröffentlicht wurde [3]. Unabhängig von dieser 
Arbeit hat K. Oswatitsch dieses Ergebnis für die linearisierten Gleichungen in Karlsruhe mit- 
geteilt. Mit Rücksicht auf die hier gewählte vorwiegend mathematische Behandlung und die 
über das Resultat geäußerten Zweifel erschien eine nochmalige eingehende Behandlung des 
Problems gerechtfertigt. Siehe auch Schwaar [4]. 


Vor dem Laufrad befindet sich ein Einlauf, der stromaufwärts in die freie Atmosphäre 
mündet. Die Strömung werde nun von dem mit dem Einlauf festen Koordinatensystem betrachtet. 
Instationäre Störungen, welche das Laufrad etwa durch eine plötzliche Änderung der Umlauf- 
geschwindigkeit aussendet, pflanzen sich stromaufwärts fort und verlieren sich in der freien 
Atmosphäre. In dieser Betrachtungsweise hat der Einlauf‘ wegen M„< 1 Unterschallströmung. 
Jene Gasteile, welche aus der Atmosphäre in das Laufrad gelangen, können von diesem keine 
wesentliche Drehung aufgeprägt bekommen und strömen folglich wirbelfrei ein. Das bedeutet, 
daß W, im System des Einlaufes Null sein muß und das Laufrad gerade eine solche Normal- 


komponente W, erzeugt, daß 9, im System des Laufrades mit dem Winkel auf der Laufschaufel 
übereinstimmt. 


3. Gitter mit Leitrad, lineare Theorie 


In diesem Abschnitt soll die quasistationäre Strömung in einem Gitter mit Leitrad behandelt 
werden. Für die Behandlung werden, entsprechend der Voraussetzung kleiner Störungen, welche 
es De zugrunde liegt, die Gitter als Strichgitter paralleler Geraden angenommen 

e * 4 “ 1 “ . = . . > } . ” ” . ” ; 
(Bild 3). Die speziellen Gittereigenschaften oder die Gittergeometrie sei bezeichnet mit: Teilung 
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Im folgenden sollen bei der 
ersten Behandlung die Randbe- 
- dingungen dadurch erfüllt werden, 
daß auf den zur x-Achse parallelen 
Laufschaufeln Quellen und Senken 
angenommen werden, was einer 


_ Randbedingung verwenden, wasdas 


on der Hinterkante der 
ezeichnen; ferner definieren wir als ‚, 


eiLeitschau- 
‘von der Grenz-Mach- 
der Verbindungslinie zwi- 


leitrad. 


= 


Erz ren 
ZIERT 


vorgegebenen v-Komponente 
dort entspricht. Bei der nächsten 
Behandlung wollen wir die exakte 


Bild 3. Koordinatensystem im Gitter mit Leitrad. 


Verschwinden der Normalkompo- 
nente der Geschwindigkeit auf der unter e geneigten Laufschaufeln bedeutet. Es wird sich zeigen 
daß zwischen der ersten, genäherten und der zweiten, exakten Erfüllung der Randbedingungen 
ein wesentlicher Unterschied besteht als beim Umströmungsproblem eines Einzelprofiles. 


4. Genäherte Randbedingungen 


Innerhalb der Linearisierung kann eine Änderung der v-Komponente als Richtungsänderung 
der Strömung, eine Änderung der u-Komponente als Druckänderung ausgesprochen werden. Auf 
den beiden Scharen der Machlinien, in Abb. 3 als Geraden unter x = + 45° eingezeichnet, gelten 
dann die Bedingungen: 


auf& = konst.: u—v = konst.; | 


aufn = konst.: u+ v= konst. f 0) 


wobei u und v die Komponenten der Geschwindigkeit W in den Achsenrichtungen sind. Durch 


dieses hier nicht abgeleitete Ergebnis der Theorie ebener linearer Überschallströmung liegt die 
Lösung bereits völlig fest. 

Die Strömungen ohne und 
mit Leitrad unterscheiden sich in- 
sofern stark, als die axialen Rück- 
wirkungen im letzteren Falle auf 
ein kleines Gebiet beschränkt sind. 
Bild 4 zeigt ein Leitrad- oder eine 
Reihe von Lavaldüsen — mit glei- 
chem Gitterabstand wie das Lauf 
rad. Natürlich sind auch andere 
Teilungen möglich, aber um die 
Übersichtlichkeit der Arbeit nicht 
zu stören, wird die Behandlung 
durchgehend auf diesen Fall be- Bild 4. Normalfall der Strömung im Gitter. 
schränkt. Die Störungen des Leit- 
rades können sich stromaufwärts von jener Machlinie nicht mehr bemerkbar machen, welche in 
das Leitrad, ausgehend von seinen Endkanten, hineinläuft. Davor herrscht ungestörte Anströmung 
Lg d%, = 0. Der mit 1 bezeichnete Teil des Laufrades steht nur unter dem Einfluß des Anström- 
zustandes; u, kann mit der unteren Gl. (4) sofort berechnet werden. Der Teil2 des Leitrades 
steht hingegen nur unter dem Einfluß von Teil 1 und ist mit der oberen Gl. (4) zu berechnen. 
Teil 3 steht wieder nur unter dem Einfluß von Teil 2 usw. So fortschreitend läßt sich in dem in 
Bild 4 gezeigten „Normalfall“ die ganze Strömung berechnen. Das einfachste lehrreiche Bei- 
spiel ist wieder ein Geraden-Strichgitter mit einer konstanten Quellergiebigkeit v, im Laufrad. 
Gemäß Gl. (4) gilt dann: 


mel,--b,, g=UuU—b, a BE ende 


Bei jeder Reflexion einer Machlinie, 
Leitrad führt, wird der u-Wert um den »tra 
positiv sein, je nachdem, ob Senken oder Quellen angenommen sind). in 
 Jinie, welche nur von einer Laufradwand zur anderen führt, bt beim 
| ine u-Änderung. cu vr Er 
ee De een Normalfall ist die einfachste Form der axialen Rückwirk 
Gitter. Die oben angegebene Lösung gilt für a = 45°, ß = 30°, e = 0 innerhalb der x T 


1 >2 >. me ao 
a 


(siehe etwa Bild 6, Formel 1: 
An der Grenze b/a = 0,536 tritt 
nämlich die von der Hinter- 
kante einer Leitschaufel aus-- 
gehende linksläufige Machlinie, 
welche einmal an der Unterseite 
der ‚benachbarten Leitschaufel 
und einmal an der Oberseite der 
Laufschaufel reflektiert wird, 
vor die Vorderkante der näch- 
sten Laufschaufel, siehe Bild 6, 
I,1. Damit fängt das Strö- 
mungsbild an, sich zu ändern. 
In Bild5 ist das Verhältnis 


5 e FE en. ; N. r Br b 1 1 
Bild 5. Strömungsbild bei 2a = 0,425; beachte die Einwirkung der Grenz-Machlinie! et 0, 495 gewählt. Die oben 


erwähnte Machlinie wird dann erst nach einigen Reflexionen an folgenden Leit- und Laufschau- 
feln schließlich in das Laufrad hineinreflektiert. Zur Berechnung der Strömung in Bild 5 geht 
man vom Zustand 0 aus und gelangt zu den Zuständen 1, 2und3. Um nun im Gitter nicht von 
einer Schaufel zur nächsten fortschreiten zu müssen, geht man am besten längs der gestrichelten 
Linien wieder zurück auf die erste Laufschaufel. Von dort geht es über 4 nach 5 usw. 

Wegen =, =», =D, = (Qunde, =v =D, = 9, konst. alt Tmun 


een Sr St SR 
Erst von 9 aus wird die Machlinie in das Laufrad hineinreflektiert. 


Man sieht aus dem Vorhergehenden, daß der Machlinienzug, welcher die Fortsetzung der 
Grenzlinie darstellt, eine große Rolle für die Strömung spielt; denn hinter jedem reflektiertenTeil 
dieser Machlinie, welcher vor den Vorderkanten der Laufschaufeln auftritt, ändert sich der Zu- 
stand. Man findet nun, falls man die Laufschaufeln zu noch kleineren Verhältnissen b/a verschiebt, 
daß auch noch zwei andere Machlinien vor die Vorderkanten der Laufschaufeln treten und also 
auch sehr stark auf die Strömung einwirken können. Diese zwei anderen Machlinien entstammen 
der von der Hinterkante einer Leitschaufel ausgehenden rechtsläufigen Machlinie und der von 
der Vorderkante einer Laufschaufel ausgehenden linksläufigen Machlinie (Bild 6: III, 1; II, r; 
Um nun diese Verhältnisse im Gitter besser überblicken zu können, wollen wir uns zunächst mit 
der „Reflexionsgeometrie‘‘ der betreffenden Machlinien im Gitter befassen. Dabei wollen wir für 
jede der drei oben genannten Machlinien untersuchen, wieviele Male sie im Rückwirkungsgebiet 
bei verschiedenen Werten von b/a reflektiert werden, ehe sie in das Laufrad endgültig hinein- 
reflektiert werden. Wir wollen dann die Anzahl der linksläufigen Ma chlinien N bestimmen, 
welche vor der Vorderkante einer Laufschaufel abhängig von b/a auftreten und welche von irgend 
einer dieser besonderen Machlinien herrühren. Zuerst müssen wir dann solche Grenzstellungen 
finden, wo eine der oben genannten Machlinien gerade vor die Vorderkante der Laufschaufel tritt 
und auch wo sie endgültig und unmittelbar in das Laufrad hineinreflektiert wird. Eine Zu- 
sammenstellung solcher verschiedenen Stellungen b/a ist in Bild 6 gezeigt. Im folgenden soll 
die Behandlung auf die Daten x = 45°, ß# = 30° und e = 0° beschränkt bleiben. 


q 
R 
Y 
ß 


Von der Grenz-Machlinie ergibt sich nun für 1 > A > 0,536... nur eine linksläufige 
a 


Machlinie vor der Vorderkante einer Laufschaufel. Bei ie u a = 0,536... beginnt 
a sin(a+ß) 
die Anzahl zu wachsen. Nimmt man zunächst einen Wert = : ne u r - 6 an, so 
cos «& sin 
findet man, falls man ö klein genug wählt, daß die Anzahl N der Machlinien vor der Vorderkante 


< 
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Re. form ein = 


b_, Sina cose _sinla-B) 
8 sin(&+E)  sin(&+ß) 


Wert bei: | Wert bei: 
45% =30° 


&=15° 


cos&-sin(&-ß) tg + 19° 
cosß-sin(a+E) tga-(3tg B+rtg a)+tge-(tga-tgß) 


3188-0) 
a sin(a+ß) 


b_ı sinta-B) ._sin(&-e) 


a 2 cosasinß sin(&+E) 


sin(&-B) _sin(«-E) 
cos &sind sin& cos& 


sin(&-B) sın(a-e) 
sin(&+ß) sin (a+&) 


(tga-tgß) - (tga-tge) 
tgB (3tg a -tge)+tga(tga+tge) 


1 sın(® -B) sin(x-E)sin (&+E) 
2 sina-cos &-sin(ß+e)[sin(&-E)+sin(&+E)] 


Bild 6. Die einfachsten Reflexionseigenschaften der drei verschiedenen Machlinien 


1 sin —Pß) 
2 cosasinß 
ergibt sich, daß aus dieserMachlinie im ganzen Gitter ein einziger Machlinienzug entsteht, welcher 
fortlaufend von einer Teilung zur anderen ins Unendliche reflektiert wird (Bild 8a). Wir können 
einfach sagen, daß die Grenz-Machlinien bei dieser Stellung durch Reflexion unmittelbar in sich 
selbst übergehen, und daß der oben gegebene b/a-Wert ein Häufungspunkt für die betreffenden 


exakt, so 


der Laufschaufel sehr groß ist (Bild 6,und 7a). Wählt man - 


Machlinien ist. Für 0,366... > 7 > 0,308... gibt es wieder nur eineMachlinie vor der Vorderkante 
der Laufschaufel (Bild 6: I,2; 1,3). 
In dem Gebiet 0,308... > s2 > 0,268... ist das Reflexionsbild viel komplizierter. Man kann 
a 


in diesem Gebiet eine Reihe von Häufungspunkten finden, und es zeigt sich ferner, daß die 
Häufungspunkte dichter und dichter liegen, je mehr man sich dem b/a-Wert 0,268... nähert, d.h., 
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s b ? R 
die Häufungspunkte haben selbst einen Häufungspunkt bei dem Wert EZ 0,268.. Diese Ver 
hältnisse sind in Abb. 7 dargestellt. Er TR 

Führt man dieselben Untersuchungen auch für die beiden anderen Machlinienzüge durch, 
so gelangt man zu den Resultaten, die in Bild 7b und c dargestellt sind. 


A 
N 
= 
Bu 2 
b 
b — 
1 2 7 = 
7 
1) 0.268° 0,208 0,366. 0,536 1 0 TaBanZ HER 
‘ b 
N 
5 
b 
1 3 
0 0153 0,1830,268 1 


G 


Bild 7a—c. Anzahl N linksläufige Machlinien im Rückwirkungsgebiet für die drei verschiedenen Machlinienzüge. 


Es ist nun festgestellt worden, daß bei gewissen Grenzen die Anzahl bestimmter Machlinien 
unbegrenzt wächst. Man kann sich auch überzeugen, daß genau bei einem solchen Wert eine 
gewisseMachlinie unmittelbar oder nach einigen Teilungsperioden in sich selbst übergeht; dabei 
gibt es im Gitter zwei Stellungen von dem ersten Typus, nämlich bei 


| ‚sn (a —P) _ 9366 

a2 cosasiın yes 
und 

ur BETTER: 

a 4 cosasinß wer 


Diese beiden Stellungen sind in Bild 8a und c dargestellt. Beispiel für eine Stellung, wo eine 
Machlinie erst nach einigen Teilungsperioden in sich selbst übergeht, findet man in der Gegend 


b 
des Wertes m 0,268... In Bild 8 ist ein solches Beispiel gezeigt, wobei 
b 2 


= (ya 1292 
S le 1) = 0,292... 


Bei einem Studium der Bilder 8, a, c und d sieht man, daß es bei diesen Stellungen geschlossene 
periodische Kurvenzüge gibt, welche unverändert ins Unendliche laufen. In Bild 8a und eist 
die Periodizität dieser Kurvenzüge eins, in Bild 8d dagegen vier. Falls man nun näher und näher 


R b n E : 
dem Wert -- 0,268... Häufungspunkte aufsucht, wird man finden, daß die Periodizität der 
geschlossenen Kurvenzüge bei diesen Stellungen zunimmt. Bei dem Wert .n = 0, 20er 
a 


findet man schließlich den in Bild 8b wiedergegebenen Fall. Hier gehen nun alle geschlossenen 
Kurvenzüge in einen einzigen über. 


Zunächst wollen wir uns mit dem Strömungsproblem bei solchen Häufungspunkten befassen. 


. Das Auftreten von geschlossenen Kurvenzügen, längs deren man laufend zu höheren oder 
niedrigeren u-Werten gelangt, bedeutet einen Widerspruch. Diese Nichtexistenz von Lösungen 
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kann man auch so ausdrücken, daß 
die Bedingungen nach Gl. (4) zu 
einem linearen Gleichungssystem 
führen, dessen Koeffizientendeter- 
minante bei gleicher Anzahl von 
Gleichungen und Unbekannten 
verschwindet. In der unmittelbaren 
Umgebung von einem Häufungs- 
punkt bekommt man nach Gl. (5) 
mit i— oo ein kräftiges Aufschau- 
keln des Zustandes, ein Zeichen, 
daß in diesen Gebieten die lineari- 
sierte Theorie nicht länger gültig 
ist. : 


Es ist bekannt, daß die line- 
are Überschalltheorie immer große 
Fehler in großen Entfernungen von 
der Störstelle gibt. Hier beim Git- 
ter haben wir ein gutes Beispiel 
dafür bekommen, denn irgend eine 
Machlinie, welche in den schraf- 
fierten Gebieten im Bild8 auf- 
tritt, wird ins Unendliche reflek- 
tiert. 

Dies ist wohl der erste be- 
kanntgewordene Fall, für den es b 
keine Lösung für die linearisierte 
Überschallgleichung gibt. Hin- 
gegen ist bekannt, daß eine Lö- 
sung der exakten Gleichung an 
„Grenzlinien‘‘ aufhört, wie erst- 
malig Tollmien zeigte. Damit 
ist nun ein Beispiel gezeigt, das 
zu Widersprüchen führt. Ein Teil 
von diesen Ergebnissen wurde 
schon von K. Oswatitsch in 
Karlsruhe gebracht. 


Z,4 und IL,1 


d 


Bild 8a—d. Geschlossene Kurvenzüge im Gitter. Periode der Reflexionen = 1 in Bild 8a—c. Periode der Reflexionen = 4 in Bild 8d 
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5. Exakte Randbedingungen 2 
i it ei indem wir den Lauf- 

Nun wollen wir uns mit einem etwas anderen Problem befassen, in 
schaufeln einen Anstellwinkel e geben. Die neue exakte Randbedingung auf den Schaufeln ist 


egeben durch 
ses ige= + 


] em, ob wir eine Expansions- oder Kompressionsströmung berechnen wollen. Das neue 
ist vor allen Da in dem Sinne interessant, daß wir mit der Erfüllung der exakten 
Randbedingungen eine engere Verknüpfung mit der in Kapitel 7 und 8 folgenden exakten 
Behandlung des Problems erwarten können. Es wird sich zeigen, daß diese Behauptung richtig ist. 
Wie früher studieren wir zuerst die Reflexionsgeometrie der Machlinien im Gitter. Wir 
wollen dabei für unsere Untersuchung einen solchen Wert von & wählen, daß wir eine möglichst 
einfache Reflexionsgeometrie bekommen. Man kann dann zeigen, falls man e so wählt, daß 


cot ß — cot a 8>88°% bei redsh, B = 30° 
cotß +3 cota 


ige >tigoa: 


ist, daß der in Bild 7 gezeigte Kurvenastinnerhalb 0,308 ... > =: > 0,268... wegfällt. Wir legen 


daher im folgenden e = 15° fest (Expansionsströmung). In Bild 9 ist wieder N als Funktion 
von bja bei = 45°, = 30° und e= 15° für die drei Machlinien dargestellt. Es ist unmittelbar 


— 
| 


Sl 

S 
[3 IT 
ya 


[) 0,316 0,423 7 ö 


Bild 9a—e. Anzahl N linksläufige,Machlinien im Rückwirkungsgebiet bei Anstellung der Laufschaufeln für die 
drei verschiedenen Machlinienzüge 


aus Bild 9 zu sehen, daß wir nun ein viel einfacheres Problem bekommen haben. Im vorher- 
gehenden Kapitel gab es Stellungen, bei denen man für 


b b 
u 
a Ali 


einen unbegrenzten Zuwachs der Machlinien erhielt, wobei (b/a); einen Häufungspunkt bedeutet. 
Besonders beschwerlich war dabei die Umgebung des Punktes 2 — 0,268. Es ist uns nun aber 
a 


gelungen, durch eine geeignete Wahl des Winkels & diesen Wert 2 — 0,268, der von e nicht 
beeinflußt ist, zu vermeiden (siehe etwa Bild 6, Formel I : 4). i 


In Bild 10 ist die Lösung für das Strömungsproblem für bestimmte feststehende Stellungen 
innerhalb einer Teilung des Leitrades dargestellt (x = 45°, ß = 30°, & — 15°). Um die Strömung 
bei dem Normalfall Bild 3a zu berechnen, haben wir für den mit 1 bezeichneten Teil des Lauf- 
rades gemäß Gl. (4) (siehe etwa Bild 4): 
vn =u+D» 
D, 


ige = — — 
= I 


oder De 


Z. angew. Math. Mech. . - ; er 
Bd. 37 Nr. 9/10 Sept./Okt.1007 Pyhming, Axiale Rückwirkungen von Überschallschaufelgittern 379 


oO 


Bild 10a—d. Die lineare Lösung des Strömungsproblems innerhalb einer Teilungsperiode 
für « = 45°, ß = 30°, e = 15° (Expansionsströmung) 
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Die in Bild 10 verwendeten Zahlen bedeuten folgende u-Werte: 
1 1+tge\? 
1 = W' —o— 2: u=u were 
er ge 
1+tge+tg!e+tg?e 
Iertere BIER) 
ea en 3:u=w ae 
' 1+1g 8)? 
3 BEL Li ur, 1 en 
° (1—tg ee)? 8 . 
2 (1+tge 
need Br N 15: = go’ 
(1—1ge)’ € e b 
€ 
5 NE 23 lo au Be 
m dige a r 
Io eN\2 In ge)! 
> u), I u ee 
a 
10, 18; u = %: er B 
1 + tg.) 
1 i-ttee Dt 
8 uam ge 19-2 u, Teer 
6 
1 /1-+-tge\? = ir 
i gl) 20: Br ne 
1 2 Fre 
10: U dpa BA = Net 
\7 
1+tge+3tg?e—tg?e u ee 
ee NL ug 
Dieses Ergebnis gilt für die Machzahl y2; für eine beliebige Machzahl bekommt man: 
Ug 
a - 
A Ztagnge 


b ; : 
Die Normalfall-Lösung ist gültig innerhalb 1 > = >0,423... In dem Bereich 
0,423... > a2 > 0,316... läßt sich die Strömung ebenfalls genau so wie früher berechnen, wobei 
2 v 
man nur die Randbedingungen v = O0 an der Leitschaufel und tg e = — Fra der Laufschaufel 


beachtet. In Bild 10b ist ein Beispiel innerhalb dieses letzten Bereiches dargestellt. 

In Bild 10 sind die konstanten u-Geschwindigkeiten in den verschiedenen von den Mach- 
linien und zum Teil von der Geometrie abgegrenzten Feldern eingetragen. Es ist sofort mit 
Hilfe dieser u-Werte zu sehen, wo irgend ein Stoß in dem Strömungsfeld liegen muß; denn über- 
all hinter einem Stoß sind ja die W-Werte plötzlich kleiner geworden. Diese Machlinien, welche 
einem Stoß entsprechen, sind in Bild 10 dicker als die gewöhnlichen Machlinien gekennzeichnet. 


b s R 
Zunächst wollen wir uns nun mit der Stellung : 0,316... befassen. Bei der früheren Be- 


handlung gab es bei dieser Stellung einen Widerspruch. Mit den in Bild 10c angegebenen Be- 
zeichnungen erhält man ein lineares Gleichungssystem, dessen Lösung wie folgt lautet: 


A U, 4 rer ee > Zone 
uU = —— : U, = u —- ——; = mel: 
—tgie 1—tge 
Up= 5 Up; g=U,—=l. 


Es besteht also hier kein Widerspruch mehr, aber das Resultat u; = u, — 0 zeigt schon, daß die 
Linearisierung für diese Stellung nicht verwendbar ist. 
Es läßt sich nun zeigen, indem man die Strömung für verschiedene Stellungen innerhalb 


b ; R een. 
Malbnn> 7 > 9155. .. untersucht, daß im wesentlichen dasselbe Resultat wie für die Stellung 
b 
an 0,316... herauskommt. Damit steht fest, daß man in diesem ganzen b/a-Gebiet die Lineari- 


sierung nicht verwenden darf. In diesem Gebiet treten nämlich ähnlich geschlossene Kurven- 
züge auf, die sich zum Teil bis ins Unendliche fortsetzen. In Kapitel 7 soll gezeigt werden, wie die 
exakte Strömung-in diesem Gebiet aussieht. Wir können dann die exakte Lösung einfach physi- 
kalisch erklären mit Hilfe der Ergebnisse, die wir oben bekommen haben. In Bild 10e ist das 


Strömungsbild innerhalb 0,155... > 5 >0,134... dargestellt. Hier treten nun keine solchen Kur- 


venzüge wie oben auf, und man bekommt eine gültige Lösung: Für 0,134... > > 0,1082 


a 


. 


lossen »n Kurvenzüge auf, und man bekommt eine Lösung in 
fear N Bere je} DR E { BEE RR TEL A On A 
Bereich, welche der Lösung für a 0,316... ähnelt. Schließlich werden innerhalb E 
Be 108... 9 >0 alle Machlinien direkt in das Laufrad hineinreflektiert; die Lösung erhält ie 

er; Bert ? | i j m 
nun einen Normalfall-Charakter (Bild 10d). | | pe j 

Bo; Es ist zu bemerken, daß wir mit der linearen Theorie überall eine Lösung bekommen, falls u; 


wir den Laufschaufeln einen kleinen Anstellwinkel geben und die darauf exakte Randbedingung 


E- ‚erfüllen. Ein kleiner Winkel kann also den früher erhaltenen Widerspruch aufheben. 5x 


IF 6. Gitter mit Leitrad, exakte Theorie » 
Das früher für kleine Störungen behandelte Strömungsproblem in Gittern soll nun hier 
mit Rücksicht auf die Endlichkeit der Störungen, d. h. ohne die Linearisierung, erweitert werden. £ 
Die Darstellung ist jedoch wieder auf ebene Strichgitter beschränkt. 2 
Das Problem soll mit der Charakteristiken-Methode, und zwar nach der drehungsfreien 
Näherung, stationär zweidimensional behandelt werden. Um diese Methode verwenden zu kön- i “ 
nen, sollen die auftretenden Stöße von einer geringen Stärke sein. Dieses ist eine Bedingung, die 
bei unserer Behandlung immer erfüllt ist, weil der Anstellwinkel der Laufschaufel hier klein 
angenommen ist (Bild 3). 
Bei der Konstruktion des Strömungsbildes ist eine geeignete Gitterstellung herausgegriffen 
worden, d. h., die relative Stellung des Laufrades zum Leitrad ist von Fall zu Fall vorgeschrieben. 
Man muß nun zwei verschiedene Fälle bei der Strömung durch das Gitter unterscheiden: 
Gitter, in welchen das strömende Medium eine Expansion erfährt und Gitter, in welchen das 
Medium komprimiert wird. 


a PN NEN 


7. Das Expansionsproblem 

Bei der Behandlung der Expansionsströmung wählen wir einen praktischen Fall mit einem 
Gitterwinkel # = 20° und einem Anstellwinkel e = 5° und studieren nun die Strömung durch das 
Gitter bei verschiedenen Stellungen. Die Anström-Machzahl legen wir mit M = 1,435 fest. Diese 
Machzahl ist wegen der Gittergeometrie konstant bis zu der vom Ende der Leitschaufel aus- 
gehenden Grenz-Machlinie. 

Zuerst wird nun eine solche Gitterstellung ausgewählt, daß das Verhältnis b/a nahe eins ist. 
Man kann dann vermuten, daß die Strömung nicht so sehr von der reinen Kanalströmung ab- 


weicht, welche für " — ] entsteht. 


Mit einer rechtslaufenden Charakteristik (entsprechend Bild 11) von dem ungestörten 
Anströmungsgebiet in dem Leitrad gelangt man auf der Laufschaufel zu einem Zustand, welcher 
sich um eine Expansion von exakt e = 5° von dem Anströmungszustand unterscheidet. Um 


a a end nun nn 


Bild 11. Strömungsebene der Expansionsströmung im Normalfall, b/a = 0,750 


festzustellen, von welchem Punkt eine Prandtl-Meyersche-Expansion ausgehen kann, 
nehmen wir zuerst an, daß ein Expansionsfächer von der Vorderkante der Laufschaufel ausgeht. 
Dieser Expansionsfächer wird zunächst an der Unterseite der Leitschaufel reflektiert, und in- 
folgedessen entsteht hier ein kleinerer Druck als auf der Oberseite der Leitschaufel. Das Druck- 
gleichgewicht über der freien Stromlinie hinter der Leitschaufel erfordert, daß sich ein rechts- 
laufender Stoß und ein linksläufiger Expansionsfächer an der Hinterkante der Leitschaufel aus- 
bilden. Damit ist indessen der Anströmzustand vor der Vorderkante der Laufschaufel schon 
geändert! Um nun auch die Randbedingung auf der Laufschaufel zu erfüllen, muß der linksläufige 
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Expansionsfächer an der Hinterkante der Leitschaufel exakt e° sein. Damit besteht aber kein 
Grund mehr für den zuerst angenommenen Expansionsfächer. 


Die Konstruktion der Strömung für 2 — 0,75 ist in Bild 11 dargestellt. Die Verwandt- 


i hier behandelten exakten Normalfall der axialen Rückwirkung und dem 
ed offenbar. Statt einer einzigen Machlinie bei der linearen Betrachtung 
treten hier Expansionsfächer auf. Berechnet man z.B. den Zustand im Feld 2, Bild 4, für 
e = 5°, & = 44°11’, so bekommt man 

yes 1+tgatge — 1,186 
Wosy Uytarı IE te 
im Vergleich zu dem exakten Wert: 


W 
— = 1,132. 
W 


0 

Für kleine e-Werte stimmen die Zustände im Normalfall mit einander gut überein. 

Dieser Lösungstypus in Bild 11 kann jedoch nur innerhalb gewisser Grenzen von bja 
existieren. Falls b/a einen solchen Wert hat, daß die erste Machlinie in dem an der Laufschaufel 
reflektierten Expansionsfächer vor die Laufschaufelvorderkante gelangt, beginnt sich das 
Strömungsbild zu ändern (Bild 12). Damit entsteht ein komplizierteres Problem. Um das 
Strömungsbild in Bild 12 zu erklären, wählen wir zunächst ein solches Verhältnis b/a, daß das 
Ende der Leitschaufel in dem an der Unterseite der Leitschaufel und an der Oberseite der Lauf- 
schaufel reflektierten Expansionsfächer liegt. Damit ändert sich der Zustand vor dem Stoß, und 
zwar so, daß wir unmittelbar vor dem Ausgangspunkt des Stoßes einen noch kleineren Druck als 
vorher bekommen. Um das Druckgleichgewicht über der freien Stromlinie hinter der Leit- 
schaufel zu erhalten, muß die Expansion an der Oberseite etwas größer als vorher werden. 

In Bild 12 ist nun die Prandtl-Meyersche-Expansion im Hodographen auf 7,5° fest- 
gelegt. Die zu diesem Expansionswinkel gehörige Stellung b/a wird danach in der Strömung in 
folgender Weise aufgesucht. Mit einer Expansion von 7,5° liegt das ganze Hodographenbild 
völlig fest, und das Aufsuchen der zugehörigen Stellung besteht nunmehr darin, die Charak- 
teristiken zu konstruieren. Dabei sind allerdings die Richtungen der Charakteristiken in der von 


III 


Bild 12. Strömungsebene der Expansionsströmung, b/a = 0,470 


Stößen durchkreuzten kontinuierlichen Strömungsebene nicht a priori gegeben. Man kann 
also nicht ohne weiteres die Stellung b/a angeben, sondern man muß einen iterativen Weg ver- 
wenden, um die Stellung aufzusuchen. Dabei spielt natürlich die Charakteristik, welche an dem 
Ausgangspunkt des Stoßes an der Hinterkante der Leitschaufel auftrifft, die größte Rolle. Auf 
diesem Weg läßt sich die Stellung b/a = 0,47 aufsuchen. Vergleicht man die exakte Lösung 
(Bild 12) mit der linearen Lösung (Bild 10b), so findet man im wesentlichen dieselbe Konfigu- 
ration des Strömungsfeldes; im linearen Fall liegen die Stöße parallel, im exakten aber nicht. 

Mit einem noch kleineren Wert b/a findet man, daß der Expansionswinkel zunimmt. In 
Bild 13 ist eine solche Stellung b/a aufgesucht, daß die Vorderkante der Laufschaufel gerade mit 
€ = 9. angeströmt ist. Dabei ist der Expansionswinkel 13,2°. Dieser Winkel ist nun so groß, 
daß der linksläufige Expansionsfächer das ganze Rückwirkungsgebiet überdeckt; denn die letzte 
Machlinie in der Prandtl-Meyerschen-Expansion endet an der Unterseite der Laufschaufel. 

Die Erklärung für diesen großen Winkel liegt darin, daß jede Machlinie nach Reflexion an 
Leit- und Laufschaufel das Rückwirkungsgebiet verläßt. Falls diese Bedingung nicht erfüllt 
wäre, würden ähnliche Verhältnisse auftreten, welchen wir bei der linearen Betrachtung begegnet 
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sind. Dort hat sich der aus der Grenz-Machlinie entstehende Machlinienzug durch mehrere 
Teilungen im Rückwirkungsgebiet erhalten. 

‚Physikalisch ist dieses Strömungsbild auch einfach zu erklären, denn bei einem kleinen 
Verhältnis b/a nimmt die Strömung den einfachsten Weg. Dieser Weg ist tatsächlich der ‚untere 


Weg“ direkt zwischen zwei Laufschaufeln. Um diesen Weg nehmen zu können, muß der Expan- 
sionswinkel groß sein. 
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Bild 13. Strömungsebene der Expansionsströmung, b/a = 0,218 


In Bild 14 ist das Verhältnis b/a so gewählt, daß der in Bild 13 auftretende Stoß vor die 
Vorderkante der Laufschaufel gelangt. Bei dieser Stellung ist der Expansionswinkel etwa 
maximal und gleich 16°. Der entsprechende lineare Fall ist in Bild 10e dargestellt. Der Unter- 
schied zwischen diesen beiden Konfigurationen ist wieder der, daß bei der exakten Betrachtung 
die Stöße zusammenlaufen. 

Mit einem noch kleineren Verhältnis b/a erhält man einen Strömungstypus, welcher der 
reinen Kanalströmung ähnelt. Diese Lösung ist mit jener der linearen Lösung qualitativ über- 
einstimmend, Bild 10d. 


Bild 14. Strömungsebene der Expansionsströmung, b/a = 0,107. 


Überhaupt haben wir bei einem Vergleich der exakten Theorie mit der linearen Theorie 
folgendes gefunden: Bei solchen Stellungen, bei denen wir mit der linearen Theorie eine gültige 
Lösung gefunden haben, sind die Strömungsbilder der linearen und der exakten Theorie einander 
qualitativ gleich. Bei solchen Stellungen, wo die lineare Theorie versagt, findet bei der exakten 
Behandlung eine „Überexpansion‘‘ des strömenden Mediums statt, welche vielfache Reflexionen 
über mehrere Teilungen von Machlinien im Rückwirkungsgebiet verhindert. Damit können wir 
die Untersuchung der Expansionsströmung abschließen. Wir werden uns zunächst mit dem 
Kompressionsproblem befassen. 

8. Das Kompressionsproblem 

Um eine Kompressionsströmung zu bekommen, müssen die Laufschaufeln einen negativen 
Anstellwinkel haben. Wir legen diesen Winkel auf & = — 5° fest und studieren nun die Strömung 
bei einer Machzahl, welche wir noch nicht bestimmen wollen. Es wird sich nämlich zeigen, daß 
die Machzahl innerhalb gewisser Grenzen von b/a hoch genug sein muß, um eine Lösung im 
Überschallgebiet zu bekommen. Der ß-Winkel wird auf $ = 30° festgelegt. 

Dieselbe Überlegung wie bei der Expansionsströmung ergibt nun sofort, daß für b/a nahe 
an eins die Strömung durch linksläufige Stöße von den Hinterkanten der Leitschaufeln um den 


Strömung im Überschallgebiet : 
lich, mit der Charakteristikentafel die Grenz- 


a 


Au Bild 15. Strömungsebene der Kompressionsströmung im Normalfall, b/a = 0,6 i m 


Jeder Stoß entspricht dann einer Verminderung der Ma chzahl um e° inder Charakteristikentafel. a 
Dabei gibt es nun drei Stöße; man kann den Mindestwert der Machzahl für eine Überschall- Be. 
Lösung unmittelbar aus der Charakteristikentafel ablesen. In Bild 15 ist die Konstruktion der 


De % 
Strömung für M = 1,775, und ir 0,6 durchgeführt. 
Das so erhaltene Strömungsbild ändert sich indessen stark, wenn b/a nur einen etwas 
kleineren Wert hat. Schon bei u 0,525 läßt sich die Überschallströmung nicht länger fort- 


setzen, denn dann tritt der an der Oberseite der Laufschaufel reflektierte Stoß vor die Vorder- 
kante der Laufschaufel. Die Stöße laufen dann zusammen und das ganze Strömungsbild ist 
verändert. Man erhält eine Unterschallströmung und die Charakteristikenmethode versagt. 


b : a a 
Die oben gefundene Lösung, gültig für 1 > z > 0,525, ist der einzige Lösungstypus, 


welchen man im Gitter bei dieserMachzahl finden kann. Falls man eine noch höhereMachzahl 
nimmt, zeigt sich, daß man für ein sehr kleines Verhältnis b/a noch eine Lösung finden kann. Der 
Charakter dieser Lösung ist mit der reinen Kanalströmung übereinstimmend. 


Für eine genügend hohe Machzahl findet man also zwei Lösungstypen innerhalb einer 
Teilungsperiode. Weiter gibt es für b/a Werte, gelegen zwischen den Gültigkeitsbereichen dieser 
beiden Lösungen, keine Überschall-Lösung des Strömungsproblemes. Im Vergleich zur linearen 
Theorie, welche für & negativ etwa dieselben Resultate wie für & positiv liefert, bekommt man j 
folgendes: Für solche Stellungen b/a, wo die lineare Theorie eine gültige Lösung geliefert hat, | 
sind die Strömungsbilder einander gleich. Bei Stellungen, wo die. lineare Theorie versagt, geht 
die Strömung bei der exakten Behandlung in den Unterschallbereich über. Vorausgesetzt ist 
jedoch, um eine Lösung überhaupt im Überschallbereich zu bekommen, daß die Anström-Mach- 
zahl hoch genug ist. 


9. Zusammenfassung 


En 


Es wird zuerst die Strömung vor und durch ein Schaufelgitter ohne Leitrad nach der zwei- 
dimensionalen Charakteristiken-Methode in dem mit den Schaufeln fest verbundenen Koordinaten- . 
system diskutiert. 


Wenn dieKomponente derMachzahl senkrecht auf den Schaufelkranz kleiner als eins ist, 
laufen Machlinien stromaufwärts in das Anströmgebiet. Dort können dann nicht mehr Strömungs- 
winkel und Machzahl vorgegeben werden, sondern nur eine Kombination dieser Größen. Weil 
aber das Gas drehungsfrei einströmt, bedeutet das, daß die Einström-Machzahl von der Umlauf- 
geschwindigkeit des Laufrades abhängt. 


Die stationäre Strömung durch ein Gitter mit Lavaldüsen oder Leitschaufeln wird bei 
verschiedenen Stellungen vom Leitrad zum Laufrad ausgedehnt behandelt. Die stationäre 
linearisierte Überschalltheorie liefert bei gewissen Stellungen keine brauchbaren Lösungen. Dabei 


ist es wesentlich, ob die Randbedingungen exakt oder nur durch Quellbelegung von Skelettlinien 
erfüllt werden. 
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R Der exakte Expansionsvorgang am Ende der Leitschaufeln oder Düsen ist dargestellt für 
jede Stellung vom Leitrad zum Laufrad. Es wird auch gezeigt, daß eine Kompressionsströmung 
im reinen Überschallbereich nicht für jede Stellung vom Laufrad zum Leitrad existiert. Schließ- 
lich wird die lineare Lösung mit jenen der exakten verglichen. 

“Der Grund, weshalb die Erfüllung der exakten Randbedingungen sowie die Anwendung 
des exakten nichtlinearen Charakteristikenverfahrens selbst bei sehr kleinen Störungen in 
Gittern wesentlich sind, ist folgender. Die lineare Lösung wird in großem Abstand von der 
Störung falsch. Eine Änderung in der Formulierung der Randbedingung oder eine kleine Ände- 
rung der Machlinienrichtung führt dazu, daß keine geschlossenen Machlinienzüge im Rück- 
wirkungsgebiet entstehen können oderMachlinienreflexionen über mehrere Rückwirkungsgebiete 
hinweg verhindert werden. Das bedeutet aber, daß Aufschaukelungen in der exakten Theorie 
des Druckes wegfallen. Nur die Kompressionsvorgänge sind in vielen Stellungen äußerst empfind- 
lich gegen Störungen. Solche führen schnell auf Unterschallströmungen und entziehen sich 
einer Behandlung mit Charakteristikenmethoden. 


10. Schlußwort 


Die hier wiedergegebenen Lösungen für ein Gitter mit Leitrad sind nur im quasistationären 
Fall gültig. Falls das Laufrad eine Relativgeschwindigkeit im Verhältnis zum Leitrad hat, 
bekommt man ein instationäres Problem mit laufenden schiefen Machfronten. In einem solchen 
Fall treten im Gitter zweidimensionale instationäre Wellen auf, welche die oben gegebene Lösung. 
im Gitter beeinflussen. In einer kommenden Arbeit werden wir uns mit dem Einfluß der Relativ- 
geschwindigkeit auf die Strömung in einem Überschallgitter befassen. 

Diese Arbeit wurde auf Anregung und unter der Leitung von Herrn Professor Dr. Klaus 
Oswatitsch ausgeführt. Einen Teil der Ergebnisse in Kapitel 2 und 4 hat er früher gefunden 
und auf der V.D. I.-Fachsitzung in Karlsruhe 1953 gebracht. Für all seine Hilfe bei der Verall- 
gemeinerung der älteren Resultate und bei der Durchführung der Arbeit im ganzen möchte ich. 
ihm hiermit danken. 
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Nonuniform Rotation of an Infinite Circular Cylinder 
in an Infinite Viscous Liquid 


By D. D. Mallick in Calcutta 
(Ghosh Research Scholar, Department of Applied Mathematics, Caleutta University) 


Es werden strenge Lösungen gefunden für die unbeständige Anfangsbewegung, die durch die Rotation 
eines Kreiszylinders in einer zähen Flüssigkeit erzeugt wird. Dabei wird sowohl der Fall behandelt, daß die 
Flüssigkeit den rotierenden Zylinder umgibt, wie auch der Fall, daß sich die Flüssigkeit innerhalb des rotieren- 
den Zylinders befindet. In beiden Fällen wird ein stationärer Zustand erreicht. Wenn die Flüssigkeit den 
rotierenden Zylinder umgibt, breitet sich der stationäre Zustand vom Unendlichen her in Richtung auf den 
Zylinder hin aus, während im anderen Falle der stationäre Zustand von der inneren Zylinderfläche aus gegen 
die Rotationsachse hin fortschreitet. Für beide Fälle wird der Ausbreitungsvorgang des stationären Zustan- 
des durch Kurven dargestellt. 


Rigorous solutions have been obtained giving the initial unsteady motion produced by the rotation:of a 
circular cylinder in a viscous liquid. Both the external and internal solutions when the liquid is outside or 
inside the rolating cylinder have been obtained. In every case a steady state is reuched. Ourves showing how 
this steadiness spreads from infinity towards the cylinder for outside motion, and from the cylindrical surface 
towards the axis for internal motion, have been drawn. 


On trouve des solutions exactes pour le mouvement initial et inconstant produit par la rotation d’un cylindre 
eirculaire dans un liquide visqueux. Ieci on traite le cas d’un liquide entourant le cylindre en rotation aussi bien 
que le cas du liquide se trouvant dans le cylindre en rotation. Dans ces deux cas on arrive 4 un Etat stationnaüre. 
Si le liquide se trouve autour du cylindre en rotation V’etat stationnaire s’Etend de Vinfini dans la direction du 
cylindre, tandis que dans lautre cas l’&tat stationnaire procede de l’etendue imterieure du cylindre vers l’axe de 
rotation. Pour les deux cas le procede& de l’extension du cas stationnaire est illustre par des courbes. 


OupeneAAloTcH CTPOTHE penieHus IA HeyCTOÜYUBOTO HAYAJIbHOTO ABMFKEHNSI, BhI3bBIBAeMOTO 
BpallieHmeM KPyTOBOTO UHIHHNPaA B BASKOH MinnRkocTu. Ilpm 3TOoM paccMmaTpuBawTcH Kar 
cıyyali, Korna Bpamarommmüca IMIMHAP HOTPY>KeH B ZRUNKOCTB, KaK N c/lyuali, KOTMA FRUN- 
KOCTb HAXONHUTCA B Bpalaroımemca ımmHnpe. B 060mx cITyyaax NOCTUTAeTcA CTAmMOHApHOe 
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COCTOAHHE. Ilpr UHJIMHAPe, HOTPY>ReHHOM B HKHIKOCTB, cTauuHOHApHOE COCTOAHHE Pac- 
IPOCTpaHAeTcaH, TIPNHÖJIMRAACh KR IAJIHHAPY U3 6eCKOHKYHOCTH, A BO BTOPOM cIIy4ae CTALUMO- 
HapHoe COCTOAHHE pacıpocTpaHufneTcH OT H3HaAHRU UHIUHNPHYeCcKoü HOBEPXHOCTUH RK OCH 
BpalleHuf. B 060nux ciyyaax ABJIeEHMeE PAacnpocTpaHeHuA cTa lMOHAPHOTO COCTOAHHA 


U306parKaeTchH RKPHUBbIMN. 
1. Introduetion 

The problem of steady motion of an infinite viscous liquid due to the rotation of an infinite 
circular cylinder about its axis has been studied by Couette. In the present paper we propose to 
study the unsteady motion when the cylinder is set rotating by an initial impulse. Two cases 
have been studied. Firstly (1) in an infinite viscous liquid rotating with constant angular velo- 
city Q about an axis, say the axis of z, an infinite circular cylinder with its axis coincident with 
the z-axis is supposed to be given an impulsive twist at time ? — 0 such that it rotates with uni- 
form angular velocity 2,, about the axis; secondly (2) an infinite circular cylinder containing a 
viscous liquid is suddenly made to rotate with constant angular velocity 2, about its axis. In 
both cases rigorous solutions have been obtained which show how the motions become steady. 
Theoretically the steadiness comes in infinite time and the steady solutions agree with the known 
Couette motions. 

The rigorous solutions contain integrals which it has been possible to evaluate only numeri- 
cally at definite distances from the centre and at definite times. In the case of the first problem 
of the cylinder set rotating by an impulse in an infinite rotating liquid, it is shown that at any 
subsequent instant the points in the liquid further from the axis are closer to steadiness than 
points which are nearer as one may expect; the solution really shows how the steadiness gradually 
spreads from infinity towards the rotating cylindrical surface. A similar study of the second case 
shows how the ultimate state of the liquid inside the cylinder rotating as a rigid body with the 
liquid is attained. In this case the behaviour of rigidity spreads from the rotating surface towards 


the axis. 
Case I 


2. Formulation of the problem 
Let us introduce cylindrical co-ordinates r, ©, z with z-axis along the axis of the cylinder, 
and let u, v, w be the components of the absolute liquid velocity in these respective directions. 
Since the motion is two-dimensional we can put w = 0. Further the motion is independent of ©. 
The equation of continuity for incompressible liquid is 


(0) 
ar Fü 
hence ru = constant = 0, so that u = 0. The equations of motion become 
v: 1 op 
OO . Ze a v 
RT, ara Ninanmlu > ee 


where o is the density of the liquid, p the pressure, and » the kinematic viscoeity. 
ES Before the impulsive twist is given to the cylinder immersed in the liquid, the liquid outside 
it is in a state of steady motion, so that equation (2) becomes 

De de v 


= eerif) 
dr? r dr r® 


of which the solution is 
A 
I a —- B RR 


From this we get for. the boundaty conditions v = 0,onr = a,andv = Q r,asr— 00, the well 
known solution 


a? 
= Oip- — 
v 2lr 2 ee 
giving the velocity at any point of the liquid. 
The boundary conditions of our present problem we take as 
on ratg, Da for alt >0 | 


as r— oo, v— Qr for any fixe 
ans y fıxedi 


| ..(4). 
at Bl), = 2) ior.r, a 
? 
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3. Solution of the problem 
We shall make use of Laplace transform to solve the problem and write 
vr, s) = e*u(r, f) di 
ö 
By this transformation equation (2) becomes on simplification 
rn a’v T dv 
dr? dr 


5 )»-—2 > 
(r&+1)5= „e GET EE (5) 


the boundary conditions (4) now become 


The solution of (5) is of the form 


Be tan) )r AK; (VE) Be  V. (7), 


where I, and K, are the modified Bessel functions of the first and second kind. The boundary 
condition (6) give 


an, 


Hence 


Ss 
a? «0,2 ) 


or # i£ Are. N: 
Ss T s 
SKK, ( v: 
v 
The solution of our problem is obtained by inverting the transformation; we then have 
cHo | 5 [ 1 
Q Sc] #1 e73 
D.= an ee [ = l — a2 ds 
2ri R/ıs s 
/ S 7 ( VE 
x en 9); 


J 
where c is a positive constant. 


K, has no zeros in the s-plane, and s= 0 is a branch point of the integrand in the finite 
part of the plane. We shall take the contour of integration to be an infinite semicirele with 
R(s) < 0, together with a cut from s = 0 to — © as shown in fig. 1. ER 


en Bar | 
| ie [ 1 | 
a 2 | st 2 wert. A d N 
of kil# ++] +fe als 
anti, sK,la y: | 
| [OR C, C, c o f v = | 


| e—=iR 


—= 2ri N residue at the poles =. 


It is clear that 


[0 as Ro 
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e er w 

s SU 1m als Sen N 1 

x (eV) «ı. SZ nasM jelg-Vs)+r-&m Sn+ıj 2 
a. i 2n+2 

Kl u ; (2a) %) 1 1 1 

7 ea boe[LayE Hy Shan +3 


| 1 
Ss 
Klr 2) 40. +5] 
20 E is 21 te) 


Irni 5 2nri | ü 
& KlaV) E s|A(a,s) + p) 


where 


r 
—r — — n 
AU 1 > 1 1 
A) => n!(n + 1) hos(4 v ) = &m I(n+1) 
Now, as the radius of the circle C, tends to zero 

est 


a2 Aut, S);e& ee a > A a? 2 
2rri 


A Zr sla@9 +5) 


On 2 s=ae’” andon ev si au) Therefore 


7 = Et Fe 
si a d ni [ 2 | | 2 Se 7%) er, 
s= - = nn 
— ni i & al & 
2 a r Ne (i «y%) RK, —iay® 


where we have used the relation 


ee 
Ka a UL 


Hence from (9) we have 


4 200,00 [lFnll3)-aeVE)neE)- 


ee SEN SETER 2: so: dx 


dh In(ay&) + n{aV2)) 


N De 9, 2a9 [ {lra) near) — lar) la) , 
r: £ 7 [ x{J:(ax) + y(ax)} 

a? a0 2 

N) 5 1=-0Qr - (@ — Q2,), showing that a 

steady state is reached on which Couette motion is established. 


When ?— 0, we have from (10) 


Den [ & \ ed al , 2a, BAG X) Yılaz) — Jılax) Yyılr 2) ir 
r ra x {J:(ax) + y%(ax)} 


—ayt dx (10). 


Ö 


It>is elear that sr am o(r- 


(N) 


er ut ee ee 
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| - 
But from the boundary condition (4) v would have to be equal to 2 ( 2). Comparing these 
r 
two results we note an interesting evaluation of the following integral 
[ Ir 2) yıaz) — (az) nlra) wa 
re ee R 
| 2{J(a2) + yraz)} B ar: Nor fra ae 


0 


In course of subsequent numefical integration this has been verified for - 2, 3,4. 
a 


The graph of the integrand of the equation (10) for vt = 0.1,a = Li,e of y e-0%12*, where 
_ AL) Yo) — A) yılr ©) 
TR) + HR} 
has been drawn in fig. 1.1 forr=2, 3, 4. The areas in both cases are negative and itis found from 
the fig. 1.1 that the total area (taking the sign into consideration) corresponding tor—=2is 
greater than the areas correspondingto r=3,andr =4. Y 
We may write equation (10) in the form 


3 aa an, 
= @lr ee r 
=-I+A+T 


where I is the initial velocity (before the cylinder is 

set rotating), A the part added to / when steadines o 
is attained, and 7 the transient part which ultima- 
tely disappears. T is negative, steadily approaches 
zero with time, and for fixed t, decreases as r increases. 
From the measurement of the areas of the curves in 
fig. 1.1 the ratio T/A can be calculated for-r—=2,3 
and 4. It is found that the numerical value of this 
ratio is less for greater r. At a given time points which -0% 
are at greater distance from the cylinder are nearer to 
steadiness than the points at lesser distance. One may 
say steadiness spreads from infinite distance towards 
the cylinder. In considering the variation of v with 
time one notes that T7 being negative, the steady 
velocity at a given point is estabilished by a con- 
tinuous increase of velocity till the total increment A -03 
is reached. Fig. 1,1 


FT 


02 | 


-06 


-08 


Fig. 2 


The graphs of the function y have been drawn in fig. 2, for r = 2,3. By Newtons inter- 
polation formula the function y for r = 2, has been represented by a polynomial of degree six in 
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the form 4 
y = 0.0029 x5 — 0.0365 x® + 0.1984 x — 0.634523 + 1.3128. 2? — 1.3448 X. 


Then evaluating the integral in (10) frr= 23, a = 1, we obtain for the velocity the expansion 


= 


30 a ae ss 
or zu a 


1 3 nr 
ee I" _ x.0.1984 ErENaTH 
En * (11) 


SE 
ee = 


which should be valid for large vf. Equation (11) shows how the transient part of the velocity dies 
out with increasing »vt. 


0 


& 
7 


=02 


Fig. 3 


In fig. 3, the graphs of the integrand forr=2, a=1 and forv»t=0, and 10 have been 
drawn. The ratio of the area for v{= 10 to that for»t=0 is about = . This shows that of the 


total increase of velocity A required to reach the steady state at distance equal to twice the 
radius of the cylinder, 94% is attained by the time » ? has the value 10. 


4. Determination of the frietional couple 
The frictional couple per unit length of the cylinder r = a is given by 


d /v 
Now Pro=urz (2), where v» = uJo; Hence 


. 


[Prol,_ = 2u — "yn(a2) Iolax) — Ylaz) Jıla2)} 
r=a 7 


Sax) + ya) 


ee’ vt 2 
0 
giving 

oo 


I 2, Ara? [ (yılaz) Islax) — yılaz) Jılax)} 
(= An u 2 = ®) 2 | ı Se Ne 0 ! —ı?V 
u \ A 7 J Jı(ax) + yilax) Fe ie ne: 


0 


The graph of the integrand in the expression for C which gives the transi 
RR : sıent part oft 
been drawn in fig. 4 for vi = 0.1,0.2,03,anda —1. 8 part of the couple has 
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Fig. 4 


: 5. Case II 
£ We assume a viscous liquid inside a circular cylinder is set in motion at imet=(,bya 
rotation Q, imparted to the cylinder. The motion as in the previous case must be independent 


= of © and as before we can take u = w= (0. The equations of motion can now be written as 
v® 1 op 
a Te ee LE EN, er RE 
le al Tabea en NEL: 


a Fr a 


The boundary conditions for this problem are 


Dar; neo, torralh deal) 
and . (16). 
ER Da for De 
Making use of Laplace transform as in Gase I we have from (15) and (16) 
De Te a TR 
Eee ee en RE rt 
v dr? u 9 en 
and 
( 
on -r=a, EINER IRRE BSR REEL: 
z 


The solution of the differential equation (17) subject to the boundary condition (18) is 


s 
2,a A| V-) 


v— VE (I), 
s (af?) 
v 


where /, is the modified Bessel function of the first kind. 
v is determined from (19) by inverting the transformation, and we have 


c+io u 
Da I; ( v: 


en = (a ÜBEN ehr rk Enns Be (ZU 
ri R 
c—io S I, a > 


I, (‚VE is not singlevalued but the integrand is. The poles ofthe integrand are given by 
v 


VER, 
Se 10» een les ER? = 
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where j„ s are the zeros of the Bessel function of the first kind and first order, J,. The residue at 
N 


the pole 
sten); is — 


vi on 
eg Tr . 
exe a 
v a 


Er 
£ ee DL 


vL er) 
Dez Tr. 
oo e r len) 


>23 1 N JG h) 


ie: 
© (2) 


and that at the pole 


Hence from (20) we have 


v=aQ|—+2 ea 


We verify that at {= 0 


r 
=( AP) a 2 re ee | = 0 s 
rl A 
since by Fourier Bessel expansion we have 
Wo 
r = A, 
= = (>), 


2a er inJolin) 
The steady state is reached with {— oo, when v— Q, ri. e. the liquid inside the cylinder rotates 
as a solid body. The following Table shows how for different times given by the different values 
ofv»t= 0.1, 0.2, 0.3, and 0.4, the unsteady part of the motion varies at different distances from 


the centre. These distances are chosen as — — 0.4, 0.8, and 0.99, the last value fixing a point 


at ninetynine percent of the radial distance from the centre i. e. almost on the’ cylinder. The 
transient part of v represented by the summation term has been calculated for each value of 
vt for the above values of r/a. It is found that the first two terms corresponding to J, and 7, give 
a fairly approximate value of the summation in (21). T(0.1), T(0.2), T (0.3), and 7(0.4) give 
the values of the transient terms for vt = 0.1, 0.2, 0.3 and 0.4, whereas S(0.1), S(0.2) etc. give the 
ratios of the transient to the steady parts of the velocity 


Sie 
a1; ea 


VIE 
Yn | vi=(.l vt=(0.2 vi=0.3 vt—=04 
r = = 
yı | %» TA) | s(0.1) | 7(0.2) | s(0.2) | (0.3) | s(0.3) | 7(0.4) | S(0.4) 
0.4 | 0.3643 | 0.1936 | 0.0825 | —0.413 | 0.0193 | 0.007 | _0.0045 | _0.023 | _0.0010 | _0.005 
0.8 | —0.2050 | —0.1583 | —0.0484 | —0.121 | 0.0109 | —0.027 | 0.0025 | —0.006 | —0.0006 | —0.002 
0.99 | —0.0102 | —0.0101 | —0.0024 | —0.005 | —0.0005 | 0.001 | 0.0001 | —0.000 | 0.0000 | 0.000 


The Table shows that at any instant of time the absolute value of the transient term T dimini- 
shes along the radius from the centre towards the cylinder and so also does the ratio of this term 


to steady value at the point. Thus one may say the steadiness spreads from the cylinder towards 
the centre. 


One can find in the usual way the frictional couple per unit length of the cylinder. This 
calculates out to 


C=Anunaye #" SE 


In conclusion the author wishes to thank Prof. N. R. Sen for his continued interest and the 
benefit of discussion with him in course of the work. 
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ae Tee Finite Bendi 


RE RER. . 
3 of a Plate into a SphericalShell_—— 
ng ofa Plate into a Spherical Shell 
By B.R. Seth, Indian Inst. of Technology, Kharagpur 


Mittels der Theorie der endlichen Deformationen wird gezeigt, daß eine Kreisplatte vom Radius d in eine 


Kugelschale verbogen werden kann, und zwar allein durch solche Kräfte und Kräftepaare, die lediglich an den 
Rändern der Platte angreifen, während die innere Oberfläche a und die äußere Oberfläche b der Kugelschale 
kräftefrei bleiben. Der Radius der Platte wird als klein im Verhältnis zu a und b angenommen, so daß (d/a2) 
. vernachlässigt werden kann. Eingehend diskutiert wird der Spezialfall, in dem die Poissonsche Zahl o gleich 1/3 
ist. Die prozentuelle Abnahme der Schichtdicke sowie die Größenänderung der Kantenradien wurden errechnet. 


By using the theorie of finite deformation it is shown that a circular plate of radius d can be bent into a 
spherical shell by forces and couples applied to ihe edges only, the inner surface a and the outer surface b of the shell 
being free from tractions. The radius of the plate is assumed to be small compared with a and b so that (d/a?) 
cam be neglected. The particular case in which the Possion’s ratio o is equal to 1/3 has been discussed im detail. 
The percentage of decrease in thickness and the change in. the edge-radii have been obtained. 


Au moyen de la theorie de la deformation finie il est demontre qu'une plaque circulaire du rayon d peut 
Etre deformee en un coupe spherique, et ceci uniquement par des forces et par des couples de forces attaqguant 
purement aux bords de la plaque, tandis que la surface interieure a et la surface exterieure b du coupe spherique 
restent sans influence de forces. On suppose que le rayon de la plaque soit petit compare avec a et b, de sorte 
qu’on puisse negliger (d/a?). On discute ü fond le cas oü le chiffre de Poisson o est egal & 1/3. La calculation 
de la diminution per, cent de l’Epaisseur de couche ainsi que calle de la variation de l’Etendue des rayons des 
bords ont ete faites. 


Ilpı moMolım TeopHuM KOHEeYHEIX MedopMauma NHOKA3bIBAeTCcH, UTO KPyToBaf INTACTUHRAa 
Panmyıa d MOKeT ObITb U30OTHYTa B YacTb ChepnyecKof 000JI0YKU, IPHUY&M CHIIBI U ITaPbI aHTU- 
mappasııelIbHbIX CHJI IPHKJIANBIBAIWTCH TOJIBKO K KpAlO IJIACTUHKH, B TO BPeMA KaK U U3HAHRAa a 
N JmmMeBas CTOpoHa b O000N04KM CHJIAM He HonBepraiotcn. Pannuyc NNACTUHKU IPUHNMaertc 
MAJIbIM IIO CPaBHEeHNHIO C Ad U 6, TAK YTO Berimymmoi (Ö/a2) npeHe6öpeyb. |leTasıbHo paccMmaTpnH- 
BaeTcaA cAyyal, Korla HoMTOAHHan Ilyaccona 6 pasHa 1/3. BeIuncaamoTtca Take IIPoleH- 
TyasIbHbIe YMEHbIlIeHNME TOJIIHHBI CIHOA U U3MEHEHHA BEJIHYMHBI PalmycoB P&öep. 


1. Introduetion 


In recent years the theory of finite deformation has received fresh impetus from a number 
of exact solutions that have been obtained for incompressible elastic bodies without assuming 
any form of the strain energy function. The exact solutions obtained for compressible bodies 
have not been very large, and in general they have been of the form obtained by B.R. Seth. 
Extensive references to existing literature on the subject are given in the Presedential Address to 
the Section of Mathematics at the 42nd session of the Indian Science Congress [1]. 

Very few attempts have been made to obtain solutions for the finite bending of plates 
into shells. The earliest exact solution on the basis of a linear stress strain tensor law was obtained 
by Seth [2] for a rectangular plate bent into a right cylindrical shell. The form of this solution 
has been used by many workers. In the pressent paper an attempt has been made to discuss 
solutions which may be used to discuss the finite bending of plates into spherical shells. A linear 
stress strain tensor law is again assumed, and this gives sufficiently good results for technical 
applications. It can be shown that on this basis all the stress equations of equilibrium and the 
boundary conditions can be satisfied by a set of displacements in cylindrical coordinates r, Ö, z 
given by 

i 

Ten ur, v=9—AD, w=z— B(r +2)? Mehr 

where a, A, B, D are known constants and c is an elastic constant. In this solution planes parallel 
to the faces of the plate go into spherical surfaces, but straight lines parallel to the axis of the 
plate remain parallel to the axis of the spherical shell. The solution can therefore apply, only to a 
plane membrane deformed into a spherical one. For a plate lines parallel to the axis of the plate 
should go into portions of the radii of the spherical shell. We show that such a solution can be 
obtained for a plate, if squares and higher powers of d/R? can be neglected, d being the radius of 
the plate and Rtheradius of the spherical shell. The resultant traction and couple to be applied 
to the edge are worked out and some numerical examples are given. 


2. Strain Components in Cylindrieal Co-ordinates 
If wi is the displacement vector the finite strain tensor is given by [2] 
2; tut u? — 29 uf et (2.1), 


where 9 and g;,; are the metric tensors corresponding to the unstrained and strained state. In cylin- 
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drical coordinates we have k 
y—=r,d, 2; 0-7 ®,z 
u=r—r, v=9-—-%, v=2—2, 
Ir ae 
;=|0 nn 0), m; 020 22) 
DE Da 


The tensor components can bereplaced by the physical components with the help of the relation[3] 


1 
6; = (gig)? &,,: 


Thus we get strain components in cylindrical coordinates. For example, &,, is given by 


; ou oe oo u a” ,„/ov\” 93 
2 u Ka ee): 


For the problem in hand we should take 
u=r— Or,d), W=d— AP, , w=z—/(R + Dur. ae 24, 


where ® and / are arbitrary functions and A is an undetermined constant. 
The corresponding strain components are 


1 ; 12492 
„sl M-2rf, w-hghr | 
1 Bar 
a er tar ee, =—-—- 909, —2rzf, Ge 
zZ > 2 
eo; = (0) 9 e,o —ı 0 
where 
Ki a (1—4%), K=1—-2(—-4A)—- (1 — 4), 
3(2.6% 
od od el are 
a Ta a a 
3. Stress-Equations of Equilibrium 
We take a linear stress-strain tensor law given by 
T,=16,4A+2ue,, RN 


where 
A=e,4+6%s+t %;- 
This is found to be quite satisfactory for technical applications and has also been used by 


Synge and Chien. 
The stress equations of equilibrium 


ee Be) 
give the two differential equations 
op, u 5 = A LE, z 
=> E- a ( —2k,-+kı = u®; A | : HN) (3.3). 
op na.fle Di 
2 Fur®; ar I 2) —0 (3.4), 
where 
p=AA—8uf fr ARTE uf Re Par (Dee (3.9) 
p = constant is a solution of (3.4) if 
BD, Rs 
Gere) (3.6), 


x being an arbitrary function of z. 
Equation (3.6) combined with (3.3) and (3.4) gives the only exact solution 


1 
a vi dr Ar w=z— B(R +22" ”_E ee 100 


nd a ah Are 
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where ‚the constants a, A and B have to be determined from the boundary conditions. This 
solution is not suitable for our problem. We require a solution in which ©(r, z) is a function of 
(r/z). This condition requires 


| ER Ka ER, 
To satisfy (3.3) we should have 
| | 1 k, ®: dB; 2° 
Ant - Deep (er) UN BR EN (3.9). 
Combined with (3.6) and (3.8) this gives h=5, 
and 
[ Tr Vi 
9=2D | ED. 2 ee er 3. (8-10), 
z+ (@ +)? 
D being a constant of integration. 
1 
ko on gives from (2.6) A = 1, k, = 1. Hence 
a 
2 2 Pa IERER 
Dean -=2D: an el, tan Sy LT), 
z+ (2 +r)2 
where tan y = (rJ2). 
The equation p = constant now gives 
- h PER DEIN 
li +kh—2(r +22) 5 k, ar hun (DB; 2] 
— 8uff?dRR—SufRff'AR— u(d} + ©) = aconstant. . . (3.12). 


As fis a function of R only and © of (r/z) equation (3.12) can only be satisfied if ®} + ©; and 
(D2/r?) are functions of R only. From (3.11) we get 
®: 5 2 4 D? ax, 
——6® a RI 3.0). 
ij Ta (z+ R)? a 
If (r/R2)? and higher powers can be neglected, as has been already assumed, the right hand side 
of (3.13) can be approximately taken equal to (D?/R?)?. Equation (3.12) now reduces to 


4 D® 
A 2 —2R? a —8u[fFaR—8u[Rff'dAR— A-+u) Ze constant (3.14). 
Differentiating (3.14) with respect to R? and putting R? = x we get 
2% 
2A +2 W2f"rAHIW)P=-20 WG  @:15)- 
Putting ce = 2u/(A + 2u), 0< e< 1, and integrating it we get 
“ D: D? (2 — ec) SR 
= Sr R d—-g)2 2) 
whose solution is 
Me 
2 i plc 37) 7 
= T = («= — DM)? —D, tan! & D ) | te.) 
Bu: : 
where D, and D, are constants of integration and 
D? 2 —c 
mn A I Er De Br Lo 
0 
The components of displacement now become 
N ee e 1 
nr any; | 
v0, dir . (3.19). 
1 9 » 2 
; . R21—9 — .D3)2 | 
w= 2— en Gun 2—D,) =D, tan? FR ı) | — D, 
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| These show that plane faces of the plate get bent into 
A spherical surfaces and lines parallel to its axis go into 
N normals to the corresponding surfaces. Also diame- 
Lu ters of the plate get bent into great circles of the 
| N surface, but any circle in a plane of the plate goes into 

a circle with its plane perpendicular to the axis of the 
| surface. In Fig. 1 strained lines are shown by thin and 

ni dotted lines and the corresponding unstrained ones by 

f | we thick lines. For example, EAF goes into E’A'F'. 


OR If we take the inner and outer radii-of the 
N spherical surface to be a and b respectively and 
N assume that the point A,z=R=ar=0) does 
N not undergo any displacement, we get the value of 
D28s 


1 


2D, [ru u Bi D>3)2 


6 


1 
De ea Dı 
ee z D)) 2.8.20), 


D, and D, are determined from the boundary: con- 
ditions. 

The decrease in the thickness 2h of the plate can be determined if a and b are known. 
From (3.19) we get the value of2 Ain terms of aand bas 


Fig. 1 


1 1 
3D 2 . _. Da? Da? 
le Bu + D,tan-ı a_ np?) 
1 1 
(3.21) 


In like manner if the radus of the plate be given by r, = a, then the corresponding circles 
in the strained state are given by 


PR 4aqaD Ei um 4:05.D 
ıTg+4D "= g+4D Br 
and their z-distances are given by 
rat =D) b (a — 4 D?) 
2 aa Fa De and a LAD: BETEN 


4. Boundary Conditions 


At any point of the boundaries of the spherical shell given by R= aand R = btthe compo- 
nents of tractions in cylindrical coordinates are 


RN,=rAAFB=EUR 9) - AD, (D, 4 2D, Se 
RN,=2QA+#W—4AuR RB eo Zoe 
R Ne = 0 . as je lu 10 el Da ar Le ee er en Er SE N . (4.3). 
With the help of (3.8) and (3.16) we can rewrite them as 
RN IFRNGL 4DE 2(3—ec) D: 
N = Die) 0 s g 
r z c E ) R2e " 1c m) er 
If the inner and outer surfaces are free from traction we get 
a 4DE 2(3—c) D: 
3— 20) ——t +- .- 
( & a2° je a (4.5), 
r 4DE 2(@8—c) D: 
ee < 
( c) ET ar ee, 0 ee 


Subtracting one from the other these give 


np 2 .2—C 212g) 
DI Dzar: 


0 


ba az® 


2(1—c). 
ii Ton 
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17®— 20) (#?’— a?) 
DZ 4 "PRü-9_ Rü-g . . . . ee ee Bar Eon 0 ee (4.8), 
11—c0@—2e) b?e — a?° 
PL et a 
D 5 Be) (ab nn a ee er (4.9) 
ı : 1 1 
For Poisson’s ratio o = E ( = 2) these constants take the simple form. 
12 ab ji 2 
D: m — — ! 2 = — 2 zz — 
r Marb’ D: bt, | D eb Be A (4.10). 
The tractions over the edge tan! (r/z) = constant are given by 
By nlEAse) FED RB Rn (4.11), 
ee ee PIERTEREN (4.19) 
Rye = 0 


With the help of (3.8), (3.13) and (3.16) these become 
Ruin 3emD> W222 772.D7@ 6) D 
=. D7 R Re ÜZOIR: tu Hra—T (4.13). 
This shows the resultant traction is T, applied at right angles to the edge. 
Its resultant for unit length of the edge is 


zZ g 


F= sr RdR— 27 en GR Do a) > (@— 0) Dilog — Di a0 — @a-0} 


Putting the values of D, and D from (4.8) and (4.9) this takes the form ve 
En i = (do) 2 ht AI eu as en Er Des og”... dB). 

For c = its value is 
N 2 wab 1og Mae N ee (4.16). 


41 


The traction T also gives a couple of moment 


2 u 1 3 2(1— c)D, 5 
ey) 3 3 9) 2 A 
M= 1 R?dR= alt 2c) (b’—.a?) 5, @ c) D?(b—.a) Bang ° (b3- a ) 
(4.17), 
which for c = n reduces to 
1 9 
M=2u(b — a) 6 (b — a)? + 11 a Ne (4.18). 


Thus we require a resultant force F and a couple of moment M to keep the plate bent in a 
spherical surface of inner radius a and outer radius db, the curved surfaces R=aand R=b 
being free from tractions. 


For a thin plate of strained thickness e we canputb =a -+ eand getfor ce = 5 5 
F=-lnas, M=- DANS. WO IAnTLL Sargek: (4.19). 
For rubber o = 0.49, c = 0.02, and we get the following approximate values 
D} — m D=c a IE RER Een (4.20) 
D: — ne E “= log - es Wrtane (4.21) 
F=6u “—_ 108 | M=6u Br Iog Arch Hat a (42). 
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The w-displacement is 1 2 Be ua re 


(i) Let us take the case of ce — re a=3,b=4. The constants are given by 


1 1 


D= n (b-+.a) (2 + ae) — 9,7965, D, = 3.130. 


3 3 
a, = b?—a? ” 
Die nal eu D-42193 
The original thickness 2 h of the plate is found to be 
2h=1:202 


Decrease in thickness — 17. per cent. 


1 
(ii) For the case c = 3° nel werget 


D=2, rD= 1%; :D 6-5 Dand 72h —2, 737 

In this case the decrease in thickness becomes 27 per cent. 

(ii) By taking a = 11, b = 11.8 we find a similar decrease. 

The force and couple to be applied to the edge in (ii) are found to be 

M=39.10% N 

The corresponding change in the radii at the edges can be found from (3.22). If a, = 6 the new 
radii is almost 12. 

For the incompressible case c = 0 


3 
DU, Do, Do, 


and 
2h='1.732.(0 4). 

In this case the decrease in thickness shoots up to 42 percent. The force and the couple to 
be applied at the edges are given by (4.22). An additional radial force has now to be applied to 
the curved surface of amount 

— 6 u (b? log b — a? log a) 


The u-displacement becomes r and hence such a state is not physically possible. 
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KLEINE MITTEILUNGEN 


On the eigenvalue Latzko's differential combination of the first three (odd) Legendre polyno- 


equation mials. More recently Durfee [2] obtained more accu- 

In the study of the problem of heat conduction in rate values by direct series expansion of the solution. 
fully developed turbulent flow in a cireular tube, In this nots, an alternative approach is given which 
Latzko [1] arrived at the differential equation gives good estimates to the eigenvalues and eigenfunc- 


tions with much fewer terms than does the expansion 


d (1 2#) N way. . (1) in Legendre polynomials. 

dx dx We first consider the solution of an auxiliary equation 
ver, sel nee che none which is similar in appearance, although not identical, 
distance, and y is proportional to the radial distribution Bihaan: 
of temperature. The eigenvalue » is a parameter which d 1 „du 5 
involves several physical constants, and is proportional dx Ve de An 


to the exponential decay of temperature in the axial 
direction. The admissible values of w are determined 
by requiring that y(0) = 0 and that y(1) be finite. 
 Latzko[l] abtained crude approximationsto the 1 ey, (6 13 \dy_ A 

first three eigenvalues by expressing y as a linear dt E 7 :)% | 497 ze) 


The change of variable x” = t converts this equation 
to the canonical form of the hypergeometric equation: 
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with the general solution 
y= AFla,b,c,t) + BU F(a—c+HL1l,) 


4), 
b—c+1l, 2-—.c,|i) | ir 
a+b =6/7 
ab = —1/49 
et 


The condition y(0) = 0 requires that A = 0, so that 
y= BiÜTF(a+1/7,b + 1/7, 8/7,t) .. (5). 
Further since 
a+b—c=0 
it follows that, at? = 1, the series will diverge ([3] Ch. 2; 
Art. 1) unless ist terminates, and this requires that 
either (@ + 1/7) or (b + 1/7) be a negative integer or 
zero. Thus admissible solutions occur when 
a+17=—n, (6) 
b+17=n+ 8/7 FE 


which gives: 
An=Tn+1)(Tn-+]), | a 
DE A € 
The resulting polynomial solutions are known as the 
Jacobi polynomials, [3] and are defined by 


In, B,2)=Fl—n,a+n,ß,2) .. . (8). 
Thus the eigenfunctions of Eg. (2) are 
Gn(x) = Pn% En (8/7, 8/7, 2?) = Ba fnl&) » » (9) 


where D„ is an arbitrary constant. The functions f, 
are orthogonal with respect to x° as weight factor: 


1 
Narr.e)n(o)dz = 07 mn, ! 
0 


as /Ao) dr so 

9 en en 
(l4n+8)[1-8-15-—-(7n+1)] n 

For application to Eg. (1) it is convenient to define a 

new set of functions ®,„(x) such that 


On Yin 
so that 
1 
= if N 
An il > Dr) Dn(x) de — Ömn |Omn _ 5 if n z ii ; 
) 
Assuming a solution of the form 
N e 
Ve Da RIERLTDN 
n=1 
and substituting in Egq. (1) we find 
N 
d d®, 
> X Hwa? =) 13 
„en Zt z) -oaT®d, (13) 
n= 


and because the D, satisfy Egq. (2), this becomes 
N 
>. 01.22. 2,@) = 02 d,K)l =U. (14). 
n=1 


We may now apply the Galerkin method to deter- 
mine the a, by requiring that the left side of Eq. 14 be 
orthogonal to the ®&, frn=0,1,...,N. Thus we 
arrive at he following system of algebraic equations 


N 
> an [dömn — © Amn) = 9 N ro 


n=0 


il 
where Ayn = S& Dun) D„(x) de. The characteristiec 
() 
equation of the system (15) is 


\dmn — ® Amal = 0 . (16) 


and the roots give approximations to the first N eigen- 
values. The @’s may then be determined (to within a 
numerical factor) to give the corresponding approxi- 
mations to the eigenfunctions. The system (15) is 
suitable for treatment by matrix methods. 

Listed below are the approximations obtained by 
using respectively 1, 2, and 3 terms: 


One term: 
An = .11429, 
©, = 8.75, 
y®=6®, 
Two terms: 
Ar = Agı = — .0055742, 


As, = .006644 , 

©) —= 8.728, y® =D, — .0516 B,, 

@, = 157.3,y®2 = .05168, + ®,. 
Three terms 


Ars = As = — .0002722, 

Ay; = Ag = — .0008766 , 

Ass = .0023512, 
o= 8.172798, y®’ = ©, — .05163 ®, — .002022 d, , 
@, = 152.8, y”? = .05120 ©, + BD, — .21236 O,, 
0; = 462.5, y® = .01225 ®, + 2117308. + ®,. 
Durfee (3) obtains the values 

©, = 8.127147, w,152.4, ©; — 435.06 . 


Comparison of these with the values obtained above 
by the three term approximation indicates that the 
use of Jacobi polynomials represents a substantial 
improvement over that of Latzko employing Legendre 
polynomials, the latter giving the following approx- 
imations: 


BSH 
© —= 164.4 , 
@, = 1700. 
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Zur Berechnung zylindrischer Behälter mit 
quadratisch veränderlicher Wanddicke 


Dieser Fall, dessen Lösung durch einfache Funktionen 
darstellbar ist, wurde bereits von K. Federhofer!) 
behandelt. Er ermittelt die Bedingung, mit deren Hilfe 
man die Differentialgleichung des Problems 


F 
eg) =) a 
in zwei gleichwertige Differentialgleichungen zweiter 
Ordnung — mit voneinander unabhängigen Lösungen — 
aufspalten kann. Die Zerfallsbedingung wird mit der 
Zylinderwanddicke 

y\2 

4 

= (z) ea ae) 
ı 


1) K. Federhofer: Anzeiger der math.-naturw. Klasse der Öster- 
reichischen Akademie der Wissenschaften, Jahrgang 1950, Nr. 11 
(Seite 275 bis 237) 


erreicht. 
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Neben diesem Weg existiert nun noch die Möglichkeit, 
ohne Zerfall die Lösung der Differentialgleichung (1) 
zu finden. Führt man nämlich die Differentiation der 
Gleichung (1) aus, dann erhält man 


+10) 0’ +20 +M)-E")=0 


h=h(x) : Wanddicke, 
b : mittlerer Rohrhalbmesser, 
oe : radiale Verschiebung, 
e', o'’ usw. : Ableitungen von o nach z, 


ZT „ er : Faktor, der die Poissonsche Kon- 
un m ) stante enthält 


Man erkennt eine Reihe, die aus zunehmenden 
Ableitungen von o besteht. Die Koeffizienten sind Aus- 
drücke von h mit abnehmenden Ableitungen. Diese 
Form legt nahe, eine Gleichung zu suchen, welche dem 
Typ der Eulerschen Differentialgleichung entspricht: 


a RM a RTL. RD... +00 =0 (4) 
Man erreicht es mit dem Ansatz für } nach Gleichung (2). 
Nunmehr lautet Gleichung (3): 
F.xi : 
Meg + BO: ag" +12: 0-0" Hate” J=0 (6) 
und die hierfür bekannte Lösung: 

o= A-4404 Aa A, ade. . (6) 
mit der charakteristischen Gleichung: 


F.a1 
MORE R— 24-0: . M. 
2. 


Damit hat man die allgemeine Lösung direkt aus der 
Differentialgleichung 4. Ordnung gewonnen. 
Substituiert man noch 


so bekommt man anstelle der Gleichung (7) eine bi- 
quadratische Gleichung: 


17 15\2. F.:x4 
Arte ey? 1 
u gu (7) He O9) 


und damit die bereits von K. Federhofer gefundene 
Lösung: 


0o= A: . cos [I(A) - In (B x)] 
+0.) . cos [ZA)- In(D.x)]. 


Die Zerfallsbedingung wird bei dieser Darstellung 
gewissermaßen erst nachträglich in der charakteristi- 
schen Gleichung untersucht. 


Dresden. 


W. Vocke. 


Nebensymmetrische Matrizen 


In der Praxis treten bei den verschiedensten Pro- 
blemen Matrizen auf, die zur Nebendiagonalen sym- 
metrisch sind, so z. B. bei den Gleichungssystemen zur 
Berechnung der Stabkräfte von Fachwerken mit zy- 


klischer Symmetrie [1], [2], aber auch bei den Glei- 
chungssystemen zur näherungsweisen Behandlung von 
Eigenwertaufgaben nach dem Differenzenverfahren bei 
Differentialgleichungen oder auch Integralgleichungen 
mit Differenzkern. Im folgenden sollen daher die 
Eigenschaften solcher nebensymmetrischen Matrizen 
untersucht und durch das zugehörige Eigenwertproblem 
charakterisiert werden. 

Zur Darstellung der Nebensymmetrie quadratischer 
Matrizen führen wir eine quadratische Matrix © ein, 
die in der Nebendiagonalen mit den Elementen I und 
sonst mit Nullen besetzt ist: 


Hanne 
= nam: el ege .. ). 
De 
Wegen 


N 


| ER 


= (re 
und © reell gilt 4 
oS=6=- en nern (2). 
Bilden wir das Produkt der n-reihigen Matrix © mit einer 
ebenfalls n-reihigen quadratischen Matrix ® = (b,£) 
CHE (d.,n +1 —$) (daB) 


; 
-(Z dans ir) Guru ®) 


so erkennen wir, daß eine Linksmultiplikation won © 
an ® die Umkehrung der Zeilenreihenfolge bewirkt. 

Entsprechend liefert die Rechtsmultiplikation von 
San 


2.0 (daB) (dx, n+1 er) 


e 
m ( 3,52» on + en r (&.,n +1-,) (4) 


die Umkehrung der Spaltenreihenfolge. 
Liegt eine Matrix ® vor, die zeilensymmetrisch ist, 


= =D. 12 ae 


so gilt 

BEshBT, Ta (5) 
Entsprechend besagt 

VD Sr, re „ (6); 


daß die Matrix 


Due 
spaltensymmetrisch ist. 
Bilden wir nun aus 


SB nr 
die Transponierte 
Scy”ce=- user‘ Ag persiy, Bi 


so erkennt man, daß Nebensymmetrie gerade dann vor- 
liegt, wenn gilt 

N =. © BIC ee 
und eine Matrix & nebenhermitesch ist, wenn die Be- 
ziehung besteht 

GIS EISEN ae) 
Wir verfolgen die Eigenwertaufgabe für die allge- 
meinere Matrix mit 


ARTE; 
Wir entnehmen diesem Zusammenhange sogleich, daß 
X und W’ zueinander ähnlich sind. Ähnliche Matrizen 
besitzen aber nach einem bekannten Satz gleiche Eigen- 
werte. Da nun die Eigenwerte von W gleich den kon- 
jugiert komplexen von X sind, die sich mit denen von 
A decken müsen, folgt daraus, daß U entweder reelle 


Eigenwerte besitzt oder, wenn einzelne komplex sind, 
gleichzeitig die konjugiert komplexen auftreten müssen. 


B=W = R-1.. 10). 


FERNE RD REITEN 
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Auch die Eigenvektoren von X und W’ hängen eng 
ee Schreiben wir die Eigenwertaufgabe in der 
orm ; 


WA er. Ad en il), 


mit dem n,-fachen Eigenwert A, und p, dazugehörigen 
linear unabhängigen Figenvektoren tx, So gilt wegen 


(10) auch 
PU Pr, eh les 
und nach Multiplikation mit ® von links 


Apr, )=ARr,) RAN 


woraus neben der schon nachgewiesenen Gleichheit der 


Menge der Eigenwerte von X und W’ noch folgt, daß 
die Eigenvektoren ),, von W gegeben sind durch 


Da ne en ee 


nn in der gleichen Anzahl wie die r, „ linear unabhängig 
sind. 

Unsere Aufgabe läßt sich mit einer allgemeinen Figen- 
wertaufgabe in Verbindung bringen, indem wir (11) mit 
® multiplizieren. Wir erhalten 

PAr.=4, Br. 
mit zwei Hermiteschen Matrizen BA und P; denn es 


gilt 

TAYI=- I P=ETN) ud VW: 
Da hierin die Matrix ® offensichtlich indefinit ist, lassen 
sich weitere Aussagen über Eigenwerte und Eigen- 
vektoren nur bei zusätzlichen Voraussetzungen über die 
Matrix ® X machen [3]. _ 

Eine Zusatzforderung an 9, welche die Transformier- 
barkeit auf Diagonalgestalt zur Folge hat, ist nahe- 
liegend und soll nun betrachtet werden. Schreiben wir, 
(10) in der Form 


HET Tr anemiehl 
so ist wegen 
PA-AT-RI PP 
die Matrix W dargestellt als Produkt zweier Hermite- 
scher Matrizen und wird dann als symmetrisierbar 


bezeichnet, wenn eine dieser Matrizen positiv definit 
ist. Da im vorliegenden Falle ® indefinit ist, soll für 


das folgende die Matrix BW’ —= AU%® als positiv definit 
angenommen werden. Eine Ähnlichkeitstransformation 
auf Diagonalgestalt ist dann möglich [4]. 

Dazu bestimmen wir in bekannter Weise die positiven 


Eigenwerte u der Matrix BW = A ® und ihre nor- 
mierten Eigenvektoren, die zu einer unitären Matrix U 
zusammengefaßt werden, so daß gilt: 


kı 0 
Banner ten 
0 2 
Wir bilden sodann eine neue positiv definite Hermite- 
sche Matrix ® 


Wal _ 
Bu IVe|. 
Ben 
deren Quadrat ® W liefert: 
KEEUNMNN)MWIW)-UYV-UDTV=- RW. 
Nun formen wir die rechte Seite von (10) um 
A-(PAW)P-RWETRWT, 
so daß folgt 


V-UdW=-W, 


TEE 
Die Matrix auf der rechten Seite ist nun eine zu A 
ähnliche Hermitesche Matrix, die durch eine unitäre 
Transformation ® ihrer normierten Eigenvektoren auf 
Diagonalgestalt gebracht werden kann, mit ihren reellen 
Eigenwerten, die zugleich die von Y sind, in der Haupt- 
diagonalen. 


Durch die Ähnlichkeitstransformation 
Zu 0 
Faıuaa-| % 
N 
\0 An 
läßt sich demnach % auf Normalform transformieren. 


Damit haben wir folgendes Ergebnis gefunden: 
Eine Matrix YV, für die 


A-BUH B-P- 9 
gilt und darüberhinaus die Matrix 
PV-AR 
positiv definit ist, besitzt als symmetrisierbare Matrix 
reelle Eigenwerte und lineare Elementarteiler. 
Nebensymmetrische Matrizen mit A=4X, BP=6& 
und © positiv definit, sowie nebenhermitesche Matri- 
zen mit B— © und A © positiv definit, gehören zu 
dieser Klasse. : 
Beispiel 1: 


ei h weh 
A=([0 1; AS=65%=[i 10). 
2.001 110 2 


A ist nebenhermitesch wegen SAYS =W. US ist 
indefinit. Die charakteristische Gleichung besitzt zwar 
reelle Koeffizienten, jedoch ist die Matrix nicht symme- 
trisierbar. 

Beispiel 2: 


zus 12.050 
=11720);, AS E70 251% 
171.20 0 1-1 


U ist reell nebensymmetrisch wegen SAYS = W und 
U © ist positiv definit. Daher ist A symmetrisierbar, 
besitzt reelle Eigenwerte mit linearen Elementar- 
teilern und es gilt 


12072070, 0.210 ,/12500 
RS IE O EIN 327201 KO Te 
0 1—1/\110/\010 


ae) 
Ze 
le] 


Hannover. K.Jaeckel. 
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Integraltransformationen mit Differenzkern, 
bei denen Kern-, Objekt- und Bildfunktion zum 
gleichen Typus gehören 


Die Gaußtransformation 
[0,0] 
5 II ee 
Al) == Be y(£) dE 
Z 
—& 


hat die Eigenschaft, daß eine Gaußfunktion 


hı en NE — 0)? 


Vr 


wieder in eine Gau ßfunktion 


y(&) = 


Y(x) eu hz „hi —)! 
Vr 
überführt wird, wobei gilt 
1 1 1 
m mtn 
Diese beiden Funktionen und die Kernfunktion 


besitzen den Flächeninhalt 1. Obige Transformation 
hat für die praktisch wichtige Frage der Zerlegung einer 


26 
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gegebenen Häufigkeitsverteilung in eine Summe von 
Einzelverteilungen Bedeutung. 
Ist nämlich 


n hi, —h,(—a,)? 
en a ee ’ 
y(8) 2, Zr 
so folgt 
N Nas — 13 ,(@— %,)? 
Y(x) = > Fi e ’ 
v=1 Yr 


wobei die Maximumabszissen «, der Einzelverteilungen 
gleich und die Streuungszahlen h,, < h,, sind. 

Gelingt es demnach, die Integraltransformation um- 
zukehren, so wird damit eine Aufsteilung der Gauß- 
funktionen erreicht und in geeigneten Fällen eine Auf- 
spaltung der Mischverteilung in Einzelverteilungen er- 
möglicht [1]. 

Die numerische Durchführung der Umkehrung ist 
indessen außerordentlich schwierig und fehlerempfind- 
lich, so daß es naheliegt, nach anderen Integraltrans- 
formationen zu suchen, die für eine Aufsteilung in 
Betracht kommen und unter Umständen numerisch 
einfacher zu behandeln sind. Hier soll nur der erste 
Teil der Aufgabe dargestellt werden. Wir fragen all- 
gemeiner nach denjenigen nicht negativen Funktionen 
K(«), die im Intervall (— 00, + 00) quadratisch inte- 
grabel sind mit dem Flächeninhalt 1, so daß mit end- 
lichen positiven Parametern h,, h, und H gilt: 


1, Kiki —a)) = | H K(He — 8) 1 Kly(E — 2) dE 
n (). 


Zur Bestimmung solcher Funktionen multiplizieren 


ei ul® — 0) 


1 
wir die vorangehende Gleichung mit es, und 
27 


integrieren über x von — oo bis + oo [2]. Auf der linken 
Seite erhalten wir 


oo 
1 fe F 
ehe h, K(h,(x — o)) WR) ge 
V2x | 5 x ) 
—o 


oder mit h,(@ — «) = o 
eo [2 


. i—o 
L = —— na do 
— 00 


und auf der rechten Seite nach erlaubter Vertauschung 
der Integrationsfolge 


oo 


1 ; r 
ne h, K{h,(E — a)) ET) as 


V2r 


are.) 
oo 


x | H R(H(&«— &)) @T9 de 
RO: 


oder nach Einführung der neuen Integrationsvariablen 


h&e—o)=0o; H«x—&)=r 
schließlich 
1 rag i Be 
Rena K(o)e 4 do 
V2r 


==i0Q0 


>) 
= 1 h : 
K (u) > | Kiez) e 72% (2) 
Zt 
u) 


können wir das Ergebnis in der folgenden Form schrei- 


ben 
u 
VErk, (%) - ar K, (") :PrK, (z) 
2 1 
und gewinnen mit 


In (V2r K,(u)) — q(u) 


RE 


Setzen wir darin e 
u h, ha, Zen 
—v; ——= 0; —— A >= 0 , 

h, v; hı 4 = H ( ) 


so erhalten wir schließlich die folgende Funktional- 
gleichung 


zunächst 


gw) = gAv) + glp(A)o) ... - - (3), 


aus der zu entnehmen ist, daß p von / abhängen muß. 
Aus dieser Gleichung sind zwei Funktionen zu bestim- 
men, nämlich g(v) und p(A), die wir als zweimal stetig 
differentiierbar voraussetzen. 

Zunächst stellen wir fest, daß g(0) = 0 sein muß, 
und können für die weitere Rechnung v # 0 annehmen. 
Die Funktionalgleichung (3) differentiieren wir nach v 
und nach A: 

yw)=Ag(Av)+pQ@)g(pv), 
0=vglAr) Hop(A)g’(pv): 

In der zweiten Gleichung können wir durch v kürzen 
und p’(A)# 0 annehmen. Denn wäre 9°(),) =0 für 
irgend eine Stelle A, > 0, so wäre nach der zweiten 
Gleichung g’(A, v) = 0, also überall g = const, so daß 
diese Möglichkeit wegen g(0) = 0 nur die triviale Lösung 
lieferte. Daher können wir aus den vorangehenden 
Gleichungen g’(pv) eliminieren und erhalten: 


PA)gW)=AP—pgUAr)... . (4). 
Hierin dürfen wir Ap’—p auf Werte ungleich Null be- 
schränken; denn wäre für eine bestimmte Stelle A, > 0 
die rechte Seite der vorangehenden Gleichung Null, so 


2 on £0 mit q’(v) =0 ein Fall 


1 
vor, der wieder nur auf die triviale Lösung führte. 

Differentiieren wir die Gleichung (4) nach v und 
nach A 


Pa)g’W) =AAp —p)g’Av), 

PA’W)=Ap’dAv)+vAp —p)g(Av), 
so läßt sich aus den letzten beiden Gleichungen q’’(A v) 
eliminieren mit dem Ergebnis 

pad) —ApPAgwW)=—ARp’(A)g(Ar) (8), 
aus dem unter Heranziehung der Gleichung (4) nun 
auch g’(A v) entfernt werden kann. Wir erhalten dann 

AP —popgd’+iAprg=0. 

Diese Gleichung läßt sich so umformen, daß die linke 
Seite nur von v, die rechte nur von A abhängt, wobei wir 


für nichttriviale Lösungen g’(v) #0 für jede Stelle 
v # 0 voraussetzen können: 


vg’w) _ __APp(A)p’(A) 
d’(e) ApPr—p)p ve 
Hieraus ergeben sich wegen der Unabhängigkeit der 
beiden Variablen die Differentialgleichungen 
vw _ A p(A) p (A) 
ı(w) men ae 
mit den Lösungen 
)=Auc+B; K=aß+Bßzeh). 
Damit für v = 0 der Grenzwert von q(v) mit dem 
Funktionswert g(0) = 0 übereinstimmt, muß B=0 


und ce > 0 gewählt werden, womit aus (3) auch « und ß 
bestimmt sind. Wir finden schließlich 


W)= Ar; pPUV)+Ae=1, c>0, 


läge wegen p’(A,) 


ER LEN! 


WETTER 
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und damit nach Gl. (2) die Fouriertransformierten 
1 edv 
V2r 


und durch das Umkehrintegral die gesuchten Funktio- 
nen 


K,o))= 


[oe] 
K(«) ZB Ar eietn. 


In der letzten eg werden nur solche Werte 
A und c zugelassen, für die das Integral existiert und 
eine quadratisch integrable Funktion K(x) liefert. 


Diese Bedingungen sind erfüllt, wenn ei” selbst eine 
quadratisch integrable Funktion ist. 

Wegen g(0) = 0 und damit /2rK,(0) =1 gilt, daß 
für jede dieser Funktionen der Flächeninhalt von 
selbst gleich 1 wird. Beschränken wir uns auf die folgen- 
den symmetrischen Ne 


And e *ll®. eier 
A [- LOB), 


— 00 
so finden wir mit a = = und c=2die Gaußfunktion 
IK) er > 
Yr 
und für a = 1; c = 1 die Funktion 
Be . 
al+x 
die aus der Potentialfunktion 
al 
ae+y 


mit y = 1 entsteht. Durch diesen Zusammenhang er- 
geben sich Möglichkeiten, die erwähnte Zerlegungs- 
aufgabe potential- oder strömungstheoretisch aufzu- 
fassen. 

Für andere ganzzahlige Werte c ist die Integration 
von (7) nicht mehr elementar durchzuführen, es ist aber 
möglich, für die gesuchten Funktionen Differential- 
gleichungen anzugeben. 

Für nichtganzzahlige Werte c > 0 ist auch dieser 
Weg nicht mehr gangbar, die Funktionen K(x) lassen 
sich dann durch Reihen mit Hilfe der Hermiteschen 
Polynome darstellen. 


Hannover. K. Jaeckel. 
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Hauptachsentransformation der quadratischen 
Form für die Streuung 


Die wichtigsten statistischen Kennzahlen einer Reihe 
von stochastisch unabhängigen Meßwerten x; ? = 1, 
2,...., n) sind Mittelwert x und Streuung s, die nach 
folgender Vorschrift gebildet werden: 

N N 
na Den Ne HE, rer ll): 
i=1 

Nimmt man an, daß die x; aus der gleichen normal 
verteilten Grundgesamtheit stammen, so läßt sich 
nachweisen, daß © und s? selbst wieder stochastisch 
unabhängig sind und einer Normalverteilung bzw. 
einer y2-Verteilung genügen. Dieser Zusammenhang 
läßt sich durch Wahrscheinlichkeitsbetrachtungen un- 
mittelbar im Merkmalsraum der x; [1] oder auf ganz 
anderem Wege durch Transformation der quadra- 
tischen Form für s? auffinden [2]. Hier wollen wir die 
Hauptachsentransformation dieser Form in den Vor- 


dergrund stellen und sie mit Matrizen explizit durch- 
führen. 
Wir betrachten die quadratische Form 


N 
x Rare .. &n) = > (m n> 2; 
nel iz 5 
und Bi sie mit 
I—n 1 A 1 x 
| min 2 selele 
1 1 1_—n Ca 
in der Gestalt 
HNO EN. Mi a (4). 


Die Matrix X ist reell symmetrisch, so daß sie durch 
eine reelle orthogonale Matrix DO auf Diagonalform 
gebracht werden kann, mit den Eigenwerten A, der 
Matrix W in der Hauptdiagonalen. 

Infolge der Variablentransformation 


ee reed) 
Y=(y Yo... Yn); DD MDEE 
rır=yVoy=Yy) 


ergibt sich für die quadratische Form 
— rd=yYyV/ADYy=Ydy mu 
+AyYH+ u (6). 
Die Eigenwerte A, der Matrix W und die normierten 
Eigenvektoren r,, aus denen die orthogonale Trans- 
formationsmatrix aufgebaut wird, lassen sich aus dem 
folgenden homogenen Gleichungssystem bestimmen 


mit 


und 


In 1 Be 1 u 
En 
1 1 1—n— A) \%n 


dessen n-reihige Determinante verschwinden muß. 

Nach Zeilenvertauschungen erkennt man die vor- 
liegende Determinante als eine zyklische [5], für die 
sich die Eigenwerte leicht berechnen lassen: 

1 0R einfach, 5 
een, n—Nfach”) =E (8). 

Für A,=0 finden wir aus dem Gleichungssystem (7) 
den zugehörigen Lösungsvektor zu 


I end Tee ie 
während wir für A2,..,n=—n nur entnehmen 
können, daß die Komponentensummen Null ergeben 
müssen. Jeder Vektor von der Form 

—1 
on (# es ser 2, s) 
i= 
ist Lösungsvektor. Die Eigenvektoren, die wir durch 
Orthogonalisierung finden, bilden die Spalten der 
orthogonalen Matrix 
Ne‘ 
Yn 


1 
ae Kar ;E 
1 
9 IR 
el m 
Sg ı 
2020005321 1 
1 
le — 


y 1 1) (2) 


N u a) 


| 
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wobei die normierenden Faktoren an entsprechender 
Stelle in der nachfolgenden Diagonalmatrix angeordnet 
sind. 

Mit der Matrix O lautet das lineare System der 
orthogonalen Transformation 9 = O’r ausgeschrieben: 


1 
atnt ten: 


%= 78 (& + %,— 2%) (10). 


yo er Darm oa to er ORT ONE OR TOERIFTHLZT 


Wir entnehmen daraus insbesondere den Zusammen- 
hang: 


N N N 
Iu= Yanyı=nı Zuy=Ny.. U). 
i=1 i=1 k=1 

Von Bedeutung ist nun die Feststellung, daß die neuen 
Variablen y7; untereinander stochastisch unabhängig 
sind, wenn es die Ausgangsgrößen x; sind, sofern sie 
als normalverteilt angenommen werden. 

Für den Nachweis nehmen wir also an, daß die x; 
der gleichen normal verteilten Grundgesamtheit ent- 
nommen sind, so daß mit dem theoretischen Mittel- 
wert ‚ı und der theoretischen Streuungszahl h gilt: 


// —h? (2, — u)? 
ulx) = an : 
Yr 
Dann finden wir die Verteilung w(y,, Y3 ---; %n) da- 


durch, daß wir in zugeordneten Räumen ()) und (r) 
die Wahrscheinlichkeiten gleichsetzen: 


. (12). 


Yırdyı YmtAym 


| .. [ wYy»Ys--- 
fr , R n \n 1 an? 
— He: Ir er a en 
IT 


R 


,Yn) dyı dys » - » 4Yn 


wobei das Integral in (x) zu erstrecken ist über den- 
jenigen Raumteil, der infolge r = Oh dem Rechtflach 
in (y) entspricht. Die Auswertung des rechten Inte- 
grals wird sehr einfach, wenn man neue Variable ein- 
führt, und zwar gerade diejenigen, die durch y = Or 
mit der Funktionaldeterminante + 1 gegeben sind, so 
daß dann die Integrationsgrenzen links und rechts 
übereinstimmen und zudem beliebig sind. Der unter 
dem rechten Integral auftretende Exponent läßt sich 
mit (11) durch die neuen Variablen einfach aus- 
drücken: 


N n N 

ln 2 = 08 ‘ 1 2 
Mi —u=- 3 -2u3u tun 
i— = 1 


1 i= 


Vz. n = n 
7 2 y—2ulny+wn= (y— Ir w? + 3 yr- 
E k=2 


Daraus ergibt sich 


w(Yy;, Du PAR nn m): H u —h2 vn 


Vr k=2 


in Produktform als Kennzeichen .der stochastischen 
Unabhängigkeit der Variablen y; und es gilt: 


h mn) 


w(yı) = Yz 


en 
DE 


Vr 
Betrachten wir nun die auf Hauptachsen transformierte 
quadratische Form (6) und beachten (11), wonach yı 
dem Mittelwert proportional ist, weiterhin aber, daß 
in (6) das erste Glied mit y, wegen A, — 0 gar nicht auf- 
tritt, sondern die Form nur noch von den übrigen 
stochastisch unabhängigen Größen y7; abhängt, 


nD#-90-&1 . ide 


so folgt daraus bereits, daß die Streuung s vom Mittel 
wert X stochastisch unabhängig ist. Die Verteilung von 
s? als Summe von Quadraten von n— 1 normalver- 
teilten Variablen ergibt sich als %2-Verteilung von der 
Ordnung n — 1." - 

Vou der bemerkenswerten Eigenschaft der neuen 
Variablen y,, stochastisch nnabhängig zu sein, wenn 
von stochastisch unabhängigen normalverteilten Vari- 
ablen x; ausgegangen wird, läßt sch noch weiterhin 
Gebrauch machen. Während die erste Formel (10) die 
Grundlage für den t-Test darstellt, liefert die zweite 
und die letzte Gleichung (10) den Zugang zu Prüf- 
größen, die für die Behandlung des Ausreißerproblems 
von Bedeutung sind [3], [£]. 
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Verteilungen von Warte- und Anschlußzeiten. 


Für gewisse Warte- und Anschlußzeiten ergeben sich 
nichtnormale Verteilungen, die im folgenden näher 
betrachtet werden sollen. Es zeigt sich, daß ein Prüf- 
test für Wartezeiten noch mit den üblichen Begriffs- 
bildungen entwickelt werden kann. Hingegen wird die 
Konstruktion eines regulären Testes für Anschluß- 
zeiten so schwierig, daß in diesem Falle ein Rangtest 
angebracht erscheint [1]. 

Wir skizzieren das vorliegende Problem an Hand 
eines praktischen Beispiels [2]: 

Zwei kleine Vorschmelzöfen geben ihr Schmelzgut 
zur Weiterbehandlung an einen dritten großen Schmelz- 
ofen ab. Die Lieferzeiten x des einen und y des anderen 
kleinen Ofens sind wegen der unvermeidlichen Streu- 
ung verschieden; ihre Verteilungen mögen gleich und 
normal sein, so daß gilt: 


h 2(r 2 2 2 
u(x) = Yr Fr ’ v(y) ER Vz En N a (1) 
i TC 


mit z dem theoretischen Mittelwert und % der theore- 
tischen Streuungszahl. Da der große Schmelzofen erst 
dann angefacht werden kann, wenn das Schmelzgut 
des später anliefernden kleinen Ofens vorliegt, entstehen 
Wartezeiten z = |® — y|. 

Zunächst ist die Frage zu beantworten, welcher Ver- 
teilung die Wartezeit z genügt, weiterhin, wie aus den 
Meßwerten einer Stichprobe die Parameter dieser Ver- 
teilung bestimmt werden können, und schließlich, wie 
ein Prüftest zum Vergleich zweier Stichproben etwa 
zur Beurteilung zweier derartiger Schmelzanlagen 
aufgebaut werden kann. 


Z. angew. Math. Mech. 
Bd. 37 Nr. 9/10 Sept./Okt.1957 


Y 


a 


BIETET AN 


 Z. angew. Math. Mech. 


Bd. 37 Nr. 9/10 Sept./Okt.1957 


Kleine Mitteilungen 405 


1. Bestimmung der Verteilung w(z) des Differenz- 
betrages z zweier unabhängiger Zufallsvariablen x und 
y mit gleicher Normalverteilung. 

Die Wahrscheinlichkeit dafür, daß das Merkmal z in 
das Intervall z, bis z, fällt, ist gegeben durch das Inte- 
grel über w(z) von z, bis z,. Andererseits ist über der 
(2, Y)-Ebene der stochastisch unabhängigen Merkmale x 
und y durch das Produkt u(x) - v(y) der Einzelverteilun- 
gen eine zweidimensionale Verteilung bestimmt, deren 
Anteile innerhalb der durch z, = jR — y,und 2 = 
e —y| gegebenen Gebiete gerade die gleiche Wahr- 
scheinlichkeit liefern. 

Im vorliegenden Falle sind die Gebiete zwei Streifen 
S, und $, zwischen den Geraden 


Yen ME Ep 
yanzı für y>x und Me, für y<x 
Demnach gilt 
EX 2a 
S w(z) de = Sf u(x) v(y) de dy. 
2 Sı4S3 


Führen wir zur Umformung dieser Integrale neue Inte- 
grationsvariable ein, in $,: 


2=y—x 

Ge 
und entsprechend in $,: 

= —y 

BE 


so erhalten wir 


Feo)da—f S wevletz) dt det FF we) vie) Ede. 
Zı ee 


A —0 


und schließlich 
we) = Su ltd Hr), 2>0, 


wobei wir angenommen haben, daß die Verteilungs- 
funktionen so beschaffen sind, daß die Integrale 
existieren und die Vertauschungen erlaubt sind. Die 
Auswertung ergibt für die Verteilungen (1) 
ven 
A ee ea ei 

w(z) = 2 (2) 
Wie man erkennt, ist diese Verteilung völlig unsymme- 
trisch; sie besitzt die Gestalt der rechten Hälfte einer 
Gaußkurve und ist durch den Parameter h bereits 
bestimmt, wobei sich die Frage erhebt, wie dieser 
durch die Meßwerte einer Stichprobe festgelegt werden 
kann. 

2. Bestimmung des Parameters in der Wartezeit- 
verteilung aus einer Stichprobe nach der Methode der 
maximalen Gesamtwahrscheinlichkeitsdichte. 

Für die Berechnung des Parameters aus Stichproben- 
werten fehlt uns für die Verteilung (2) ein Zusammen- 
hang, der den Formeln für Mittelwert und Streuung 
bei einer Normalverteilung entspricht. Da diese die 
Bedingungen für eine maximale Gesamtwahrschein- 
lichkeit eines Wurfes von Meßwerten darstellen und 
eben auf diese Weise als deren charakteristische Maß- 
zahlen Eingang gefunden haben, beschreiten wir für 
unsere Aufgabe den gleichen Weg. j 

Vorgegeben seien n Meßwerte 2,, 23,...,2n. Die 
Wahrscheinlichkeitsdichte, daß der Meßwert z, als 
auch der Meßwert z, usw. in die sie umgebenden Inter- 
valle fallen, ist gegeben durch 


w(z) - w(2,) . . . Win) = = hre 


die in h ihr Maximum erreicht für 


Mit (3) haben wir eine Vorschrift zur Berechnung des 
Parameters h gefunden und können nun der Frage 
nachgehen, welcher Verteilung diese Variablenver- 
bindung genügt. 


3. Verteilungsfunktion der Variablen 3 — 3 Er 
y al 
mit stochastisch unabhängigen Variablen z, der Ver- 
teilung (2). 

Für die Berechnung der Verteilung dieser Größe 
ist zu beachten, daß der Form nach für die z, zwar eine 
Gaußverteilung vorliegt, daß diese jedoch gegenüber 
dem üblichen Fall auf positive Argumente beschränkt 
und dafür mit dem Faktor 2 überhöht ist. 

Für n = 2 betrachten wir 

3=-iti: 10a; 


und bilden in bekannter Weise 


3 , 

S W;(s,) ds, = S S w(z,) w(2,) dz, dz, 
sL G 

2 


über ein Gebiet @, das aus einem durch s und sy 
bestimmten Viertelringstreifen besteht. Führen wir 
neue Variablen ein durch 


21 = 5,0089 
2, = $,8nY, 
so ergibt sich schließlich 


2 


W,(s,) = s,; S W (s, 608g) w(s, sin o) dp 
Ö 


h? 
En) Sa aa 
a 


d. i. für A? s$ eine y?-Verteilung der Ordnung 2. 
Durch vollständige Induktion finden wir für 


=-1t2+4: +2 
die Verteilung 


Wal) = 


v—ı 
>| 2 ): h\n _— 
Van)! (2) Bu 


und stellen damit fest, daß h? 3 z? einer y?-Verteilung 
ee 


mit dem Freiheitsgrad n genüst. 

4. Prüftest zur Beurteilung zweier Stichproben von 
Wartezeiten. 

Gegeben seien zwei Stichproben z und Z vom Um- 
fang n, und n,. Es soll beurteilt werden, ob beide aus 
einer Grundgesamtheit stammen können, oder ob 
zwischen beiden wesentliche Unterschiede bestehen. 
Wir bilden 


ho m >; N, = 
Fi 2,2; N. 22% 
€ EN Ipy= ri 
F= Ha NE >12...) 
Fre 
Mn ne 


und erkennen, daß gerade dabei der unbekannte Para- 
meter herausfällt. Die Prüfung der Stichproben auf 
ihre zufällige oder wesentliche Abweichung läßt sich 
demnach in üblicher Weise mit der F-Verteilung mit 
den Freiheitsgraden (n,; n,) durchführen. 

Nach Erledigung der Fragen für die Wartezeiten 
sollen noch die für die Anschlußzeiten betrachtet 
werden. 

In unserem Beispiel der beiden Vorschmelzöfen soll 
uns nun nicht deren Wartezeit interessieren, sondern 
die Anschlußzeit i für den Arbeitsbeginn des dritten 
Ofens. 

Wir erkennen, daß 


i=y ist;weın yz2%, 


t=% 1:6, wenns Yo. 
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Bezeichnen wir die gesuchte Verteilung von t mit w(t), 
so gilt entsprechend den früheren Überlegungen 


tz 

S wit) d = N) v(y) da dy , 

Hl 

wobei hier das Gebiet B ein rechtwinkliger Streifen ist, 
bestimmt durch 


Mey: wenn yz% 


= Yp» 
und 

h=%, b=%, «wenn ER 
ist. 
Wir können das Doppelintegral zerlegen 


iz 


Foto) dt = foty) | Sula) ie) dy 


th h 
t, tz 
+ fl) (Lot a) 


woraus durch Differentiation nach t, mit , >t,=1 
bei Voraussetzung stetiger Integranden folgt 


t t 
w(t) = vit) [ ulx) de + ult) S vly) dy 
) —oo 
und für den Sonderfall w(t) = v(t) 


w(tl) = 2 ul) f ua) dx. 


Wählen wir «(t) wie in (1) als normal verteilt, so gilt 
21 un IR 
wit) = — ee "it. — | a 
Vr Vr 
—@ 


Zur Charakterisierung dieser interessanten Verteilungs- 
funktion bilden wir die ersten Momente 


© 
; 1 
(t) h(t— u) d = —— 
wid 
—6,0) 
[0,6] 
i 1 
| w(t) h? (t — u)? di = > 
ee 


bzw. 
[0.6] 


i 1 
| w(t) \ Bar er )a —=( 


[0 0) 


| it) \t 2 I 1 
’ a Il = nn 
wet) Be Ver dt a 


—00 

womit der Mittelwert und die Streuung der Verteilung 

w(t) bestimmt ist: 

l 1 
Vz Y2h 
ar—1l 1 


a ro 
7 23h? 


t=u- 
. (6). 


Wir entnehmen daraus, daß die Streuung von wit) 
geringer. ist als die von (t), und daß der Mittelwert 
rechts der Maximumabzsisse u liegt. Die Stelle t, 
des Maximus von (ft) ergibt sich mit 


1 a 1 2 
Ii— u) =2 oa)= et; De) ns e "dr 
Ir \z. 
aus der Gleichung is 
Fl) 
a 


so daß wir mit 2m = 0,358 finden 
1 
zu er 
im u + 0,5 V2% 


Die Parameter der Verteilung lassen sich aus einer 
Stichprobe im Prinzip wieder mit Hilfe des Verfahrens 
der maximalen Wahrscheinlichkeitsdichte berechnen, 
jedoch liefern die hierbei auftretenden transzendenten 
Bestimmungsgleichungen die gesuchten Parameter 
nicht explizit. 

Eine Möglichkeit, sie durch einfache Formeln fest- 
zulegen, bieten die Beziehungen (6), aus denen. wir u 
und h ermitteln können, wenn dabei t und o? durch die 
Stichprobennäherungen ersetzt werden. Die Vertei- 
lungen dieser Größen für die Durchführung eines Prüf- 
testes streng zu ermitteln, ist bei der vorliegenden 
nichtnormalen Grundverteilung außerordentlich schwie- 
rig. Es ist daher vorteilhaft, für die Verteilung 


t 
w(t) = 2 uft) S ulx) de 

—o 

die Summenfunktion zu bilden 


Fett) de— | Fuß) au) 


und daraus 


Kae de = | u dr (7). 


—c&9 


Da (x) eine Normalverteilung ist, findet man für 
praktische Zwecke völlig ausreichend die Parameter 
aus der Auftragung auf Wahrscheinlichkeitspapier. Die 
Zufallssrenzen für die Lage der ausgleichenden Geraden 
zu bestimmen, ist allerdings nur näherungsweise mög- 
lich. Die Frage eines regulären Prüftestes für zwei 
Stichproben ist demgemäß für die Anschlußzeiten viel 
schwieriger zu beantworten als für die Wartezeiten, so 
daß man in diesem Falle einen Rangtest heranziehen 
wird. 
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An Analogy Between the Flexural Vibrations 
of a Cone and a Disc of Linearly Varying 
Thickness. 


Summary 


The following investigation concerns an analogy bet- 
ween the lateral vibration of a conical bar and the 
axisymmetrical vibration of a circular disc of linearly 
varying thickness. The vibration of a conical bar has 
been investigated very completely both as regards the 
caleulation of the periods and the nodal positions, and 
it follows that the results are applicable to a disc having 
the same boundary conditions. 


Analysis 
‚ Consider the conical bar generated by rotating the 
line y—= A x about the x-axis. From the Bernoulli- 
Euler equation of bending, the differential equation 
governing the vibration of the bar is 


ern. Pa 
zie n ze) + re 2 m 
or 
02 R ow 40 ow 
2ie 3) Er at nd 


where w is the displacement and o the density of the 
material. 
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VE 


D=D,® where D, = Eh}j12 (1 —v*) 
after some simplification 


REN et we) 
is the slope of the middle surface and h is the 
ss at radius x. As might well be imagined, so- 


gi ons to this equation are generally diffieult to ob- 
 tain [2]. However, for a disc whose thickness varies 


linearly with the distance from the center we write 
and obtain, 


7] w w ow 
> ar 


Bel) 0% 


Taking Poisson’s ratio v as 1/3, which is a good value 
for several materials, the above equation may be 
written as 

@l[ Mu 32o „dw 

a ae) tat er! 


. (4) 


which is identical with equation (1) if o in the latter 


equation is replaced by 8 A? 0/3 hd. Thus the problems 
of the vibration of the cone and dise of linearly varying 


 thickness are analogous if the boundary conditions 


* 


Be Sl u aa Sin na ln zn li a a nn 


Na \ 


have the same form in each case. This is so in the case 
of clamping, where the boundary conditions are w = 
ow/öx = 0 for both the cone and plate but it is not 
generally the case. 

The vibration of a complete cone of length a clamped 
at the base has been investigated in a very complete 
manner by Wrinch [3], who showed that the trans- 
cendental equation for the periods is 


J,(a) I,(«) = Jı(a) I;(«) 


To Fa. air RER) u 
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and writing h,a = H, the outer thiekness of the dise, 
we have a : 


2m _ 32m af )/2 (5.906, 9.197, 12.402) 


> 


v H 3E 
32a 1/20 04 Ex 
= 1/2 (34.88, 84.59, 153.81)-1 
H Vs m 59, 153.81) 


The corresponding solution for the solid disc of constant 
thickness H is well known, the periods for vibration 
in the first three axisymmetrical modes being written 
for comparison purposes in the form 

2n  32natı/20 =; 
: z SE (40.84, 159.12, 355.60) (8) 
It follows that the ratios of the periods of the variable 
to the constant thickness disc in the first three modes 
are 1.17, 1.88 and 2.31 respectively. The following state- 
ment, repeated by analogy from the investigation of 
Wrinch [3], may also be made: In the axisymmetrical 
vibrations of the above disc, the nodes in the nth 
mode occur at intervals of order 1/n in the neighbour- 
hood of the outer edge, and at intervals of order 1/n? in 
the neighbourhood of the center. 
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tisch bestimmter und unbestimmter Stabtragwerke. 
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Der hier wiedergegebene Vortrag, gehalten vom Verfasser im Februar 1956 in der Sektion Maschinenbau 
der Deutschen Akademie der Wissenschaften zu Berlin, zeigt die Aufgaben des Strömungsingenieurs im 
Maschinenbau sowie die Möglichkeiten ihrer Lösung durch mathematische oder experimentelle Methoden 
an aaglellei: Beispielen. Die Beispiele sind aus den Arbeiten des Institutes für angewandte Strömungs- 
lehre der Techxiischen Hochschule Dresden gewählt. 


Eine enge Zusammenarbeit zwischen Ingenieur und Angewandtem Mathematiker wird empfohlen. 


Nach einer gegebenen Übersicht über neuentwickelte Methoden zur Untersuchung und Gestaltung von 
Flügelgittern werden einige Probleme der Staubabscheidung, des Hydrozyklones der Diffusorauslegung 


des Einsatzes von Strahlpumpen besprochen. 
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Das Buch ist einem der interessantesten und modernsten Gebiete der Geometrie gewidmet. Der Autor be- 
handelt die Flächenverbiegung großenteils mit den Methoden der Differentialgeometrie und bringt sie mit 
der Theorie der Polyeder in Verbindung. Untersuchungen über die Verbiegbarkeit geschlossener und offener 
Flächen mit den mannigfachsten Regularitäts- und Randbedingungen werden durchgeführt. Im Mittel- 
punkt des Interesses stehen die eindeutige Bestimmtheit von Flächen durch die Metrik, stetige Verbiegsam- 
keit im Großen, die infinitesimale Verbiegbarkeit erster Stüfe. Ferner sind die merkwürdigen Eindeutig- 


keitssätze bei Flächenstücken mit einem Flachpunkt berücksichtigt. Eine Zusammenstellung-in deutscher 


Sprache fehlte bisher. 
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